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ABS TRACT 


In  this  work  the  coaxial  waveguide  antenna  is  treated  by  the 
Wiener-Hopf  technique  and  the  transient  radiation  from  a  cylindrical 
monopole  is   developed  in  the   light  of  the  rigorous   results  obtained 
from  the  Wiener-Hopf  analysis.     Analytic  expressions  are  derived  for 

(1)  the  electromagnetic  fields  in  the  feed  line  and  (2)   the  far  zone 
radiation  fields  of  the  coaxial  waveguide  antenna,  with  time  harmonic 
excitation  voltage.     Complete  characterization  of  the  transient  beha- 
vior is   also  found  for  (1)   the  fields  interior  to  the  feed  line  and 

(2)  the  radiated  fields   for  excitation  voltages  arbitrary  in  their 
time  dependence  to  the  extent  that     kb,ka  <<  1    .     This  corresponds 
to  the  case  of  a  thin  antenna  and  excitation  voltage  with  a  non-zero 
rise  time,  specifically  chosen  so  that  frequencies  violating  the 
restriction     kb,ka  «  1     are  negligible. 

The  transient  radiation  from  the  cylindrical   monopole  is 
developed  in  a  closed  analytic  form  which  is   relatively  easy  to 
interpret  and  apply.     The  expressions  found  offer  an  alternative  to 
transient  analysis  by  conventional  methods   requiring  numerical    tech- 
niques involving  extensive  computer  calculations.     They  are  also  the 
basis  for  an  uncomplicated  procedure  to  synthesize  a  desired  behavior 
of  the  transient  radiation   from  cylindrical  monopole  antennas. 
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1.      INTRODUCTION 

This   dissertation  presents  an  analysis   of  the  coaxial  waveguide 
antenna  by  the  Wiener-Hopf  technique  and  the  determination  of  transient 
radiation  from  the  cylindrical  monopole  antenna.     This  work  is  the 
first  successful   treatment  of  a  semi -in  finite  coaxial  structure  by  the 
method  of  Wiener-Hopf  and  the  first  rigorous  determination  of  transient 
radiation  emanating  from  an  open  ended  semi-infinite  waveguiding  struc- 
ture.    In  addition,  a  new  and  simple  method  for  calculating  the 
transient  electromagnetic  radiation  of  a  cylindrical  monopole  antenna 
is  presented. 

The  study  of  electromagnetic  boundary- value  problems  by  the 
Wiener-Hopf  technique  began  in  the  early  1940's.     Among  the  most  notable 
of  the  early  workers  were  Schwinger  [1]  and  Copson  [2],     The  early 
applications  of  this  powerful  method  were  inthe  rigorous  solution  for 
the  fields  within  a  bifurcated  parallel  plate  waveguide  and  diffraction 
of  plane  waves  by  a  conducting  half-plane.     The  method  was   later  used 
to  solve  the  problem  of  radiation  impinging  on  or  emerging  from  the 
open  end  of  an  unflanged  circular  pipe  and  a  pair  of  parallel  plates. 
Much  of  the  early  history  of  application  of  the  Weiner-Hopf  method  is 
detailed  by  Bouwkamp  [3]  and  an  authoritative  account  of  the  method  was 
published  in   1958  by  Noble  [4].     In  this  work  the  integral   equation 
approach  of  the  early  workers    is    not  stressed  .       A   more  direct  simple 
method  due  to  Jones  [5]    is    emphasized.     It  is  the  application  of 
Wiener-Hopf  technique  due  to  Jones  that  will  be  employed  within  this 
report. 
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There  have  been  several   recent  investigations  of  the  transient 

behavior  of  antennas   [6-15].     The  focus  of  many  of  these  investigations 
is  on  the  time  behavior  of  the  radiated  electromagnetic  field.     Refer- 
ences  [6,8]  and  [14]  present  theoretical   as  well   as  experimental 
results.     In  most  of  the  prior  research  efforts,  the  problem  is   formu- 
lated in  terms  of  circuit  theory  and  the  antenna  is   treated  as   a  lumped 
circuit  element.     The  feedline  and  the  exciting  sources  are  replaced  by 
their  Thevenin  or  Norton  equivalents.     The  formulation  of  the  problem 
in  terms  of  circuit  theory  yields  complex  frequency  domain  fields  which 
must  be  Fourier  transformed  if  the  time  behavior  is  of  interest.     Such 
calculations   require  numerical  methods   due  to  their  extreme  complexity. 
The  work  of  Schmitt,  Harrison  and  Williams   [6]  is  a  determina- 
tion of  the  transient  radiation  and  reception  performance  of  a  thin, 
finite   length  cylindrical  monopole  over  a  perfectly  conducting  infinite 
ground  plane  with  a  coaxial    feed.     In  the  circuit  analysis   formulation 
(frequency  domain)   the  antenna  was   represented  by  an  impedance  and  the 
radiated  field  related  to  the  voltage  across  the  antenna  impedance 
through  a  proportionality  factor,   called  effective  height.     The  time 
dependent  behavior  is  obtained  by  Fourier  transforming  the  frequency 
domain  field.     A  similar  procedure  was  used  to  obtain  the  properties  of 
the  antenna  in  reception;   the  induced  voltage  on  the  receiving  antenna 
was   related  to  the  incident  electric  field  intensity  through  the  effec- 
tive height.     The  theoretical  work  of  Schmitt,  Harrison  and  Williams 
was  supplemented  in   1971  by  Abo-Zena  and  Beam  [9]  to  include  the  time 
behavior  of  near  zone  field  intensities  and  in  1970  by  Palciauskas  and 
Beam  [10]  to  include  the  far  zone  time  behavior  at  latitudinal   angles 
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other  than  90°,  which  was  the  only  angle  considered  by  Schmitt  et  al . 

King  and  Harrison  [8]  formulated  a  frequency  averaged  reflec- 
tion coefficient  to  theoretically  obtain  the  reflected  electromagnetic 
fields  on  the  coaxial   feed  line  of  an  infinitely  long  cylindrical 
antenna  or  for  the  time  interval  before  the  excitation  pulse  can 
reach  the  tip  of  a  finite  length  antenna.     Use  of  frequency  averaging 
is  necessitated  by  the  fact  that  the  antenna  input  admittance  is 
frequency  dependent,  which  constrains  the  conventional    transmission 
line  reflection  coefficient  to  be  frequency  dependent. 

The  transient  current  on  a  circular  tubular  infinitely  long 
antenna  that  is  excited  by  a  voltage  which  is  a  step  function  of 
time  was  determined  in   1961  by  T.  Wu.     He  later  (1969)   furnished  a 
corrected  analysis  for  this  current  in  Reference  [11].     An  earlier 
work  (1960)  by  Brundell   [12]  treated  the  more  general  problem  of 
determining  the  space  and  time  behavior  of  the  electromagnetic  field 
for  such  an  antenna.     This  work  was  also  supplemented  in  1970  by 
Latham  and  Lee  [13]  who  found  the  early  and  late  time  behavior  of 
the     field  radiation  from  a  hollow  infinite  cylindrical   antenna 
excited  by  a  voltage  that  is  a  step  function  of  time. 

Lamensdorf  [14]  determined  that  the  time  behavior  of  the 
coaxial  cone  antenna  in  reception  is  primarily  a  time  derivative  of 
the  incident  electric  field  intensity  falling  on  the  aperture.   This 
work  also  contains  a  heuristic  discussion  of  the  time  dependent  elec- 
tric field  radiated  by  a  cone  antenna.   Chang  [15]  analyzed  the  tran- 
sient reception  characteristics  of  the  annular  slot  antenna  using  the 
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circuit  theory  approach.     References   [14]  and  [15]  are  recent  works, 

1970  and  1971   respectively,  and  have  application  as  part  of  a  receiv- 
ing system  for  measuring  scattered  impulse  responses  of  various 
objects  or  as   a  receiving  antenna  on  an  aircraft  or  guided  missile. 

In  t  he  systematic  study  of  radiation  from  the  coaxial   aperture 
antenna  which  follows,  analytic  expressions   for  the  electromagnetic 
fields  of  the  coaxial   aperture  antenna  are  given.     These  expressions 
account  for  waves  traveling  in  both  the  forward  and  reverse  directions 
on  the  feedline  and  the  radiation  fields   for  harmonic  time  dependent 
excitation.     Complete  characterization  of  the  transient  behavior  is 
found  for  the  fields  interior  to  the  feedline  and  the  radiated  fields, 
for  excitation  voltages  which  are  arbitrary  to  the  extent  that 
kb.ka  «  1    .     This   corresponds  to  the  case  of  a  thin  antenna  and 
excitation   voltages  with  non-zero  rise  times,  specifically  chosen  so 
that  frequencies  which  violate  the  restriction     kb,ka  «  1     are  neg- 
ligible. 

The  annular  slot  antenna  is  considered  to  be  a  special   case  of 
a  cylindrical  monopole  antenna  (h  =  0).     The  fields   radiated  from  the 
annular  slot  antenna  are  approximated  by  well  known  methods  in  which 
the  incident  field  on  the  aperture  is  considered  to  be  the  aperture 
field.     The  radiated  fields  so  determined  are  interpreted  in  light  of 
the  rigorous   results  obtained  from  the  coaxial  waveguide  antenna  and 
a  model  is  obtained  for  a  cylindrical  monopole  antenna  in  which  the 
height  of  the  radiating  element  is  greater  than  zero  (h  >  0).     The 
transient  radiation  field  intensities  are  obtained  from  this  antenna 
model.     The  determination  of  the  transient  radiation  through  the 
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procedures  developed  in  this   research  offer  an  alternative  to  comput- 
ing transient  behavior  by  numerical  methods  where  extensive  machine 
calculations  are  required.     The  procedures   reported  on  herein  can  also 
be  used  as  the  basis  of  an  uncomplicated  synthesis  procedure  to  produce 
a  desired  time  behavior  in  fields   radiated  from  the  coaxial  waveguide 
and  the  cylindrical  monopole  antennas. 

The  results  of  this   research  should  prove  to  be  \/ery  useful  in 
the  areas  of  1)  determination  of  antenna  properties   (driving  point 
impedance  and  radiation  pattern)  through  time  domain  measurements  [16], 
2)   identification  of  and  discrimination  between  conducting  bodies 
through  radiation  scattering  from  their  surfaces  and  3)   the  calibration 
of  equipment  designed  to  measure  the  effectiveness  of  shields  designed 
to  protect  equipment  or  installations  from  intense,  short  rise  time 
electromagnetic  pulses.     In  each  of  the  above  areas  a  precise  knowledge 
of  the  time-behavior  of  the  radiated  electromagnetic  field  is  required. 
Additionally  the  ability  to  synthesize  a  particular  time  behavior  in 
the  radiated  field  is  highly  desirable  and  yery  important. 
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2.     THE  COAXIAL  WAVEGUIDE:      AN  APPLICATION  OF  THE  WIENER-HOPF  TECHNIQUE 

In  this  chapter  the  frequency   (spatial)   domain  electromagnetic  in- 
tensities for  an  open  ended  coaxial  waveguide  will   be  derived.   This  will 
be  done  using  the  methods  of  Wiener  and  Hopf  [4],     Use  of  this  technique 
will  provide  a  determination  of  transverse  electromagnetic  mode  (TEM) 
as  well   as  higher  order  mode  intensities  within  the  structure. 

We  take  the  coaxial  waveguide  structure  to  be  that  of  Figure  2-1. 
The  inner  conductor  is   a  hollow  tube.     Both  outer  and  inner  conductors 
are  assumed  to  be  semi-infinite  in   length,  infinitesimally  thin  and 
perfectly  conducting.     The  cross  section  is  circular  providing  for  a 
coaxial   annular  region  between  the  conductors.     We  also  assume  that  the 
waveguide  is  embedded  in  a  homogeneous,  isotropic,  dielectric  medium 
having  the  constitutive  parameters  of  free  space. 

Specification  of  the  Problem 

We  begin  the  specification  of  the  problem  by  listing  the  Maxwell 
equations: 

n   „  t  & 


-*- 


V*H     -     !♦§» 

V  •    D     =     p 

V  •   B     =     0  (2.1) 

The  fields  which  satisfy  the  above  equations  must  also  satisfy  the  con- 
stitutive relations: 
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(rA+) 


Figure  Z-l.     The  coaxial  waveguide  antenna 
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B     =     uQH  ,  D     -     et 


where 

J 
c 


2  =  Vo 


Under  the  assumption  that  the  coaxial  waveguide  has  been  excited  in 
TEM  mode,  there  is   a  traveling  electromagnetic  disturbance  proceeding 
in  the  positive    z    direction.     The  time  origin  is   taken  to  be  the 
point  in  time  when  the   leading  edge  of  the  TEM  mode  exciting  field  is 
incident  upon  the  aperture   located  in  the     z  =   0     plane. 

The  specific  method  used  to  initially  excite  the  coaxial   struc- 
ture and  the  field  distribution  prior  to  the  time  origin  are  not  of 
interest  in  this  work. 

It  is  well   known  that  the  TEM  mode  field  of  a  coaxial  waveguide 
is  independent  of  angular  variation  [17,  p. 326]  in  the  cylindrical 
coordinate  system.     The  coaxial  structure  also  does  not  introduce  any 
boundary  conditions  which  have  angular  dependence.     From  these  observa- 
tions we  conclude  that  the  scattered  field  must  also  be  characterized 
by  the  same  symmetry  (i.e.,  independent  of  angular  variation  in  cylin- 
drical  coordinate  system). 

In  cylindrical   coordinates  the  field  components   can  be  written 


as 


and 


E     =     E  a     +  E  ,a.  +  E  a 
p  p        $  4        z  z 


H     =     H  a     +  H  .a.   +  Ha 
P  P         <J>  <J>         z  z 


where     a   ,a,,a       is  the  triad  of  unit  vectors  along  the  coordinate  axes 

P     4>    z 
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We  now  Fourier  transform  the  Maxwell  equations  to  change  from  a 

time  domain     to  a    frequency  domain  representation.     The  Fourier  trans- 
form pair  [18]  is  defined  as   follows: 


F(co)  = 


f(t)  = 


1 


fin 

1 

fin"   J 


oo 

r 


f(t)  eia)t  dt 


F(co)   e  doj 


When  the  first  two  Maxwell  equations  are  multiplied  by         -e11 


operated  on  by  the  integral  operator,  the  results  are: 


fin 


1 

fin 


V  x  E  eiu)t  dt 


1        (  3B     ioot    .. 


and 


f  V  x  H  eiwt  dt 


fin 


fin 


rt  j.  3Dt  „iwt    .. 
[J  +  H]  e        dt 


Interchanging  the  curl   operator  and  the  integral  sign  on  the  left  hand 
side  and  integrating  by  parts  on  the  right  hand  side  yields  equations 
in  the  frequency  domain.     To  form  these  results,  use  has  been  made  of 
the  fact  that  the  field  intensities  must  be  zero  as     |t|  ■+  °°  . 

V  x  E(co)     =  +  iu  B(w) 

V  x  H(cj)     =     j(o))   -  iw  D(w) 


It  similarly  follows  that 
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V  •  D(tu)     =     p(u) 

V  •    B(u>)      =     0 


B(w)   =  u0H(w)        and       D(w)   =     £QE(a3) 


When  the  field  components  in  cylindrical  coordinates  are  substituted  into 
the  frequency  domain  Maxwell  equations   for  a  source  free  region   (i.e., 

-y 

J(oj)     is  set  to  zero),  we  obtain  the  following  relations  between  the 
field  components: 

i.Bp  =   -^  (2.1.1) 

l.i,.    -t*-^  (2.1.2) 

9E  9E7 

and 

-iu>  D     =   -  |-H.  (2.1.4) 

p  92      <p 

HA       3H. 
-iu  D    -     -i+  J  (2.1.5) 

Z  P  dp 

9H         9H 

-"•V  ^-jf  (2-'-6) 

By  using  the  constitutive   relations  and  combining  equations  (2.1.3), 
(2.1 .4)  and  (2.1 .5),  we  derive  the  scalar  wave  equation  for    H     , 


This  equation,   upon  simplification,  becomes: 
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3p  p        p  3z  c  r 

where     ye     =  -y    .       Si  mi larly  equations   (2.1.1),   (2.1.2)   and  (2.1.6) 

c 
produce 

[72^-V  ?7  +  7]E*(u)  =  °  (2-3) 

3p  M  p  9Z  C  r 

From  the  foregoing  equations  it  is  apparent  that  the  solution  of 
(2.2)   and  (2.3)   for    H     and     E     is  enough  to  completely  determine  all   field 

components.     The  work  incorporated  in  Appendix  A  shows  that     E      must 
be  zero  for  the  coaxial  waveguide,  therefore  the  only  field  intensity 
components  present  are     H   ,  E     and     E     •       The  rest  are  identically 
zero.     The  analysis  has  been  reduced  to  solving  the  boundary  value 
problem 

2  2 

[3       t  1_  +  1    *        1     +  k2]  H   (p.z.co)     =     0  (2.4) 

3ZZ        8pZ       p   9p       pZ  * 

H  ,(0,z,u))     =     0    ;  -»  <  z  <  od  (2.4.1) 

E   (a,z,w)     =     E   (b.z.w)   -     0   ;     z<0  (2.4.2) 


The  Sommerfeld  radiation  condition, 


9H 
lim     r(T-2-  -  ik  H.)  =     0     must  also  be  satisfied. 


At  this  point  it  is  useful   to  point  out  that  in  the  above  equations 

kZ  =  \  (2.4.3) 

c 
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is  a  scalar  quantity  and  equation    (2.1.5)   can  be  written  as 

This  will  be  of  use  later. 

The  well  known  expressions   for  the  TEM  mode  fields  for  the 
coaxial  waveguide  [17]      will   be     given      here  to  set  the  notation.     We 
take  the  field  intensities  incident  on  the  aperture  from  the  negative 
z  direction  to  be 

Ei     =  jd«L_    e^L 

p        inb/a        p 

$         \l\s     £n  b/a      p  p  ' 

v(u>)     is  the  frequency  domain  voltage  between  the  inner  and  outer  con- 
ductors. 

The  space  containing  the  waveguide  is  divided  into  three   regions, 
portrayed  in   Figure  2-2.   Equation  (2.4)  applies  to  each  of  these  regions. 
The  field  satisfying  this  equation  in  region  one  is  designated 
H     (p,z,to)    .     Similarly  we  have     H,    (p,z,w)     for  region  two  and 

<J>1  q>2 


H,   (p,z,o))     for  region  three.     The  equation  for  each  of  the  regions  is 

<j)3 

a  partial   differential  equation  in  the  variables     p     and     z  .     By 
Fourier  transforming  in  the  spatial   dimension     z   ,  an  ordinary  differen- 
tial equation  in  the  variable     p     is  obtained. 

Following  the  notation  of    [4]   ,      the  spatial    Fourier  transform 
pair  is   defined  as   follows: 
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wmwmmmutiw/j*. 


2b 


ZMW?/M7/y7AP/fo7/ZW7, 


Region  I 


'//////////////////////////. 


Region   1  b  <  p 

Region  2  a  <  p  <_  b 

Region   3  p  <  a 


Figure  2-2.     Division  of  the  space  surrounding  the  coaxial  waveguide 
antenna  into  regions 


F(o) 
f(z) 


1 


fin     J 


1 

/2~7    J 
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f(z)  eiaz  dz 


TaZ 


F(a)  e  da 


(2.8) 


(2.8.1) 


The  "plus"  and  "minus"  functions   are  defined  as 


F+(a)=     ' 


•Sir  i 


f(z)  e1az  dz 


(2.8.2) 


where 


F"(a)  =  -J-        f(z)  elaz  < 


F(a)   =  F  (a)  +  F'(a) 


(2.8.3) 


(2.8.4) 


a     =  a  +  ix         is  complex  . 

Using  well   known  theorems  involving  Fourier  transforms  in  the 
complex  plane  [4],   the   regions  of  analyticity  of     F  (a)     and     F~(a) 
will  be  found.     Suppose  that  we  are  given     x     <  x  <  x+  ,  all    real 
numbers.     In  addition  suppose  that 

|f(z)|     <    A  exp(x_z)        as  z  ■*  °° 


where     A     is  a  constant.     Then 

M 


|F+(a)|    < 


1 


>/2¥ 


f(z)  elaz  dz 


0 


v^fT 


A  e 


( i  a  +  x   )  Z 


dz 


where     M  <  °°     is   chosen  large  enough  for  the  inequality  to  be  valid. 
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The  last  integral   of  the  above  expression  converges  if  and 


only  if 


(ia+  t  )z 

e  ■*  0 


as     z  •* 


or 


(ia  -  x  +x_)z 

e  ■*  0 


as     z  -*■  °° 


This  last  expression  is   valid  provided 


T    >    T 

By  similar  reasoning,  it  can  be  shown  that  if     |f(z)|   <  B  exp(x  z) 
as     z  -*■  -»  (B  constant)     then     |F~(a)|      is  convergent  provided 

T    <    T+ 

Therefore     F  (a)     and     F~(a)     are  convergent  in  the  strip 
t_  <  x  <  x+     of  the  a  plane. 

We  also  note  that  if     |F  (a)|     and     |F~(a)|      are  convergent 
in  a  region  of  the  complex  plane,  then  the  magnitude  of  nth  ordered 
derivative   (n  <  °°)   of     F  (a)     and     F~(a)     with   respect  to    a     also 
exists  in  that  region.     In  fact  such  functions  are  nth  order  continu- 
ously differentiable  in  the  region. 


^(F+(a)) 
da 


da 


Fr(a)|   e 


i<Ha) 


UlL 


SB 


znf(z)e1az  dz 


0 


.n 


M 


v^rf 


znf(z)  elaz  dz 


♦  -L  I  z"  . 


(ia  +  x   )z 


dz 
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where     M  <  °°     is  properly  chosen  to  ensure  the  inequality.     Concen- 
trating on  the  term  of  significance,  we  write 


n     (ia  +  TJz  n 

zn  e  dz  =  zn  e 


(ia  +  T   )z 


M 


M 


la  +  t 


n      (ia+  T    )z 
zn_l   e  dz 


k 


After    n     integrations  by  parts,   the  above  becomes 

?     n!   (z)P(-1)n-p     j1atTJ2 
P*!1*     (ia+T)n"P 


M 


„  „/„   i\        i      (ia  +  x   )z 
+  (_i)n  niDdLLLLii  e 
/ .  %n 

(ia  +x   ) 


which  is  finite  for  finite     n     provided 


T    >    T 

Again  by  a  similar  discussion,  it  can  be  shown  that     F"(a)     is  nth 
ordered  continuously  differentiate  in  the  open  region     t  <  t+     provided 
|f(z)|   <  BeT+z     as     z  +  -»  . 

Since  all   that  is   required  for  analyticity  in  a  region  of  the 
complex  plane  is   that  the  first  derivative  of  the  function  exist  in  that 
region,  we  are  assured  that     F+(a)     and     F~(a)     are  analytic  in  the 
strip     x     <  t  <  t+     of  the  a  plane  provided  of  course  that 

t  z 


|f(z)|    <  A  e 


and 


t,z 


|f(z)|    <  B  e  + 


Upon  multiplication  of  equation  2.4  by    e  a      and  operat- 

ing  with  the  integral         dz  ,   the    following  equation  in 

—CO 

p     emerges. 


»£?    J   dz 
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2  2 

-yH^p.z.oOe        +  [-T+  -^ 


9p 


^y+  k2]  H^p.a.u) 


=  0 


Integrate  ng  by  parts  twice  and  imposing  the  conditions  that  H  (p,z,oj) 
and  -j—  H  ,(p,z,co)  -+ 0  as  jz|  ■*■  °°  ,  produces  an  ordinary  differential 
equation  in     p   , 


dp         K     K       p  T 


(2.8.5) 


2         2.2 

Y     =  a     -  k 


(2.8.6) 


It  is  noted  at  this  point  that  the  incident  magnetic  intensity  does 
not  satisfy  the  differential  equation  (2.8.5)  since  it  does  not  converge 
in  the  limit  as     z  -*•  -°°  .     This  will  be  explored  in  detail  a  little 
later. 

For  region  1: 


r— 

a  2 

dp 


P  dp 


Y  ]  H,    (p,a,u>)     =     0 


"♦ 


(2.9) 


H,    (p,a,co)  ■*  0  as       p  -> 


'♦ 


1 


(2.9.1) 


-r-  [p  H"   (p.a.u))] 


dp 


'♦ 


=     0 


p=b 


(2.9.2) 


Equation  (2.9.1)   follows   from  the  Sommerfeld  radiation  condition  and 
equation   (2.9.2)    follows   from  equation   (2.5)   and  the  boundary  condition 
(2.4.2). 
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For  region  3: 


"    +  11- 
T      p  dp 


"Y  -  Y  ]  HA   (p.a.u) 

P 


4>, 


=     0 


(2.10) 


H     (0,a,a))     =     0 
<p3 


[p  H"   (p.a.a))] 


dp  <}) 


=     0 


p=a 


(2.10.1) 
(2.10.2) 


Equation   (2.10.1)   and  (2.10.2)   result  from  imposing  (2.4.1)   and  (2.4.2) 

For  region  2: 

In  this  region  the  magnetic  intensity  will  be  written  as  a 
linear  combination  of  the  incident  TEM  magnetic  intensity  and  another 

term    \p     (p.z.w)    . 

(f>2 


H<t,2(p'Z'w)   =     HJ,  +  \b'z^ 


(2.11) 


The  term     ij>     (p.z.w)    represents  obviously  all   of  the  reflected  fields  that 
exist  within  the  feedline,      on      z  <  0   .     For  the  other  extreme  of 
z  >  0   ,     ty.    (p,z,co)  must  account  for  the  fields  which  exist  in  free  space 
from  sources  within  the  feedline.     It  must  be  representable  as 

iK    (p,z,o))  =  {-H1.  +  other  field  terms)    ,     z  >  0 

since     h"!     cannot  be  present  in  the  half  space     z  >  0     because  it  vio- 
lates  required  boundary  conditions. 
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Since  the  intensity     H,    (p,z,co)     must  be  a  solution  to  the  dif- 

ferential  equation   (2.4)   the  terms     H       and     ty     (p,z,w)     must  satisfy 

the  following: 


and 


V 


3z 

f 

9z2 


[- 


3p 

aP2 


—  | — -  +  k2J  h!     -     0 

p   3p  1  c|> 


18  1  .2-,    ,      ,  v 

p  ^F"7     k  ]  %(p>z'w) 


=    o 


The  incident  field     H1^  =  !M-elkz     identically  satisfies   the 

4>        P 

first  of  the  above  equations.     Therefore  we  only  need  to  solve  the 
second  equation.     The  magnetic  field  intensity  for  region  2  will  be 
completely  known  after  this   is  done. 

The  above  differential  equation  for      ty.    (p,z,w)  will  now  be 
transformed  by  multiplying  through  by  e  a      and  operating  with 

r°°  /2T 

dz   . 


vEtF 


2  2 

9z2     h  dpZ     p  8p    P2  ^ 


Integrating  by  parts  twice  and  requiring  that     \\>.   (p,z,u)     and 

-j— ip     (p,z,w)  ■>  0     as     |z|  -»■  °°  ,  we  derive  the  differential  equation 
az     $2 

for  region  2. 


r    d         18         1  2-,    .      /  » 

dp  P  t 


dp 


(p     (Jj        (p  ,01,(1))) 


=     0 


p=a,b 


(2.12) 


(2.12.1) 
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The  condition  expressed  in   (2.12.1)   follows  when  it  is   recognized  that 


iK>  =  ° 


(2.12.2) 


and  use  is  made  of  the  boundary  condition  expressed  in   (2.4.2)and  (2.11). 

It  was  previously  stated  that  the  incident  magnetic  intensity 
does  not  satisfy  equation (2.8.5).     In  the  derivation  of  (2.8.5)   it  is 
essential   that  the  function  and  its  derivative  vanish  as     |z|   -*  °°     . 
It  is  obvious  for  slightly  lossy  media  that     h"!     does  not  vanish  as 
z  ■*  -°°  ,   therefore  it  cannot  satisfy  equation   (2.8.5). 

Upon  multiplication  of  equation   (2.11)  by    e  a       and  operat- 


ing with 


Jlv 


dz   ,  the  following  results: 


H.   (p,a,u>)  =  —     f  (H1.)  e1az  dz+<K  (p.a.aj) 


(2.12.3) 


Substituting  the  value  of    h"!     in  the  above  integral   produces: 


M(co) 


e1kZ  J<*  H7 
e         dz 


(2.12.4) 


This  integral   does  not  converge.     If  it  is   assumed  that  the 
medium  is  slightly  lossy,  we  can  write  the  quantity     k     as   a  complex 
number 


k  =  k,  +  ik2 


k,  >>  k?   ,       k,,k?  >  0 


The  loss  term     k?     will  be  allowed  to  go  to  zero  later  in  the  problem 
to  yield  results   for  a  lossless  medium.     The  value     k     is  specifically 
chosen  in  this   form  to  ensure   convergence  of  the  TEM  magnetic  intensity 
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as     z  -*■+<»  .     This  choice  of    k     ensures  convergence  of  (2. 11. 4) at  the 

upper  limit  (z  -*>  +°°)    .     Therefore  it  is  concluded  that  we  should  operate 

r  r 

on  (2.11)with  the  integral        dz     rather  than         dz   .     When  this  is 
done,  the  following  equation,  which  will   be  usTd  later,   results  : 

Ht  (p.a.u)   =  +  ^ +  /  (p.o,u)  (2.12.5) 

h  /2Tp(a  +  k)         *2 


In  the  differential  equations  for  the  various   regions  of  the 
problem  we  have  made  use  of  a  quantity 


2        2.2 

Y     =  a     -  k 


The  square  root  is 


Y 


=  Jo-  -  kZ  (2.13) 


which  is   a  multivalued  function  of  the  complex  variable     a     defined  on 

a  two  sheeted  Riemann  surface.     The  value  of    y    must  be  uniquely 

defined  so  that  only  the  proper  branch  of  this  double  valued  function 

is  used.     The  branch  cuts  for  this   function  are  determined  in  Appendix 

B  and  depicted  in  Figure  2-3. 

The  necessity  of  satisfying  the  Sommerfeld  radiation  condition 

requires  that  the  asymptotic  frequency  domain  behavior  of  the  magnetic 

intensity  for  the  space  around  the  coaxial  waveguide  antenna  (exclusive 

of  regions  2  and  3  with     z  <  0)     be  that  of  an  outward  traveling  or 
evanescent  spherical  wave. 

4ik|r| 


H,  n,  f 


* 


r| 
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Figure  2-3.     Branch  cut  for  multivalued  function     y 
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|r|  -JF77 

If  the  value  of    z     is  such  that     z  »  p     the  asymptotic  behavior 
specializes  to 

Jklzl 

H      (p,z,o>)     ^     —r—, —  ;  Z  ^  ±00 

<P|  |Z| 

eikz 
\\)     (p.z.w)     ^    ;  z  +  +0° 

f>2  z 

eikz 
H^  (p.z.u)     ^    —j-  ;  z  ■+  ■*» 

In   region  2  with     z  <  0   ,  the  magnetic  intensity  is  known  [17, 
p. 327]  to  consist  of  the  TEM  mode  and  higher  order  E       modes.     Only 
E      modes  are  present  due  to  the  symmetry  requiring  no  variation     with 
the  angle     <j>  .     The  only  field  intensity  components  allowed  are    H     , 
EQ     and     Ez  . 

From  equation   (2.11)  we  have 

ce-ikz 


%(p'Z'w)     =  p 


pi ■ — 2~ 

I  Con   iJlKnP)Yo^na)-Yl^nP)Jo^na^e"1Z  "  "^ 


where  x/k2-   (y^)2  =     i/(Yon)2-  k2     ,     n=l,2,3,. 


y'       is  the  nth  ordered  root  of  the  auxiliary  equation 
on 


J  (v1  b)  Y  (y1   a)    -    Y  (y*  b)  J  (y1   a) 
owon  '     ovlon   '  owon   '     owon   ' 
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C     and    C        are  constants   relative  to  the  spatial   variables 
on  r 

p     and    z  .     Jv(z)     and    Yv(z)     are  the  ordinary  Bessel   function  of  the 
first  and  second  kind,   respectively  of  the  vth  order. 

In  region  3  with     z  <  0   ,  the  magnetic  intensity  is  also  known 
[17,  p. 322]  to  be  the     E      modes  of  a  cylindrical  waveguide. 


T3  n=  1 


pu„„  .1zyk2.  fa* 


on\  „      J  "a 


where  7k2  -  (^f-)2  =     J C~f-)2-  k2         ;  n=l,2,3,..-     ,     CQn     is  a 

constant  relative  to    p     and     z     and     u        is  the  nth  ordered  root  of 
the  auxiliary  equation, 

J0(u)     =    0 

Given    k  =  k  +  ik     ,      and    the  asymptotic  behavior  of  the  field  in 
region  1,  it  is  apparent  from  the  earlier  discussions  that     H     (p,a,w) 
is  analytic  in  the  strip 

-k2  <  x  <  k2 

of  the  complex    a    plane. 

For  regions  2  and  3  we  must  consider  two  cases,  the  first  case  is 

when  the     E        modes  within  the  region  are  evanescent  and  the  second 
on 

case  is  when  they  are  not. 

Case  1:  k  <  y'm  region  2 

u 
k  <  -^  region  3 
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Case  2:  k  >  y'  region  2 

u 

k  >  -20.  region  3 


For  Case  1   the     E„„  modes  in   regions  2  and  3  are  evanescent, 
on  3  ' 

therefore  the  asymptotic  value  of  the  magnetic  intensity  as     z  ■*  -°° 
is  the  TEM  mode  portion  of    ty,       within  region  2  and  zero  within  region 
3.     We  conclude  that  in  the  event  of  Case  1, 


iK    (p,a,a))       and     H     (p,a,w) 

are  analytic  with  the  strip     -k2  <  t  <  k2     of  the     a    plane. 

For  Case  2,   the     E      modes  are  not  cut  off  and  the  asymptotic 
forms  behave  as 


i—     ~2  /        ~Z 

-izk/l  -  £j  -izCk^  ik2)yi  -  \ 

e  =  e  ,     z  -»■  -°° 

This  behavior  produces  a  strip  of  analyticity  that  is  at  least  given  by 


-k2  <  t  <  k2  Re(7l   -   (|)2) 

The  symbol     Re     means  the  real  part  of  the  quantity  that  it  precedes. 
Although  at  this  point  there  is  no  apparent  difficulty  with  this   result, 
the  inclusion  of  values  of     ka     and     kb     such  that 

on 
and 

kb     =     u 


om 


where     u         is   the  mth   ordered  root  of     J   (u)  —  J   (r-  u)   =  0     requires 
om  o  o  b  ^ 

special   attention.     It  will   be  shown  later,   that  the  field  intensity 
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outside  the  waveguiding  structure  is  zero  at  these  frequencies. 

In  light  of  the  above,  we  restrict  ourselves   to  the  fundamental 
mode  by  imposing  the  conditions 

kb  «  1  and  ka  «  1 

when  we  are  considering  a  time  variation  other  than  harmonic  time 
dependence.     Under  this   restriction  Case  2  does  not  occur,    and 
the  functions 

H^   (p.a.u)      ,       ty     (p,a,co)      and     H     (p,a,u) 

are  analytic  in  the  strip 

-k2  <  t  <  k„ 

of  the  complex  a  plane,  see  Fig.  2-4. 

When  we  consider  harmonic  time  dependence,  the  strip  of  analy- 
tic!* ty  is   changed  to 


k2  <  t  <  k2  Rejl  -  (l)2 


"2 

where     6  =  y'Qn     or    -~     for  region  2  or  3,    respectively,  n=l  ,2,3,  •  •  • 

The  value  of    6     is  the  magnitude  of  the  cutoff  wave  vector  of  the 

highest  order     EM     mode  which  is  non-evanescent  in   reqion  2  or  3,   as 
on  oi 

appropriate.     Also     ka  ¥  u         and     kb  4  u 

on  om 

To  obtain   a  unique  solution   to  the  problem,   it  is  necessary  to 
determine  the   algebraic  behavior  (as     a  -»■  °°)   of  the   factors   in  the 
completely   factored  and  decomposed  Wiener-Hopf  equation   [4,  p. 37]. 
This   determination  is  made  by  considering  the  edge   conditions   [19]. 
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I'T 


a  plane 


Figure  2-4.     Strip  of  analyticity  for    H     (p,a,w), 
\l>.    (p,a,aj)      and     H     (p.a.w) . 


*c 


♦, 
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For  the  infinitesimally  thin  conductors  utilized  in  this  problem,   the 


edge  conditions   are 

-1/2  + 

Ez(p,z,u))  <*  z     '  as         z  ->-  0     at     p  =  a  and  b 

Equation   (2.5)  states  that 


1    r) 
Ez(p,z,o>)     <*    -^[pH,(p,Z,iii)] 

The  point  of  interest  is  that 

-^  [p  H^p.z.w)]-  z"1/2  (2.14) 

as     z  -►  0       at    p=aorb.     z  ->-  0       means  that     z  =  0     is  approached 
from  the  direction  of  positive  values  of    z  . 


Development  of  the  Wiener-Hopf  Equations 

In  this  section  the  Wiener-Hopf  equations  will  be  developed.   The 
solution  to  these  equations  will    lead  to  the  spatial   transform  of  the 
desired  field  intensities. 

We  begin  by  writing  solutions  to  equations   (2.9),   (2.12)  and 
(2.10)  in  that  order. 

H.    (p,a,u>)     =     A(a,oj)   K,(yp)  (2.15) 

(p.,  i 

iji     (p,a,u>)     =     B(a,u>)    I^yp)  +  C(a,u>)   K^yp)  (2.16) 

H.    (p,a,w)      =      D(a,u))    I,(yp)  (2.17) 

The  form  of  the  solution  to  (2.9)  expressed  in  (2.15)  was  chosen  so 
that  the  asymptotic  behavior  would  satisfy  (2.9.1) ,   thereby  satisfying 
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the  Sommerfeld  radiation  condition. 

The  solution  to  (2.10)  expressed  in   (2.17)  does  not  include  the 

first  order  modified  Bessel   function  of  the  second  kind  since  it 

becomes  infinite  as     p  ■>  0     a  condition  that  violates  equation  (2.10.1). 

The  function     I,(yp)  -»■  0     as     p  -*•  0   ,  so     H     (O.a.oo)  =  0     and  (2.10.1) 
i  q>3 

is  satisfied. 

The  coefficients     A(a,oo),  B(a,co),   C(a,w)     and     D(a,aj)     are  un- 
known at  this  point.     It  is   the  task  of  finding  these  coefficients 
which  leads  to  the  Wiener-Hopf  equations.     Each  of  the  coefficients 
will  be  expressed  in  terms  of  two  heretofore  unmentioned  functions, 
W  (a,u))     and    H   (a,w)    ,  which  will  be  found  by  employing  the  Wiener- 
Hopf  procedure.     To  this  end  we  will  make  use  of  the  relationships 
which  exist  among  the  magnetic  intensities  and  their  derivatives   (with 
respect  to    p  )     evaluated  on  the  boundaries  separating  the  three  regions 

At    p  =  a    we  have 

H ,   (a,z,u>)  =  H     (a,z,a))         for        z  >0 

Therefore, 

H*   (a,a,oo)  =  H*  (a,a,oj)  (2.18) 

From  the  analogous   condition  on  the  magnetic  intensity  at     p  =  b     for 
z  j>  0   ,  we  have 

H^   (b,a,w)   =  H^   (b,a,co)  (2.19) 

<P-|  <P2 

These  equations   relate  the  magnetic  intensities  at  the  bound- 
aries and  will  be  used  at  a  later  stage  in  the  development  which  leads 
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to  the  Wiener-Hopf  equations. 

It  was  assumed  earlier  that  the  coaxial  waveguide  antenna  struc- 
ture was  made  of  a  perfectly  conducting  material.     Such  material   does 
not  support  a  tangential  electric  field,  therefore 


and 


E   (a,z,co)     =     0 


E,(b,z,w)     =0  if 


z  <  0 


E  (p,z,to)     is   a  continuous   function  of    p  ,  thus     E  (p,z,co)     has  the 
same  value  at  a  given  value  of    p     regardless  of  the  direction  of 
approach,  whether  from  smaller  or  larger  values  of    p   . 
Equation  (2.5)  states  that 

E,     -    IjL  (ph  ) 

Z  p    dp  <(>' 

When  this  is   combined  with  the  preceding  discussion  about  the  tangen- 
tial electric  field,   the  results  are 

£-[pHA   (p.z.o.)]     ■     4-[pH.    (p.z.o))] 


dp  LM  <J) 


dp       ^      tyr 


at    p  =  a  .     By  a  spatial   Fourier  transform  (see  equation   (2.8))   of 
these  equations,   it  is   found  that 


[pH,    (p,a,u))] 


^Lpn4> 


[pH,    (p,a,co)] 


p=a 


3p"  Lpn* 


p=a 


When  we  consider  that     E  (a.z.w)  =  0     for     z  £0   ,  this  equation  re- 
duces to 
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Hp"  tpH^Cp^^)] 


d      r     ,+ 


P=a 


dp  CP^   (p.ot,u)] 


(2.20) 


p=a 


Note  that  we  have  used  equations   (2.11)   and  (2.12.2)   to  obtain  (2.20). 
By  very  similar  arguments  with     p  =  b    we  also  find  that 


dp-  Cp"J  (P"0»u)] 


P=b 


^p-[p^2(p,a,u))] 


p=b 


(2.21) 


We  now  define    W   (a,u)     and    H  (a,w)     as   follows 


W    (a,co)      =     -g-     [pH      (p,a,co)] 
H   (a, go)     =     -g-  [pH     (p,a,to)] 


p=a 


p=b 


These  functions  of    a    are  analytic  in  the  upper  half  plane  given  by 
x  >  -k2   . 

Returning  to  the  spatial   frequency  domain  field  solutions  given 
by  (2.15),   (2.16)   and  (2.17),  we  multiply  each  of  the  equations  by     p 
and  take  the     p     derivative  of  the  resulting  equations  [20,  p. 376]. 


4r  [pH+  ]  +  4-  [pHT  ]  =   -YP  A(o,u)   K  (YP)  (2.22.0) 


dp  LM  ^J       dp  LK  <))-| 


h  Cp*J  ]  +  3p~  [p^  ]  =  wCB(o.m)    Iq(yp)   -  C(o.»)   Kq(yp)]     (2.22.1) 


i  K,]  +  ^  [pv]  =  yp  D(a'w)  !o(yp) 


(2.22.2) 
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By  taking  the  last  of  the  above  equations  and  evaluating  it  at     p  =  a  , 
it  is  seen  that 


dp  LK  4>3J 


p=a 


=  W   (a, to)     =     ya  D(a,u))    I    (ya) 


or 


D(afW)     =     ^-^ 


a7 


(2.23) 


Note  that 


dp  LH  (J),,-1 


=     0 


p=a 


since, 


E   (a,z,tjo)  =  0         for        z  <_  0 


Recall   that 


Bp"  [pVj(p»z»")J  =  ^  CpH^tp.z.to)]  =  %-  [p*^(p.2,»)] 


at     p  =  a    . 


(This  may  be  easily  seen  to  stem  from  equation 
(2.12.2)  and  the  continuity  of  the  tangential  electric  field     E  ). 
Since  the  tangential   electric  field  is  zero  at     p=a  or  b     for     z  <.  0, 


and 


3p-[p^2(p,z,u))] 


^  [p*'  (p,a,o»)] 


0         for        z  <  0 


P=a,b 


p=a,b 
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Combining  this  result  with  equation  (2.20)  and  equation   (2.22.1) 
evaluated  at    p  =  a     leads  to  the  following: 

W+(a,u>)  =  Ya[B(a,a>)   IQ(Ya)   -  C(a,co)   KQ(Ya)] 

or 

W  (a,co)   -  YaB(a,(Jo)   In(Ya) 

c(^>    ■" WkJTZ —  (2-24) 

Substitution  of  the  values  of     D(a,w)     and     C(a,ca)     given  by  equations 
(2.23)   and  (2.24)  into  equations   (2.17)   and  (2.16),   respectively,  yields 


W  (a, to)  -  B(a,to)Ya  I   (Ya) 


and 

\p     (p,a,u)   =  B(a,w)    I1(Yp)-[ ya  K   (   a) ]   ^(Yp) 

Evaluation  of  the  irmiediately  preceding  equations  at     p  =  a     and  use  of 
the  Wronskian  relationship  produces  the  following  intermediate  expres- 
sions: 

+                                             W+(a,oa)   I1(Ya) 
H^(a,a,o>)  +  H^(a,o,u)  ■    ya  1qM 

+  B(a  u)        M  (a»u)  K-jCya) 

♦^(a.o.u))  +  ^(a.a.u)     =  Ya  KJ[yi)  "      Ya  kQ(Ya) 

By  subtracting  the  lower  of  the  two  equations  immediately  above  from 
the  top  equation,  we  find  that 
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H^  (a, a, to)  -  \p     (a, a, to)  +  H~  (a, a, to)  -  ijT  (a, a, to)  = 

O  €L  O  C- 

W    (a,(o) 6(01,0)) 

(Ya)2I   (Ya)  K  (Ya)  "  ^  KQWa) 


The  Wronskian  was  used  to  achieve  some  simplification  in  the  above 
equation  . 

If  we  evaluate  (2.12.5)  at    p  =  a    we  find 


H     (a, a, to)  -  ty.    (a, a, to)     = 


JM(to) 


h  h  m  a(a+k) 

But  it  is  known  from  (2.18)  that 


HA  (a, a, to)  =     HA  (a, a, to) 
Therefore, 


I,  U|U|IU  II, 

4»3  <D2 


H,   (a, a, to)  -  \\).   (a,a,to)  = 


iM(to) 


h  m  a(a+k) 


and  it  follows  that 


iM(to)  ,,-    ,  \        ,-   /  \  W  (a^to) B(g,to) 

lili^ —  +  H  .   (a,a,to)   -  \\i.    (a, a, to)  =  *       '   ' 

ffi  a(a+k)         ^3  h  (Ya)^I0(Ya)K0(Ya)      Ya  KQ(Ya) 

(2.25) 
Evaluating  (2.22.0)   at    p  =  b    we  see  that 
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Note  that  the  boundary  condition 
invoked.     This  allows  substitution  of 


Ez(b,z,k))  =  0  for  z  <  0  was 


W^ 


p=b 


into  equation  (2.22.0). 

When  (2.22.1)  is  evaluated  at  p  =  b  it  becomes 


or 


H+(a,a))  =  Yb[B(a,a))  I^yb)  -  C(a,u)  KQ(Yb)] 

H+(a,w)  -  B(a,u>)  yb   I  (yb) 

C(a,co)  = — .,  /  .  > 

Yb  KQ(Yb) 


(2.27) 


If  we  replace  the  coefficients  in   (2.15)   and  (2.16)  by  the  ex- 
pressions given  in   (2.26)   and  (2.27)   respectively,   then  it  is  apparent 
that 


H,    (p.a.oj)     = 


H   (a, to)  ■/   /      \ 

YbTJTb)    Ki(Yp) 


and 


ty     (p,a,u) 

q>2 


=  B(a.«)   I^yp)  "  [H+(a>UiYD'BK(^jYb  ^^^(yp) 


By  evaluation  of  the  preceding  equations  at     p  =  b     and  simplifying, 

we  get 

+                                                       H+(a,oj)   K^Yb) 
H^(b,a,u>)  +  H^(b,a,a))  =  -       yb  Ko(Yb) 

_     B(a.M)  "+{^   K1W 


^  (b.a.w)  +  ^(b.a.w)   =  ■  yb  kJYb) 


Yb  ko(Yb) 
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Subtraction  of  the  latter  from  the  former  yields 

H^(b.a.«o)  -  ♦^(b,a,Wl+H^(b.ai»)-^(b.aiU)  =  -  ^f$y 
It  is  seen  that 

H+  (b.a.u)  -  i,    (b.a.u)  =  -^ 

♦l         *2  ^iF  b(a+k) 

by  evaluating  (2.12.5)   at    p  =  b     and  using  the  equality  given  by 
(2.19).     Substitution  of  this   result  leads  to 

iM(a))         +  HT   (b,a,U)   -  #:  (b.a.co)   ■  -  StJ^E]  (2-28) 

m  b(a+k)         *1  h  Yb  Vyb) 

Equation   (2.28)  is   analogous  to  (2.25).     These  two  equations  are  the 
result  of  matching  the  tangential  electric  intensity  (-°°  <  z  <  °°)     and 
the  magnetic  intensity  (z  >^  0)  at  the  common  boundaries  between  the 
regions  into  which  the  space  was  earlier  divided.     These  equations   are 
also  combined  to  form  one  of  the  Wiener-Hopf  equations  by  eliminating 
B(a,u))  between  them.     The   resulting  Wiener-Hopf  equation  involving 
W  (a, id)  is 

x       [H~  (b,a,u>)   -  \\)~   (b.a.w)] -ya  K  (ya)[H~  (a.a.w)  -  \K  (a,a,w)] 

■(f).  (Pp  u  ro  ro 

(2.29) 

Since     C(a,w)  in  this  problem  is  unique,  the  two  values--one 
given  in  expression   (2.24)   and  the  other  given  in  (2.27)  must  be  equi- 
valent.    When  these  expressions   are  equated  the  results   relate  the 
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two  variables     W  (a, to)     and     H+(a,oo)    , 

+  a  UYa)      +  I.(yb)       I   (ya) 

W  (a'w)  =  b  KQ(Yb)  [H  («•«)] -B(«.«)[yYBT-  y^iy^a  KQ(Ya) 

The  terms     W  (a, to)     and     B(a,to)     appearing  in  the  above  expression  are 
eliminated  by  using  expressions   (2.28)   and  (2.29),  giving  a  Wiener- 
Hopf  equation  involving  the  function     H+(a,to)    , 


1tM<"?    [I0(yi)  -  I  (yb)]  =    H+<a^ 

v?F(a+k)       °  °  Yb  K.(Yb) 


+  yb   I0(Yb)[H^  (b,a,w)   -  f   (b.a.to)] 


-  Ya  IQ(ya)[H~  (a, a, to)   -  f.   (a, a, to)]  (2.30) 


In  the  coaxial  waveguide  antenna  we  usually  have     b-a  « 


a+b 


2 
to  ensure  TEM    mode  excitation  at  microwave  frequencies.     For  this  case 

the  field  intensities  will  be  rapidly  varying  when  compared  to      p    ,and 
the  fields  in  the  annular  region  will   approach  those  of  a  parallel  plate 
waveguide  [17,  p. 329]  of  width   (b-a). 

For  the  parallel  plate  waveguide  it  is  well   known  from  the  pre- 
servation of  symmetry  that  for  a  TEM  mode  excitation  the  current  dis- 
tribution on  the  two  parallel  plates  are  equal   in  magnitude  and  oppo- 
site in  direction  [21,  p. 126].     For  excitation  by     E      mode  fields, 
the  currents  on  the  two  parallel  plates  will  be  equal   in  magnitude  and 
in  the  same  direction  for  odd  n,   and  in  opposite  directions   for  even  n. 
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Up  to  this  point  no  approximations  have  been  made  and  equa- 
tions  (2.29)   and  (2.30)  are  exact.     It  is  recognized,  though,  that  the 
form  of  these  equations  is  such  that  they  cannot  be  decomposed  into 
two  expressions,  one  analytic  in  the  upper  half  plane  and  the  other 
analytic  in  the  lower  half  plane  with  an  overlapping  strip  in  which 
both  functions  are  analytic.     This  is  made  apparent  by  observing  that 
the  coefficients  multiplying  the   "minus"  functions  on  the  right  hand 
side  of  each  equation  are  not  identical.     These  same  coefficients  also 
contain  a  branch  point  and  are  therefore  not  of  the  form  where  an 
exact  solution  is  possible  [4,  p. 153]. 

We  note,  however,   that  if    b-a  «  ^*—    then  the  field  inten- 
sities approach  those  of  a  parallel  plate  waveguide  and  if  TEM  mode 
excitation  is  impressed,  the  currents  on  the  cylindrical  waveguide 
conductors  are  approximately  equal  in  magnitude  and  opposite  in  direc- 
tion.    The  mathematical  expression  of  this  condition  applicable  to  the 
structure  under  consideration  is 


2irb[H.    (b,z,u>)  -  H.    (b.z.w)]  *  2iTa[H,    (a.z.w)   -  H     (a.z.w)] 

9-|  h>2  ?3  *2 

for    z  <  0   .     We  now  recall  equations   (2.11)   and  (2.7): 

H.    (p,z,w)     =     H1.   +  ifi     (p.z.u)  (2.11) 

i     m     M^ikz  (2jr) 

<J>  P 

Substitutions  of  these  quantities  in  the  above  equations  yields 


-39- 


bEH^  (b.i.fti)   -  ^  (b.z,«)   -^-eikz] 

-  a[H     (a, z,03)   -  f     (a,z,u>)   -  ^Meikz]  for     z  <  0 


Under  the  assumption  that  the  approximation  is  very  good,  we 
can  equate  the  two  currents  and  obtain 

b[H     (b,z,u)  -  ty     (b,z,o))]  =  a[H     (a.z.co)  -  ip     (a,z,u>)] 

9-j  q>2  4>3  $2 

for     z  £  0 

When  this  equation  is  spatially  Fourier  transformed  the  results  are 
precisely  what  is  needed  to  allow  decomposition  of  the  two  Wiener-Hopf 
equations   (2.29)   and  (2.30),  i.e., 


b[H.    (b,a,to)   -  i|i~   (b,a,w)]  =   a[H     (a,a,u))   -  \\>,    (a, a, to)]  (2.31) 

9-.  (f>2  <Pt  $2 


when 


b-a«-^-  (2.32) 


It  is  also  noted  that  the  expression   (2.31)   is  analytic  in  the 


lower  half  plane  given  by     t  <  k2     or    I^Ret/l  -  ($/k)   )     depending 
on  the  excitation.     It  is  convenient  to  define     S"(a)     such  that 


S"(a)  =  b[H"   (b,a,u))   -  iff   (b.a.u)]  =  a[H"   (a,a,oj)-^"   (a, a, to)] 
9-j  ^2  ^3  ^2 


(2.33) 


To  derive  decomposable  Wiener-Hopf  equations  we  substitute 
(2.33)   into  (2.29)   and  (2.30).     The  resulting  equations  are: 
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(2.34) 


iM(w)                             W+(a,u)) 

i    c°*  / ,- 1  ^        n 

v^?(a+k)       Y2aI0(Ya)[Ko(yb) 

-  KQ(Ya)] 

iM(o>)                              H+(a,w) 

0 

+  S"(a)  =  0 

(2.35) 
v^?(a+k)       y  bK0(Yb)[I0(Yb)  •    IQ(Ya)] 


It  is  interesting  to  note  that  if  the  radius  of  curvature  is 
large  we  may  imagine  the  conductors  in  Figure  2-1   to  be  cut  along  the 
z  direction  and  flattened  into  parallel  planes.     Therefore  it  is 
reasonable  to  expect  that  equations   (2.34)  and  (2.35)  will  become  the 
Wiener-Hopf  equations   for  two  parallel  planes  in  the  limit  as   radii     a 
and    b     approach  infinity.    The  planes  would  be  b-a  distance  apart  and  the 
axis  of  symmetry  would  be  displaced  a  distance     a^b  -*■  °°     from  the 
coordinate  z  axis.     By  taking  the  limit  as     a  and  b  ->  °°     and  using  the 
asymptotic  forms  for  the  modified  Bessel   function,  equations   (2.34) 
and  (2.35),   respectively,  become 

ifia) w!m  [ n_Ya -,  .  s-(a)  =  0 


v£?(a+k)  Ya  -y±r^-  Y(b_a) 

e         fc       sinh  -^ — '- 


and 


iMM +  H+M    [ ,,       .   Yb ]   -   S"(a)    =  0 

ma+  k)  y2b  _Yibzai 

e  sinh  J-i2 — *- 

It  can  also  be  shown  from  the  definition  of     "  (a,w)     and     H   (a,w) 
that  for  a  parallel  plate  case     H+(a,u))  =  -  U  (a,w)   .     Therefore  the 
equations   (2.34)   and  (2.35)  become  identical   in  the  limit  as  the  radii 
approach  infinity  and  each  is  the  Wiener-Hopf  equation  for  the 
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radiation  from  an  open  ended  parallel  plate  waveguide  with  TEM  mode 

excitation  [4,  p. 107  or  21,  p. 128]. 

The  solution  of  equations   (2.34)   and  (2.35)   for     ^(a.to)     and 

H  (a,w)     allows  a  complete  characterization  of  the  magnetic  intensity 

in  the  spatial   frequency  domain.     The  foregoing  work  already  contains 

the  equations   relating     A(a,w)     and     D(a,u>)  to    W   (a,w)     and    H  (a,oo), 

These  coefficients  are  given  by  (2.26)  and  (2.23),   respectively.   To 

obtain  the  remaining  coefficients     C(a,oo)     and     D(a,w)   in  terms  of 

W(a,oo)     and    H  (a,u))    ,  we  return  to  equation   (2.22.1)   and  evaluate 

it  at    p  =  a  and  b   .     When  these  results  are  considered  with  (2.20), 

(2.21)   and  the  relation 


■^  [p*^(p.a,to)] 


=  0 


p=a,b 


we  get 


W  (a.oj)  =  ya[B(a,u)Io(Ya)   -  C(a,w)Ko(ya)] 


H+(a,w)  =  Yb[B(a,a3)I0(Tb)   -   C(a,w)K0(Yb)] 


The  solutions  to  these  simultaneous  equations  yield     the  coefficients 

we  seek, 

H+(a,o))Kn(Ya)       W+(a,oj)Ko(Yb)  1 

] [I   (Yb)K  (Ya)-I   (Ya)K  (Yb)] 


B(a,co)  =  [- 


Yb 


Ya 


C(o,w)  =  [- 


H+(a,w)I   (Ya)     W+(a,co)I   (Yb) 


o' 
(2.36) 


Yb 


"va  ] CIJYb)K  (Ya)-I   (Ya)K  (Yb) ] 


o' 
(2.37) 
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To  complete  the  list  we  add: 


«->   ■  ^T  (2-23) 


We  return  to  the  earlier  discussion  that  certain  values  of 
ka.kb  are  excluded.     Namely 


and 


ka  =  u™ 
on 


kb  =  u 

om 


where  u_  is  the  nth  ordered  root  of    J   (u)  =  0 

on  o 

and  u_  is  the  mth  ordered  root  of    J  (u)  -  J  (£-  u)  =0 

om  o  o  b 

There  is   an  important  mathematical    reason  for  excluding  these  particu- 
lar values  of  ka  and  kb  that  shows  up  in  the  factorization  procedure 
of  Appendix  D.     The  physical    reason  for  their  exclusion  is  that  at 
these  frequencies  the  magnetic  intensity  everywhere  outside  the  struc- 
ture is  zero. 

If     ka  =  u       ,  then  the  magnetic  intensity  in  the  region 
p  <  a   ,  z  <  0     must  be  zero.     This  value  of    ka     corresponds  to  the 

cutoff  wavelenqth  for  the     E        mode  of  a  circular  waveguide  of  radius 
3  on 

a  .     Exactly  at  cutoff,  the  magnetic  intensity  for  the     E        mode  in 
the  ci  rcular  waveguide  is  zero  [17,  p. 322]. 

Now  considering  that     b-a  «  ^*—    the  structure  has   fields   that 
behave  approximately  as  those  of  a  parallel  plate  waveguide.     In  this 
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limit  we  have  shown  above  that  the  field  exterior  to  one  of  the  plates 

is  zero.     By  symmetry  the  field  that  is  on  the  exterior  side  of  the 

other  plate  must  also  vanish.     By  continuity  it  is  concluded  that  the 

field  is  zero  everywhere  exterior  to  the  waveguide  channel. 

For  the  other  case,     kb  =  u       ,  we  examine  the  asymptotic 

behavior  of     I  u     I    as     m  ■*  °° 
1    om' 


J0<u>  "  Jo(F  ">  *Jh cos(u  "  T>  \H~  cos(^  u  "  t) 


TTg-U 


«/§j(cos(u  -})  -  cos(f  u  -\)) 


with     a  and  b  of  the  same  order  of  magnitude  as     u  ■*■  °°       .     Through 
the  trigonometric  identity 

cos  x  -  cos  y  =  -   2  sin(^-)  sin  -^- 
it  is  determined  that  the  zeros  are  approximately  located  at 

for     u      ■*■  °°     or    m  -*■  °°  . 
om 

The  asymptotic  values  of     u        given  above  are  recognized  as 
the  asymptotic  values  of  the  cutoff  frequencies  as     m  ->  °°    for  the 

E        modes  of  a  circular  waveguide  of  radius     -5-    and  the  even 

om  c 

ordered     E        modes  of  a  parallel  plate  waveguide,  with  a  channel 
om 

width  of    b-a   ,   respectively.     The  second  value  of     uQm     also     cor- 
responds to  the  asymptotic  value  of  the  even  ordered     EQm    modes  of  a 
coaxial  waveguide  [22].     In  the  event     b-a  «  ^-,  we  may  approximate 
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the  structure  as  a  circular  waveguide  of  mean  radius     ^-    or  a 
parallel  plate  waveguide,  depending  on  the  excitation.     At  the  cutoff 
frequencies  for  such  structures,  the  electromagnetic  fields  are  iden- 
tically zero.     Therefore  we  have  excluded  the  use  of  such  frequencies, 
since  they  invoke  no  response. 

Solution  of  the  Wiener-Hopf  Equations 

The  two  Wiener-Hopf  equations   (2.34)   and  (2.35)  are  solvable  in 
the  conventional  Wiener-Hopf  fashion.     The  various   terms  are  factored 
into  the  product  of  a  "plus"  and  "minus"  function,  where  required,  and 
the  resulting  equation  is  decomposed  into  the  sum  of  "plus"  and  "minus" 
functions.     [The   "plus"  ("minus")   functions  are  regular  in  the  entire 
upper  (lower)  half  plane.]     Analytic  continuation  arguments  are  then 
invoked  to  yield  the  solution. 

To  preserve  the  continuity  of  the  work,  Appendix  D  was  devoted 
to  the  factorization  [23]  of  the  kernels 


{l0(ya)[K0(Yb)  -  Ko(Ya)]}_1  =  x(a)  =  x+(a)   x"(a) 


and 


{K0(yb)[Io(Yb)  -  I0(Ya)]}_1  -  y(a)  =  y+(a)  y"(a) 


In  this  section  we  will   use  the  expressions     x  (a),  x"(a),  y  (a),  y"(a), 
Their  specific  functional   form  is  found  in  Appendix  D. 

We  rewrite   (2.34)   and  (2.35)  using  the  newly  defined  symbols 
for  the  kernels  that  must  be  factored, 
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iM(w)_    +     W^}_  [x+(o)   x-(o)]  +  s-(o)   __  Q 


v£7(a+k)  y  a 


38) 


1M(u)         +i^i-    [y  +  (a)y-(a)]  +  S-(a)   =     0  (2.39) 

»£¥  (a+k)  yb 

Recall   that  the  value  of    y       is  given  as 

2        2.2 

Y     =  a     -  k 

2 

When  this  value  of    y       is  substituted  in  the  above  equations,  we  can 


reduce  them  to 


and 


iM(o))(a-k)         +  W  (a,co)  x  (a)  +  S"(a)(a-k)  _  Q 
^2?(a+k)   x~(a)  a(a+k)  x"(a) 


iM(h))(a-k)  +  H   (g,M)  y   (a)    ,.  S~(a)(a-k)   _  Q 

v^?(a+k)  y"(a)  b(a+ k)  y"(a) 


Examination  of  the  region  of  analyticity  of  the  terms  in  the 
above  equations  reveals  that  the  middle  term  in  each  equation  is 
analytic  in  the  half  plane     t  >  -k2     and  the  last  term  in  each  equa- 
tion is  analytic  in  the  half  plane     t  <  k~   .     The  first  term  in  each 
equation  has  a  pole  at     a  =  -k  .     If  this  pole  were  not  present  it 
would  be  analytic  in  the  half  plane     t  <  k„   .       By  removal   of  the 
residue  of  this  pole,  we  obtain  the  completely  decomposed  equations: 
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W+(a>co)  x+(a)  _  i    IZ      kM(a)) 
a(a+k)  ™  (a+k)  x+(k) 


=  .  ^S~(a)(a-k)  +       iM(a))     j-a^k     +     2k     y{2AQ) 
x"(a)  ^2?(a+k)     x"(a)       x+(k) 


H+(a>a))  y+(a)   _    .    !Z_      kM(u) 
b(a+k)  ^   (a+k)  y+(k) 


=  .    |S'(a)(a-k)    ,        iM(hi)      j-a  -  k     +     2k     ]](2.41) 
y'(a)  v^2?(a+k)     y"(a)       y+(k) 


In  the  above  expression  we  have  used  the  result  from  Appendix  D,  that 

x  (a)  =  x~(-a) 

and 

y  (a)  =  y"(-a) 

The  terms  in  (2.40)  and  (2.41)  are  arranged  so  that  the  left  hand  side 
is  analytic  in  the  half  plane  x2  >  -kp  and  the  right  hand  side  is 
analytic  in  the  half  plane  Tp  <  kp  or  kp  Re(Jl  -  (§/k)  )  depend- 
ing on  the  type  and  frequency  of  the  excitation.  The  important  point 
in  this  discussion  is  that  there  is  a  strip  of  overlap  in  the  regions 
of  analyticity,  that  is,  both  sides  of  (2.40)  and  (2.41)  are  analytic 
in  a  common  strip 


-k,  <  x  <  k9     or        k?  Re  J~\   -  (3/k)2 
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We  therefore  conclude  that  each  side  of  the  expressions   (2.40) 
and  (2.41)   represents  the  analytic  continuation  of  the  other  side. 
Since  these  expressions  are  regular  throughout  the  complex  plane  they 
are  an  entire  or  integral    function. 

If  we  can  show  that  each  side  of  (2.40)  and  (2.41)   is  bounded 
as     j ot |   ■*■  °°     in  the  half  plane  for  which  that  side  of  the  expression 
is   regular,  then  we  can  use  Li  ouvi  lie's   theorem  to  set  each  expression 
equal   to  a  constant.     Li  ouvi  lie's   theorem  states  that  the  only  bounded 
entire  functions  are  constants. 

To  determine  the  limiting  behavior  as     |a|   -*  °°     for  the  terms 
in  (2.40)   and  (2.41),  we  examine  the  edge  conditions  given  in   (2.14) 
as 

as     z  ■*  0+     at     p  =  a  or  b   .     Taking  a  slightly  more  general   case,   the 
Fourier  transforms 


1 


and 


1 


zn  e1az  dz 


zn  elaz  dz 


determine  the  limiting  behavior.     In  these  integrals     n  <  0     to  have 


absolute  and  uniform  convergence  of  the  improper  integrals. 


It  is  also  well  known  that  if    |  f(z)|  <e 


■kpZ 


as     z   ► 


then 


lim  f(z)     =     lim     iaF  (a) 

z+o+ 


(lim  ia)(   lim  F  (a)) 


a-*-00 


a- 
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k2z 


for    x  >  -l<2     where     a  =  a+ix   .     Similarly,  if     |f(z)|    <     e  as 

z  -*•  -°°     then 

lim     f(z)     =     (lim  ict)     (lim  F"(a)) 


z+0" 


a  -»•-«> 


for    t  <  k2   .     Thus   the  limiting  behavior  as      |a|   -*•  °°     in  the  region 
of  regularity  of  the  spatial   transforms  is  determined  by  the  behavior 
of  the  space  domain  function  in  the   limit  as  the  origin  is   approached. 
Using  the  gamma  function   [18,  p. 183] 


r(x)  =   f  e_t  tx-1dt     ;     Re  x  >  0 


and  contour  integration,  we  find 


sn 


zn  enaz 


fca  +  J_rj£jleif<^ 


v^Ttt     a 


nTT 


and 


L%'^2  =  .i^;^" 


v^tt     a 


with     -1   <  n  <  0     in  both  integrals.     The  lower  extremity  on  the 
range  of    n     is  set  at     -1     to  ensure  the  existence  of  the  gamma 


function 


1 


and 


From  the  edge  condition  we  have     n=-  j    and  determine 

U+/       \  -1/2 

W   (a,co)     ^     a 

H+(a,u))     ^     a]/Z  as      |a|    *«     with     x  >  -kp 
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From  the  definition  of    S~(a),  it  can  be  seen  that  it  is  the 

spatial  transform  on  negative  values  of    z     of  the  current  normal  to 
the  edge  at     z  =  0     on  either  the  inner  or  outer  conductor,  exclusive 
of  that  current  contributed  by  the  incident  field.     The  magnetic  in- 
tensity is  continuous  at  all   values  of    p     for     z  =  0   ,  thus  the 
current  in  each  conductor  must  approach  zero  as     z  ->  0"  .     This  is 
apparent  from  the  fact  that  the  current  in  each  conductor  is  propor- 
tional  to  the  discontinuity  in  the  magnetic  intensity  across  the  con- 
ductor.      The  discontinuity  disappears  at     z  =  0   .     Because  the  total 
current  approaches  zero  as     z  ■*  0"   ,  the  quantity  which  transforms  to 
S~(a)     must  approach  that  of  a  TEM  mode  field  intensity  traveling  in 
the  minus  z  direction.     This  quantity  should  just  cancel   the  incident 
mode  field  intensity  at     z  =  0   .     Thus 


2Tra[H.   (a.z.w)   -  i|>.   (a,z,u>)]  =  2Trb[H     (b,z,a>)  -  i>     (b,z,u)] 

=s>     M(u>)  e"lkz       as     z  ■*■  0" 


The  transform  of  this  quantity  is 

0 
M(w)  e  e        dz  =  *— *•  (— -r-) 


2TrS"(a)    ^ 


i^TT       J 


/Ztt 


From  this   result  we  conclude  that     S"(a)  ^  ( — r-)     as     |a|   ->  °°    with 


t  <  +k2   . 


and 


Also, 


and 
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Collecting  some  results   from  Appendix  D,  we  write 

x+(a)     -v.     |a|1/2 

+  1/2 

y   (a)     %     |a|  as      |a|   -*°°    with     t  >  -k„ 


x"(a)     *     |a|1/2 


1/2 
y  (a)     ^     |a|  as     | ex]   ->  °°    with     t  <  k? 


An  examination  of  (2.40)   and  (2.41)  in  the   light  of  the 
asymptotic  behavior  of  the  terms  shows  that  each  side  of  these  ex- 
pressions goes  to  zero  as     |a|  -*■  °°    in  their  half  plane  of  regularity. 
Liouville's  theorem  is  therefore  applicable  and  each  side  of  (2.40) 
and  (2.41)  must  be  identically  zero.     The  quantities  of  interest  are 


W^g.Q))   x+(ct)        .    IZ    k  M(u))  _  Q 

~^^  ^  (otfk)  x+(k)  ^ 


and 


H+(g,g))  y+(a)       i    fl    k  M(a>)  =  Q 

b^  ^  («Hc)  y+(k) 


Therefore 


W+(a,.)     =     i  M    +ka  "M  (2.44) 

V7T  x+(a)   x+(k) 

u+(       \         i    /?    kb  M(cj)  /9  ,c\ 

^*y+(a)  y  (k) 
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and  the  coefficients     A(a,u)),  B(a,w),   C(a,co)     and     D(a,w)     have  been 
determined,  see  equations   (2.23),   (2.26),   (2.36),  and  (2.37).       The 
spatial   domain  field  intensities  are  now  written: 

u     /           \           ,•    /2       k  M(o))         K1(yp)  l0   .,, 

H     (p.a.w)   =   -  l    /- *-* t 1 —  (2.46) 

♦l  ^     YK0(Yb)     y+(k)  y+(a) 

u     I           \  ■      i  I*      k  MM         IT(yp)  r?  /m 

H     (p,a,a>)   -     1  /-    »-J r r (2.47) 

*3  v/  u     yl  (Ya)     x+(k)   x+(a) 


and 


,     /  \   _     ,•    /2  kM(o))    r 1 -, 

♦^(p.a.»)  -     ij-    —  [I0(Yb)K0(Ya)  -  i0(Ya)K0(Yb)] 


I1(YP)K0(ya)  +  K^YpJl^Ya)       I^YPjK^Yb) +K1(YP)I()(Yb) 

y  (k)  y  (a)  x  (k)   x  (a) 

(2.48) 
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3.  Frequency  Domain  Description  of  Fields  Interior  to  the  Coaxial 

Waveguide 

In  this   chapter  we  shall   derive  the  frequency  domain  fields 
interior  to  the  coaxial  waveguide.     We  will   completely  characterize 
the  TEM  mode  as  well  as  the     E         (q=l  ,2,3,- •  • )  mode  fields.     As  pre- 
viously  stated  no  other  modes  are  allowed  due  to  the  symmetry  of  the 
structure. 

Let  us  begin  by  recalling  equation   (2.11) 


H     (p,z,to)     =     H\  +  iK    (p.z.oj) 

q>2  cp  ^2 


The  incident  field  is  known  to  be 

M/    x     ikz 
H1     =     M(^)  e 

*  P 


Our  task  is  to  find    \\iA  (p.z.w)   .     From  the  definition  of  the  spatial 
transform  we  write 


\ 


ii     (p,z,co)     =     


00 

4>.    (p,a,w)  e  da  (3.1) 

cj>2 


Before  attempting  this  integration  we  will  examine  the  integrand 

and  all  singularities  will  be  located.     To  do  this  we  make  use  of  the 

relations   [20,  p. 375]  between  the  modified  and  ordinary  Bessel  func- 
tions that  follow: 

1  1     , 

-  j  vth  2"  irl 

I   (z)  =  e  Jv(z  e         )    ;  (-ir  <  arg  z  <  j ) 

Kv(z)  =  i  ri  e^  H^Hze^       );     (-tt  <  arg  z  <  J) 
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We  also  make  use  of  the  analytic  continuations   [24,  p. 80] 

I   (z  emi)  =  emv7Ti   I   (z) 
v  v 

K  (z  emi)  =  e"mv7Ti   K  (z)   -     i  *™J™  I   (z) 

VV  '  Vv     '  Sin    V7T        v^     ' 

By  these  formulas  it  can  be  shown  that 

Io(Yb)Ko(Ya)  -  I0(Ya)Ko(Yb)  =  - J{Jo(Y,b)Y0(T'a)-Jo(Y,a)Yo(Y,b)] 

n    2 

where     y  =  -\  ^jk  -  a    =  -iy     • 

It  is   clear  that  expression   (2.48)  does  not  have  a  branch  point 
at    a  =  +k   ,  since  we  may  replace    Y     by     -Y     and  not  change  the  value 
of  the  expression.     However,   (2.48)  has  a  branch  singularity  at 
a  =  -k     due  to  the  presence  of     x  (a)     and    y  (a)    .     From  the  above, 
we  conclude  that  the  only  singularities  of    i>     (p,a,o))     which  exist  in 
the  upper  half  plane  are  poles,   located  at    a  =  k     and     a  =     a 
(q=l,2,3,---)  with 


^'J*-iy'„f-*K-* 


where    Y'       is   the  qth  ordered  root  of 
1  oq  n 

Jo(Y'b)Yo(Y'a)   -  Jo(Y'a)Yo(Y'b)     =     0  (3.2) 

The     Y'       determined  by   (3.2)  are  the  eigenvalues  determining  the  cut 
off  frequencies  of  the     E        modes  of  a  coaxial  waveguide  [22]. 

To  evaluate  the  integral   (3.1),  we  integrate  along  the  closed 
contour  in  the  upper  half  plane  as  shown  in  Figure  3-1.     On  this 
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a  plane 


Figure  3-1.       Integration  Contour 
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contour  the  integral   converges  for     z  <  0     and  the  contribution  along 

the  semicircular  path  at     |a|  -*■  °°     is  easily  seen  to  be  zero  by  re- 
placing the  modified  Bessel   functions  in     ty     (p,a,w)     by  their 

TO 

asymptotic  values  and  using  Jordan's   lemma  to  evaluate  the  resultant 
integral   over  the  portion  of  the  contour  at     j otj   -*•  °°. 
By  the  residue  theorem,  we  have 

\l>     (p,z,co)  =  2tti'     I     residues  of  { if).    (p,ct,w)e~1aZ} 

at  the  poles  interior  to  contour 

of  Figure   3-1  (3.3) 

The  residue  of  the  pole  at     k     is 


h.<     ..irikM((D)e"1az-,r 1 i 

l«(o.k)[-JLJ ][  M   _  1  ] 


■laz 
a+k  ~o"~'"oxl~'       ~ov,~'"o' 


I^YPjKotYaJ+K^YpJI^ya)       I^YP^Yb)  +  K1(yp)I0(y>>) 

y+(k)  y+(a)  x+(k)  x+(a) 


=  ik  M(a))e"ikz       1im[j^EE] 
77  a-*  k  »/a+F 

[y+(k)]'2(^  Jin  ya+k)   -  [x+(k)]~2(-  ^nYb  +  Lg 
An  b/a 

■  srW.^"""2  -  r*+(k)r2)  ^  0.4) 


where  limiting  forms  of  the  Bessel  functions  have  been  used  to  charac- 
terize their  behavior  for  small   arguments. 
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The  residue  of  the  pole  at    a       is 

q 


lim  (a-a  ) 


ikM(h))  e 


-naz 


a-*- a. 


ivy{-  ^0(y'bU0(y'a)  -  JQ(y'^0(y'b))2} 


I^YpJ^fya)  +  K1(Yp)J0(Ya)       I^ypjK^yb)  +  K^yp)  JQ(yb) 


y  (k)  y  (a) 


x  (k)   x   (a) 


Using   the  L'H6pital      rule  this  becomes 
11M  [.  l^Me -, 


a-*  ol 


IT 


i^C 


!rJl(Y,p)Y0(Y.,a)  -^(y'pJJ^y'a)       J^y'p)  Y^y'b)  -  Y^y'pJJ^y'b) 


y  (k)  y  (a) 


+         + 
x  (k)   x  (a) 


] 


^[Y'(-f)(J0(Y'b)Yo(y-a).J0(y'a)Y0(y'b))] 


(3.5) 


n — i 

Observe  that    y*  =    i/k  -  a       and 


y; 


=  v k  -  a     .     By  invoking  relation 
oq  q 

(3.2)    the  derivative  in  the  denominator  of  the  above   relation  is   readily 


reduced  to 
d 


^{Y'[-f][J0(Y'b)Y0(Y'a)   -  Jo(Y'a)Y0(y'b)]} 

=  -fY'L^(J0(Y,b)Y0(Y,a)   -  J0(y'a)Y0(y'b))] 


Jy,[^r][b(J1(y,b)Y0(y,a)  -  J^y'ajY^y'b)) 


+  a(J0(Y'b)Y1(Y'a)   -  J^y'ajY^y'b))] 
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Putting  this  value  in   (3.5)  yields  the  residue 


-izVk2-  (Y'    )2 
i     k  M(<o)  e  Y«T 


J   (y*    b) 
o   'oq   ' 


o   'oq 


y+(k)y+(/k2-  (y'q)2)     x+(k)x+(Vk2-  (Y;q)2) 


[b<  Jl  Kqb>  YoKqa)  "  W>  Yl  (Ti,b))  +  a( Jo(Y;qb)  Y]  (YQqa)  -  J^a)  Yo(y-qb))j 

(3.6) 

J   (y'   b)        Y  (y«   b) 
The  relation     j°/'?qa\  =    vV^a)       was  also  used  to  simplify  (3.6).   This 
oVYoq   '         oVToq  ' 

relationship  follows  from  (3.2). 

iji.    (p.z.w)     for     z  _<  0     is   found  by  substituting  (3.4)   and  (3.6) 
into  (3.3).     The  result  is 

•iz/k*-(Y'    )Z 


R  _-ikz  ^'oq' 

z  <  0 


i,     (p.z.w)   =  -M(u)  Je"M  -     I     MU)Rf(p)e 
<P2  P  q=1  q  q 


R  =  jAj^ii/My2  -  [x+(k)r2]  (3.7) 


Jin  b/a 


Rq  \/^7 
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2k 


J   (y'   b) 
owoq   ' 


y+(k)y+(7k2-   (y^)""  Wga>    x+(k)   x+(  ^2.(^)2) 

(3.8) 
fq(p)  =  J1(y;qP)Y0(Y;qa)  -Y^Y^pJJ^a)  (3.9) 


The  magnetic  intensity  for  the  interior  of  the  feedline  is  seen  to  be 

H.    (p.z.oj)  =  ^^ +  iK   (p.z.u)      ,       z<0  (3.10) 

h>P  P  h>2 
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4.     The  Radiated  Electromagnetic  Field 

To  find  the  magnetic  intensity  in  region  1,  we  must  take  the 
inverse  spatial   Fourier  transform  of  the  spatial    frequency  domain  mag- 
netic intensity  found  in  Chapter  2.     The  integral   to  be  evaluated  is 

.     7       k  M(w)   K1(Yp)e-iaz 

H.   (p.z.w)  -  -  £         ^ t da  (4.1) 

*1  *  1       YK0(yb)  y+(k)  y+(a) 

for    p  _>  b   .     This  equation  was  obtained  from  equation   (2.46). 

Examination  of  the  integrand  of  (4.1)  using  the  analytic  con- 
tinuations of  the  Bessel   functions  and  the  value  of    y     on  both  sides 
of  the  branch  cut  reveals  that  the  integrand  has  branch  points  at 
a=±k.     Any  attempt  to  evaluate  the  above  integral   in  closed  form  by 
contour  integration  would  have  to  account  for  integration  along  a  path 
alongside  the  branch  cut.     If  the  contour  is   closed  in  the  upper  half 
plane  (for    z  <  0)     then  the  contribution  from  the  path  alongside  the 
branch  cut  in  the  upper  half  plane  must  be  calculated.     If  the  contour 
is  closed  in  the  lower  half  plane  (z  >_  0)  ,  then  the  contribution  from 
the     path  alongside  the  branch  cut  in  the  lower  half  plane  must  be 
calculated.     In  the  special   case  of    z  =  0     the  integration  around  the 
branch  cut  in  the  lower  half  plane  resulting  from  closing  the  contour 
in  a  negative  sense  equals  the  integration  around  the  branch  cut  in 
the  upper  half  plane  which  results  from  closing  the  contour  in  a  posi- 
tive sense.     There  are  no  poles   located  inside  the  contours   under 
consideration. 
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The  integrations  around  the  branch  cuts  are  quite  complicated 

and     are  not    evaluated  in  closed  form.     Since  our  primary  interest  is 
in  the  radiated  fields  (i.e.,  those  with  1/r  behavior), it  is  sufficient 
for  our  purposes  to  approximate  the  value  of  the  integral    (4.1)   for 
y  -*■  oo    9 


-jr-3 


Q     +   Z' 


In  the   limit  as     r  -»- °°  ,  we  replace  the  Bessel    function 
dependent  on     p       by  its     asymptotic  value  for  large  arguments  [20, 
p. 378].     Using  Figure  4-1,  we  readily  have 

p  =  r  cos  x  and  z  =  r  sin  x 

where     |x|    <.  tt/2     .       Using  the  asymptotic  value 

Myp)    *    I?-    e'YP 
1TK  V2yp 

the  integral    (4.1)  becomes 

-(yr  cos  x+  ictr  sin  x) 


ha  (r.e.o,)  ■  -iM("|k 

H  v^ffy  (k)  _l     YK0(yb)y  (a)[/yr  cos  x] 

as     r  ■*■  °°  ,     p  >.  b 


da 


(4.2) 


The  variable     9    has  been  introduced  at  this  point  since  this 
is  the  variable  in  which  we  wish  to  have  the  final  expressions.     It 
is  obvious   from  Figure  4-1   that 

e    ■    J-  x 
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(r,0,*>) 


Figure  4-1.       Coordinate  Axes 
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Define  a  change  of  variable  by  introducing  the  transformation 

a  =     k  sin  z  (4.3) 

z     is  the  complex  variable     x+iy   .     Under  this  change,  the  function    Y 
becomes 


Y 


=   \/k2sin2z  -  k2  =  -ik  cos  z  (4.3.1) 


It  is  necessary  to  specify  precisely  the  value  of    z  ,  since  it 
is  given  by  the  multivalued  function 

z  =  arc  sin(j-) 

which  has  infinitely  many  branches.     The  precise  specification  will 
allow  a  one  to  one  mapping  from  the  complex     a    plane  to  the  complex 
z  plane.     Since 

„iz      -iz 

g  =    sin  z  =     e     "  e 


k         °,M  '  2i 

we  may  write  the  equation 

e2iz  _  2i   «eiz  -   1     =     0 

This  quadratic  is   readily  solved  for     z 

z  =  -i  ln[1  £+  [l-(f)2]1/2] 

In  this  solution,  the  principal   branch  of  the  logarithm  function  has 
been   used  and  the  plus  sign  has  been  used  with  the  square  root  function 
[25],       By  the  branch  previously  defined  for    y  ♦  tne  above  equation 
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simplifies  to 

z  =  -i   £n[i  ^^-]  (4.4) 

The  equation  just  developed  is  useful  in  determining  how  the  two- 
sheeted    a    plane  maps  into  the  z  plane. 

Under  the  transformation  (4.4)   the  numbered  regions  of  Figure 
4-2  map  into  the  correspondingly  numbered  regions  of  Figure  4-3.     The 
path  of  integration  for  the  integral   of  (4.2),   depicted  in  Figure  4-2, 
maps  into  the  path  shown  in  Figure  4-3. 

As  examples  consider  two  points  on  the  path  of  integration  in 
the     a    plane 

1)  a  ->  -°°    with     |x|   ■+  0 

2)  a  -*•  +«>    with     |  t |   ■*  0" 


For  case 

t_  r/    i_i  ,  :i_i\2     .2-.1/2  _     „_  ri    i/i         k2 

71 


2 

Y  =  11m  [(-kl+i|x|)2-  k2]1/2  -     lim  [|a|(l   -  -*-_)  -  1|t|] 
a  ■*■«>  a-»-°°  2|a| 

\t\+0+  |t|+0+ 


and 


.i     lim    &n[1(-|ql+lr1|T|  +T)>  -1'    lim  [>  TTtM  --J  +  1« 

M-o+ 


For  case  2: 


2 

Y=l1m    [(|o|  -  i  | t|  )2-  k2]1/2    -iib[|o|(1  --5-y)  -1|t|] 
|a|-x»  |a|-»*»  2|a| 

|t|-»0"  |t|-^  0- 


and 
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(I) 


LY 


(4) 


Re  y  >0 

E3  Imy  >0 
I     I  Im  y  <  0 


Re  y  <  0 


Figure  4-2.     Integration  path  in    a    plane 
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■W 


A 


® 


© 


f 


© 


z  plane 


© 


© 


© 


7T 


Figure  4-3.       Integration  path  in     z     plane 
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2  =  _i     lim  £n[i(H   -  ;'Ti+^)]  =  -i.lim    £n[i  %L]  =  J.  .  1( 
|(t|-*« 
[t|+o" 


|a|-x»  *  |a|-Ko 


Upon  substituting  (4.3)  and  (4.3.1)   into  (4.2),  the  integral 
under  examination  becomes 


H.    (r.e.o.)    =    M!*UL- 
*1  ^F  y  (k) 

oo 

x     f  exp(ikr(cos  x  cos  z  -  sin  x  sin  z))  dz  ^4  g, 

_£  y  (k  sin  z)K  (-ikb  cos  z)   /-i rk  cos  x  cos  z 


The  integral   is  in  the  form  amenable  to  the  saddle  point  method  of 
integration,  since     r   ■*  °°    . 

Define  the     z     dependent  quantity  in  the  exponent  as 

g(z)  =  i(cos  x  cos  z  -  sin  x  sin  z) 

=  i   cos(z  +  x) 
It  follows  that 

jj    g(z)    =     -i  sin(z  +  x) 

The  saddle  points     z     are  the  points  at  which     g(z)     is  a  maximum. 
Therefore  from  the    above  we  know  that 


z     =  m  -  x  >  n=0>  ±2, ±4  ••• 
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Concentrating  on  the  point    z     =  -x  ,  we  note  that     g(z  )  =  i    . 

Since  a  steepest  descent  path  is   a  path  along  which  the  imaginary  part 

of    g(z)     is  constant,  the  steepest  descent  path  through     z     =  -x     is 

given  by 

Im  g(z)  =  1 

Since     g(z)  =  i    cos(x+iy  +  x)  =  i(cos(x+x)ch  y  -  i   sin(x+x)sh  y)    , 
the  steepest  descent  path  through     z     =  -x     is  given  by 

cos(x+x)ch  y  =  1 

This  path  is  shown  in  Figure  4-4.     We  are  assured  that  this  is  a  path 
of  steepest  descent  rather  than  one  of  steepest  ascent,  since  the 
Re  g(z)     monotonical  ly  decreases  along  this  path  as  we  move  away  from 
the  saddle  point. 

By  the  residue  theorem,  we  deform  the  path  of  integration  of 
Figure  4-3  into  the  steepest  descent  path.     The  major  contributions   to 
an  integral   along  a  path  of  steepest  descent  come     from  points  in  the 
neighborhood  of  the  saddle  point.     The  slowly  varying  terms  of  the 
integral  are  replaced  by  their  value  at  the  saddle  point  and  the  func- 
tion    g(z)     is   replaced  by  the  leading  terms  of  its  Taylor  series  about 
the  saddle  point. 


g(z)       *     i   -  j<z+x): 


M(u>)k  eikr   LDexp(-i  ^z+x)2)  dz 
H<,  (r,e,u>)  - 


JzH  y  (k)y  (-k  sin  xV-irk  cos  x  K  (-i kb  cos  x) 

(4.5) 
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Figure  4-4.     Steepest  descent  path 
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Along  the  steepest  descent  path 


and 


ch  y  = 


dy 

dx 


1 

cos(x+x) 


tan(x+x) 


cos(x+x)sh  y 


fi 


As  we  approach  the  saddle  point  in  the  manner     x+-x     and 
y  -*■  0     ,   the  slope  is 


11m    *■ 

dx 

y+o+ 


-i 


As  we  approach  the  saddle  point  in  the  manner    x  ->-x     and 
y  ■+  0"   ,  the  slope  is 

ii.  £  -  -i 

x+-X 

y  +0- 
Therefore  the  angle  that  the  SDP  makes  with  the  real   axis   at  the 
saddlepoint  is  either 


—      or      -  - 
4        or  4 


Examination  of  Figure  4-4  shows  that  at  the  saddle  point  the  SDP  and 
the  negative  real   axis   form  the  angle     3tt/4     in   the  second  quadrant 
and     -tt/4     in  the  third  quadrant.     This  observation  permits   us  to 
express  < 

z  +  x  =  u  =    ' 


u|   e 


i   3tt/4 


u  e   SDP     and     y   >  0 


u|    e"1   ^      '-       u  e   SDP     and    y  -  0 
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It  is  also  observed  that 


du      dlul     -1  tt/4        .    ..  ..,         .    . 

-rr  =     i.  '   e         '         at  the  saddle  point 


Since  equal   contributions   to  the  integral   of  (4.5)  are  received  as  we 
move  away  from  the  saddle  point  in  the  two  directions  of  increasing     u, 
we  write  the  integral  which  appears  in   (4.5)   as 


exp(-i  ^(z+x)2)dz    -    2 


kr,    ,2         _.  tt 
e      c  e       *     du 


SDP 


(4.6) 


Using  the  well   known  integral 


J-  2  r 

e  ax        2a 


the  value  of  (4.6)  is  seen  to  be 


2     (e 


^|u2| 


•     TT  .     TT 

4N   ,,..,       Q_1   4     /2T 
"  e  Vk7 


1  d|u| 


(4.7) 


Equation  (4.5)  becomes 


ikr 


H^  (r.e.w)   =  M(o>)  ^r 


i-T— r 


y  (k)y  (-k  sin  x)KQ(-ikb  cos  x)cos  x 


-}       (4.8) 


To  simplify  the  above  expression  we  use  the  identity  relating  modified 
Bessel    functions  with  ordinary  Bessel    functions: 


K  (-ikb  cos  x)  =  i  yHV'fkb  cos  x) 

O  c      0 
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From  Appendix  D,   the  expression 


■1 


y(a)  =[i  £  H^CKbJCJ^Kb)   -  JQ(Ka)]J 

is   formed,  where     K  =  +iy   .     In  expression   (4.8)  we  have 

i  jH^(kb  cos  x)  =  1  J-H^tkb  sin  e) 

If  in  the  expression     y(a)    ,  K     is  set  equal   to     k  sin  0     and     a     is 

2         2 
set  equal   to     k  cos   6   ,  which   follows   from  the  fact  that    y     =-K     = 

2       2 
a  -  k     ,   then 

y(k  cos   0)  =  y   (k  cos   9)  y~(k  cos   6) 

=  i  |  H^(kb  sin  0)[JQ(kb  sin  6)  -  JQ(ka  sin  e)] 

When  this  equation  is  solved  for     i  S-  hr   '(kb  sin  9)     and  the  result 
^  2     o 

substituted  into  (4.8),  the  far  zone  magnetic  intensity  is  seen  to  be 

ikr      J  (kb  sin  9)   -  J  (ka  sin  9)       +,.  , 

H.    (r.O.o.)  -  M(U)  *—    C~ iiT^V 3(y  {k+  C°S   9)) 

^1  r  sin   9  y   (k) 

(4.10) 

The  relation   (4.10)  has  the  restriction  that     p  _>  b   .     In  the 
region     z  <  0     the  equation   (4.10)   is  only  valid  in  the  region     p  _>  b; 
however,   for     z  >  0   ,  the     H         can  be  analytically  continued  from 
region   1   to  2  and  from  region  2  to  3.     Therefore  the  restriction  re- 
quiring    p  >  b     can  be  removed  in  the  region     z  >_  0     and  equation 
(4.10)  becomes   an  expression   for  the  magnetic  intensity  for  all   values 
of    p     in  the  far  zone   (r  ->  <*>)  exterior  to  the  feedline. 
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It  is   recognized  that  equation   (4.10),  except  for  the  factor 

y  *-r  C0S — -  ,  exactly  equals  the  results  obtained  from  assuming  that 

y  (k) 
the  incident  field  is  the  total   aperture  field  and  that  no  contribution 

is  made  to  the  radiated  fields  by  currents  which  flow  on  the  waveguide 

structure. 

The  new  field  factor  contribution     ^  \  cos — '-    to  the  coaxial 

y  (k) 

waveguide  antenna  first  obtained  in  this  work  modifies  the   field  pat- 
tern  obtained  by  the   conventional   solution   to  the  problem  of  radiation 
from  this  type  of  antenna  [19,p.302].       It  will  be  seen  later  that  this 
factor  has   a  significant  effect  at  high  frequencies   and  modifies  the 
field  pattern  to  direct  the  radiation  in  a  more  forward  direction. 

By  equation   (2.11)  we  expressed  the  magnetic  intensity  in 
region  2  for  all     z     as  the  sum  of  the  incident  field  and  another  term 
which   represented  the  reflected  and  higher  order  mode  field  intensities. 
By  the   radiation  condition,  we  know  that  the  incident  field  cannot 
exist  in  region  2  for     z  >  0   . 

To  show  that  the  incident  field  intensity  is  canceled  in   region 
2  for     z  >  0   ,  we  evaluate  the  integral 


oo 


*.   (p,z,w)  =  —  I  f .    (p.a.u)  e"laz  da  ,     z  >  0  (4.11) 

— oo 

by  using  the   contour  of  Figure  4-5,     iJj.    (p,a,u)     is   given  by   (2.48). 

<j>2 

The  value  of  the  contour  integral   is   calculated  by  the  residue  theorem. 
The  value  of  the   integral    (4.11)   is   the  negative  of  the   contribution  from 
the  path  around  the  branch  cut  and     -2-rri      times  the  summation  of  the 
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Figure  4-5.        Integration   contour 
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residues  at  the  poles  of  the  integrand  located  on  the  interior  of  the 

contour     c.     The  contribution  along  the  semicircular  arc  at     |a|   ■+»  is 

zero. 

The  integral    (4.11)   is  not  evaluated  in  closed  form  due  to  the 

complex  integration  on  the  path  alongside  the  branch  cut.     However,  we 

shall  show  that  the  contribution  to  the  value  of  the  integral    from  the 

circular  path  around  the  branch  point  (k)  is  just  enough  to  cancel   the 

incident  field  intensity. 

i  fi 
On  the  circular  path  around  the  branch  point    a  =  -k  +  te       , 

e  ■+  0     and    0  <  8  <_  2tt  .     The  integral   of  interest  becomes 


lim 
e  +  0 


L      JH 


fc 


(p,a,w)  e  da 


] 


-k 


(    auf    \  Jkz    /     -Yp  in  ya  +  —       -yp  in  yb+  — 
-     lim      _lkM(w)  e         f     r_  yp  yot  da 

c-O^tt     in  | 


1P_]  aa 
J        yT(k)  yT(a)  xT(k)   xT(a)  * 


lim 
e+O 


i-i 


k  M(h))  e 


ikz 


it  in 


■vm 


2tt  ^ 

1      -,      {  da      ■ 

^y]   J  X  / 

o    Y  p  ; 


M(hi) 


p    £il¥ 


F  [yTkT  "  xTFJ 


•]eikz 


We  have  made  use  of  the   limiting  forms  of  Bessel   functions   for    small 
arguments   and  the  following: 
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x+(-k)   =  x"(k)        ;         x+(k)   x"(k)   =   x(k) 

y+(-k)  =  y~(k)        ;         y+(k)  y"(k)   =  y(k) 

To  see  that  this  is  just  enough  to  cancel   the  incident  magnetic  inten- 
sity, note  that 


1  1       _  0n  b 

yTFT  "  3^kT  '  £n  i" 


which  is  easily  found  from  the  definition  of    y(k)     and     x(k)     given 
in  Chapter  2  and  the   limiting  forms  of  Bessel   functions   for  small 
arguments. 
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5.     Time  Behavior  of  the  Radiated  and  Reflected  Electromagnetic  Fields 

The  time  behavior  of  the  field  intensities     is    found  by  comput- 
ing the  inverse  time  transform  of  the  frequency  domain  expressions.   To 
find  the  time  behavior  of  the  radiated  magnetic  intensity  emanating 
from  the  waveguide  antenna,  the  inverse  transform  of  equation   (4.10) 
must  be  obtained.     Because  of  the  complicated  functional  behavior  of 

the  factor    y  I     cos   $>     the  computation  of  the  integral 
y  (k) 

i       7            Jkr    J  (kb  s""n  9)-J„(ka  sin  6) 
h  (r.e.t)  =  -^   j  M(U)  V  t- liFT 3 


y+(k  cos  e)  e-1o,t  daj  (5J) 

y  00 


is  extremely  difficult 


For  harmonic  time  dependence  the  exciting  voltage     V(t)     is 
ico  t 
Re(A  e         )   (see  equation   (2.7)).     Then 


M(w)  =  Re[^J^^6(oJ-.0)] 


where     6(x)     is  the  Di  rac  delta  function.     Therefore,  the  case  of 
harmonic  time  dependence  (5.1)   reduces  to 


,+(?cose)    „-'K(t-r/c)] 

H   (r,9,t)   =  Re 


y^o  £n  b/a    y+ 


b/a    y^(^) 


wo\  r  sin 


X 


a)  b  co  a 

[Jo(-|_sin  e)  -  J0(-f-sin  e)]J  (5.2) 


+  /    0 
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The  value  of    k     is     u>/c     as  given  by  (2.4.3) 

The  factors  of  equation   (5.2)  except      y  ^~c~  cos     j     are  we]- 

understood,  and  are  not  in  need  of  further  simplification. 
From  Appendix  D  we  transpose  the  function 

k 


y  (k  cos  6)  _ 
y+(k) 


n  + 


n 

m=l 


m 


I02'(u)  £n 


k  cos  9 

k 
vk  -  u 


ik  cos  6 

I  H n- 


,  4^1 -cos   e)(i^iisin  6 
2m  '  fi    tt 


1  + 


'1  + 


a+b 


k(|g)(cos  6-  l)(l-C+l  J) 


exp(k(f^cos  9[-£n(k(|^-)cos  9)  +  |^-£n(2  cos  6) 

.  k  (a^.Ji^)^b 

i2tt  tt        a+b  ' 


(5.3) 


In  this   formula  we  have  used  the  fact  that     if      a  =  k  cos  6   ,   then 


§*((+  HI) 

VY-  (a-k)' 


o    r-ik  sin  9  +  k(cos 
*nL-ik  sin  9  -  k(cos 


^j]  ■  H0]  -  +* 


using  the  principal  branch  of  the  logarithm  function.     Additionally, 


K(»{«).i[i    *  (_ — '   — ,1 

*L        *ub    j2(ub)+y2(ub)  J 
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+  1/2 

It  is   also  known  from  Appendix  D  that    y  (a)  *x»  |a|  as 

a  ■+  °°  .     From  this   asymptotic  behavior  we  see  for  very  high  frequen- 
cies  (A  «  b-a)   that 

lim      y+(k  cos  9)     %   |cos   0|l/2  (5<4) 

k  ■*  °°       y  (k) 


In  Chapter  2,   there  was  imposed  a  restriction  that     kb.ka  «  1 
for  non-harmonic  time  dependence.     In  keeping  with  this   restriction, 
equation   (5.4)  should  only  be  used  for  the  case  of  harmonic  time 

dependent  excitation.     For  values  of     ka.kb  «  1     the  value  of 

y+(k  cos   9) 

tt-t is  more  complicated  but  instructive  to  determine.       To 

determine  the  value  of  this  field  factor  in  the  low  frequency  limit 

(X  >>  b-a)     we  use  the  series  expansion  of  the  infinite  product 

1  +  W  )  i-^^d-  cos  e) 
,  '   m1     I        2m 

m-i  1771  cos  ll 

m=  1    l  \  +  t\ — i —  I 
1   m1        ; 

.  ,  .  |«i*)(1-a»  6)    J  [_1         .L]  ♦  0(ka)2 

"■1  i<v,K-z-> 


The  infinite  series  is  the  difference  of  two  divergent  series  but  is 
itself  convergent  [26,  p. 33].     This  is  apparent  since 

1  J_      „  1  1        _       1 

i      i    a+t 
1    m1      2 


(m  -  j)tt  4um 


1 
as     m  ■*  °°     and  the  infinite  series       [      -y      is  known  to  be  absolutely 

m=l     m 
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convergent.     The  following  integral 


■     K(2)(u)   £n 
6-K)+ 


1 

+     k 

A2+ 

2 
u 

1 

k  cos 

6 

s/^rj 


du 


also  appears  in     *■  ■■    c .     This  integral  will  now  be  shown  to  be 

y  (k) 

convergent  and  its  approximate  value  will  be  determined  in  the  limit 


kb,ka  «  1    . 

Upon  substitution  of  the  value  of     Kvt/(u)     and  simplifying, 
the  above  integral  becomes 


-(2), 


bn-^(-     ] 


6^L'--b^(ub)-Y2o(ub) 


)]  An 


\u    ]  -  cos  6 ldu 

L     cose+/i-(^)2  J 


Introducing  a  change  of  variable,     ub  =  z  ,  we  get 


6  ->0 


1    !       2_r  1  ] 

*  L         ™     J2(z)+Y2(z)  - 


in 


1  + 


1    -    COS    6 


cos  e  + 


A7^ 


2     J 


dz 

(5.5) 


To  examine  the  behavior  of  (5.5)   for  small   values  of     z   ,  we 


consider 


TT  J 


[1 


1 


«S->0 


*z  l+4?,n2z 

TT 


■]  Jin 


1  + 


1    -    COS    0 


COS  8  +  /l    +    (£r)2 


dz 


W 
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This  integral   converges  as         [-( » *— )]  dz     converges 

■  l+^-£n  z 


Now 


0 


f[l(-         ] 


6 


1+^2  ^n2z 


)]  dz 


dz 


TT 


2? 


i  z  £n^(z) 


d['  £iTT]    "  •  ~ 


4    £n  z 


This  quantity  is  seen  to  be  convergent  for  all     e  ^  1    .     We  can  also 
remove  the  principal   value  sign  and  change  the   lower  limit  to  zero 
without  affecting  the  value  of  the  integration. 

At  the  other  extreme  for  large  values  of    z   ,  we  consider 


1 

TT 


1-  2    r       ] 
tT7  l7TXI7 


H-'(z) 


]}  £n  fl  + 


1  -  cos 


cos 


1i% 


)2-! 


dz 


for    M  »  1    .     The  asymptotic  value 

H(1)(z)    *      (—)yZe  *  [1+0(1)]     for     z  + « 

0  TTZ  z 


is  now  substituted  with  the  result 


±      [1 


1 

TT 


[l+0(£)] 


O(^)  £n 


-]     £n 


1  + 


1  -  cos 


1  + 


cos  0-1 


1     -    COS    0 


/<p 


dz 


V  (fe)2  "  » 


-  1 
dz 
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The  logarithmic  function    iln(l+x)     can  be  expanded 


2         3 
An(l+x)  =  x  -  -^-+  4-  - 


,     Ixl   <  1   ;     x  1  -1 


When  this  expansion  is  considered  along  with   the  fact  that 


1  + 


1   -   cos   6 


=e>    1  + 


cos  0  -   i 


*JW 


as 


z  ■*  °°       (with     z     real    and     |A|    finite),  then 


1      0(1)  «n[  1  ♦ 


1     -    COS 


COS    6 


-'/ft)^ 


]dz 


1 

IT 


0(1)  (|)    dz 


M 


This  integral   converges. 

The  integrand  of  (5.5)   has   an  integrable  singularity  at     z = kb 

when     9  =  tt/2   .     By  the  above  we  conclude  that,  for  all   values  of     9 

and     k     the  integral   of  (5.5)   is   convergent. 

y+(k  cos   0) 
Since  our  interest  is  in   the  value  of      7~ TR1 ^or 

ka,kb  «  1    ,  we  next  endeavor  to  find  the  value  of  (5.5)   in  the  limit 

as     kb  +  0   .     We  will   use  this   limiting  value  as  an  approximation  for 

the  integral   appearing  in  equation   (5.3). 
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lim    1  f  [1 
kb  *°  *  6 


~(2      ]     2       )]Wu       ]  -  cos-9— 


dz 


kb 


■     "■    M     [l"|i<-5 ^—z )]*nfl  + 

kb  ^0  ^  i  uz  jf(z)+v?(z)  \ 


1  -  cos 


0(z)+Y;(z) 


cos  e  +  7i  -  (j^-)2 


dz 


+     lim    If  [1-^ 
kb+0  *  /L 


"J5(2)+Yt(z) 


1-5 )]£n(l  + 


1  -  cos  e 


cos  e  -     i 


/W^ 


dz 


=     0 


Note  that  the  quantity 


lim     [JLn(l  + 
kb+0 


1  -  cos 


cos  e  -i  J^)Z-  1 


-)] 


From  the  above  work  we  conclude  that  as     kb     goes  to  zero,  the  integral 
(5.5)   goes  to  zero.     Thus,   for  low  frequencies,  we  approximate  the 
function 


exp 


6-i+ 


K(2)(u)     in 


1  + 


[J- 2" 

vk  -  u 


1  + 


k  cos  e 


du 


vk  -   u  > 


appearing  in  equation   (5.3)  by  unity, 


If  in   (5.3)  we  make  use  of  the  limit  that     x£nx  =  x     as     x-*0, 

then  the  field  factor    *   (j[  cos  °) 

y  (k) 


can  be  approximated  as 
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^Mfji- (i  - ik(^)[1  - cos  9]  Vn^'-N(lil2) 

y   (k)  m=l    I  otm[   ~2~ 

x  exp  I  i  -^i  sin   6  +  ik(|^-)[C(l   -  cos   e)]-  fj^cos  e  s.n(2  cos  6)  -  £n  2]] 

for     ka.kb  «  1  (5.6) 

The  time  behavior  of  the  far  zone   radiated  magnetic  intensity 
for  a  non-harmonic  time  dependent  incident  field  intensity  with 
ka,kb  «  1     is  given  by  (5.1)  when  equation   (5.6)  is   used  for 

y  *     cos — -  .     In  the  far  zone  the  value  of     r     in  equation   (5.1)   is 

*  M  a+b 

very  much  larger  than     b     or    -g—  .     To  the  approximation  used,   the 

exponential   function  in   (5.6)   is   replaced  by  one  and  for  that  reason 

it  does  not  contribute  to  the  phase  factors  in  integral    (5.1).      In 

addition  the  term  of  (5.6)   containing  the  infinite  sum  can  be  dropped. 

This  is  apparent  since     ka,kb  <<  1     and  the  Fourier  transform  of  this  term 

will  be  of      order     l/c  ,  which  will  make  its  contribution  wery  much 

smaller  than  the  contribution  from  the  remaining  term.     In  view  of  this 

y  (k  cos  e)     a    1 

y+(k) 
and 


H  (r,9,t)  ■  — 


./   s   -i  (cut-  kr) 

!^^ . — 5 [J  (kb  sin  e)  -  J  (ka  sin  e)]  dw 

r  sin  0         L  o  o  ' J 


for    ka.kb  «  1  (5.7) 
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When  the  incident  magnetic  intensity  is  harmonically  time 
dependent,   then  the  value  of  (5.7)  is  easily  obtained.     In  this   case 


M(oi)  -  Re^2T  /—  — - —  5(w-a>J 

yo  In  b/a 


and  for  low  frequencies   (A  »  b-a)     the   radiated  fields  are 


co  a 


fe~     cos(o)  (t-  -))  go  b 

Vr'ejt)   ~~  ViT      r(ln°b/a)s1"  9  ^oHf"  sin  6)  "  JoHh  sin  °)] 


Ee(r,9,t)     =     y^H^(r,6,t) 


(5.8) 


For  non-harmonic  time  dependence  the  radiated  field  can  be 
found  by  direct  integration  of  (5.7)   or  by  a  time  domain  convolution. 
The  integrand  of  equation   (5.7)    can  be   considered  as   the  product  of 
the  Fourier  transforms  of  two  time  domain  functions.     Therefore  this 
integral   can  be  evaluated  by  the  convolution  theorem  and  the  impulse 
response  is   readily  identified  from  the  result. 

In  applying  the  convolution  theorem,   use  is  made  of  the  fol- 
lowing well   known  integral   [27]: 


(    JTR 


y^ 


i  2     3. 
Vz  -  t 


J   (az)e"         da  ={ 
o 


;     o  <  |t|  <  z 


t    =  z 


;       otherwise 


(5.8.1) 


Upon  evaluation  of  equation   (5.7)  by  the  convolution  theorem  we  get 
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th  e  result  which  follows: 


t' 


— oo 

where     f(x)     is 


(5.9) 


V£s1n  e)2-  t2 


f(x)  =< 


M  sin  e)2-  t2     /£sin  e)2-  t2 


b     . 

;     x  < sin 

c 


;  --  sin  6  <_  t  <_  -  -  sin  0 


-  sin  e  <_  t  <  jjr  sin 


■Msin  9) 


1 2" 

-    T 


a  b     . 

;     -  sin  9  <_  t  £  -  sin 


S     t  >-sm 


(5.9.1) 


and 


t'  =  t  -  - 

1  c 


The  response  to     an  impulse  of  voltage  excitation  is  therefore 


IS.    f(t  -  V 

H^r.G.t)   :      /       iTr  sin  e  Jin  b/a 


(5.10) 


From  the  character  of  the  impulse  response  we  see  that  at  an 
observation  point  in  the   far  zone,   the   radiation   lasts  a  total   time  of 
(2b/c)sin    8  .     By  ignoring   for  the  moment  the   frequency   restriction  of 
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ka,kb  «  1     and  allowing  an  excitation  with  zero  rise  time,   a  step 
function  of  voltage  .   then  by  integration  of  the  impulse  response  we 
obtain  the  step  response, 


Vr'e't)=JiT( 


l 


r) 


o    ttt  sin  6  Jin  — 

a 


J+  sin"1^ 
Z  b  sin 


.     -1  TC 

Sln  bW~l 


sin 


-1 


TC 


a  sin  6 


J+  sin-l  _^ 
2  b  sin  G 


where     t 


_r 
c 


;     t  _< sin  6 


sin   G  <  x  < sin   6 

—      —      c 


—  sin  6  <  t  <  -  sin 
c  —      —  c 


■I  sin  0  <  t  <  —  sin  9 
c  —      —  c 


t  >  —  sin 
—  c 


(5.11) 


It  is  advantageous   to  have  the  step  response  for  numerical   cal- 
culations since  its  magnitude  is   finite  at  all  points  and  the  far  zone 
fields   can  be  found  by  the  convolution  of  the  time  derivative  of  the 
excitation  and  the  step  response.     This  procedure  does  not  violate  the 
low  frequency  restriction. 

Because  of  the  low  frequency  restriction  the  input  excitation 
should  not  contain  any  high   frequency  components.     To   render  the  high 
frequency  components  in  the  excitation  negligible,   the   rise  time  of 
the  excitation  pulses  must  be   long  compared  with   the   time     b/c   .      If 
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the  input  is  of  the  form     [1   -  e"  '    ]  u(t)     then  the  time  constant     T 
must  be  chosen  so  that     T  >>  b/c   .     This  will  preclude  any  violation  of 
the  low  frequency  restriction. 

For  such  inputs,  equation   (5.9)   contains  convolution  integrals 
of  the  type 


1   -  e 


JtkL 

T 


-  ^sine  v^sin  e)2-  x2 


dx,  with sin  6  <  t  <  -  sin  9 

c  —      —  c 


and    t  <  t' 


-    t-T 


(1  -  e 


)   dx 


(1   -  e     T)   dx 


,-ifi/T  y 


x  dx 


b  I 


in  6  /(|  sin  e)Z-  tZ 


,(!-. "I  t'l/T) 


dx 


It'l/T 


;in  e   Mlin  e)d  -x2 


JWsUi  e)2-(t')2 


0<|t'|    <  ^-sin  e 
—  —  c 


Since     T  >>  b/c    we  neglect  the  second  term  and  conclude  that 
the   response  to  an  input  of  finite   rise  time  with  time  constant 
T  >>  b/c     approximately  equals   the  step   response  multiplied  by 


Exci  tation 
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Figure  5-1.     Radiated  magnetic  intensity  of  the  coaxial  waveguide 

antenna  for  (a)  step  function  and  (b)  exponential  excitation 

Ordinate:     H,  =  H.(/|i  /l     r  v,n  -  sin  e)   (volts) 
(p        (j)       o     o  a 

Time   constant  for  (b)     is     cT/(b  sin   e)   -    .1 


Antenna  Dimensions:      a 


8b 
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A  basic  assumption  used  in  Chapter  2  is  that  the  total   current 

in  the  inner  conductor  is  equal   and  opposite  to  that  in  the  outer  con- 
ductor.    This  assumption  implies   that  the  contribution   to  the   radiated 
fields  from  all  parts  of  the  structure  outside  of  the  aperture  are  neg- 
ligible.    Examination  of  the  impulse  response  shows  that  the  radiation 
appears   to  originate  from  the  edges  of  the  open  ended  structure 
located  in  the   aperture  plane.     The   outer  edge  gives   a  positive   con- 
tribution  and  the   inner  edge   gives   a  negative   contribution.     The 
radiation  starts  when  the   incident  current  reaches   the  edge   and  lasts 
a  time  proportional    to  the   time  it  takes   light  to  traverse   the  diameter 
of  the   conductor.      The   difference  in  sign   of  the   contributions   is   due 
to  the  difference  in  sign  of  the  current  on  the  conductors. 

We  now  turn  to  the  problem  of  determining  the  time  behavior  of 
the  field  intensities  for     a_<p£b;z_<0.     The  magnetic  intensity 
in  this   region  is  given  by 


H     (p.z.t)   = 


1 


/2T 


e-ia,t  M(u>) 


eikz       Re-ikz 


P  P 

-iz7k2-   (V)2l 
-     I     R   f  (p)   e  q 

q=l     q  q 


dw         (5.12) 


-iuj  t 


For  harmonic  time  dependent  excitation   voltage     Re(Ae     °   )    , 


M(w)   =   Re[/2T 


o       A 


o  lr\ 


b     6^-uo)] 


and  (5.12)   becomes 
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(T      .       cos[w  (t  -  £)]       cos[u>  (t  +  £)] 
p,z,t)  =  /-°  _£-[ 2 9 o c_£  R(     } 


+     J,    fq(P)    RqK)e 


0) 


^>2-  (viq)2 


(5.13) 


where     R(oi  ),   RQ(^0)   and     fQCp)     are  the  functions   (3.7),   (3.8)   and 
(3.9)     evaluated  at     w  =  go     ,  if  applicable. 


w      should  not  corres- 


pond to  the  cut  off  frequency  of  one  of  the     E      modes  of  the  coaxial 

waveguide,  since  our  formulation  is  not  valid  for  these  frequencies. 

We  also  note  that  as     ui    -*■  °°  ,  R(oj  )  -*■  0   ,  thus  the  TEM  mode   reflected 

o  o 

field  tends   to  vanish   for  very  high   frequency  excitations.       If 
ka.kb  «  1    ,  it  is  obvious  that  for  large  values  of     z     the  higher 

order     E        modes  are  negligible  since  they  are  evanescent  in  this 

on  3     3  "* 

frequency  range.  Since  the  time  behavior  of  the  incident  field  is 
given,  we  shall  only  be  concerned  with  evaluating  the  integral  for 
the   reflected  TEM  mode   fields   or 


H4Jp'z,t)  reflected 
c  TEM 


1 


^ 


\  MM   Re-ikz-ia)t  dw 


R-  [y+(k)]"2  -  [x+(k)] 

a 


-2 


(5.14) 


From  the  values  of    y   (k)     and     x  (k)     determined  in  Appendix  D,  it 

2  2         +         -2 

can  be   adduced  that  to  the   order  of  (ka)     or  (kb)    ,   [y   (k)]       =  0     and 
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[x+(k)]"2  -  -in  ^  (ka)]  {l  -  2ika    l  (-r-ir-  -  I-)} 

a     °  n=i    |aja       niT 


_  SJca.  (i   _  c  +  JLn  2)        ., 

x      e     17T  e"lka 


x     exp[2  I       K(1)(u)    {£n(1   +         lkl      )}du]  +  0(ka)2   (5.15) 
6+0+  7k2-  u2 


By  examination  of  the  integral   in   (5.15)   at  the  extremes  of  very  small 
and  very  large  values  of     u   ,   it  is  shown  that 


I       K(1)(u)    Un(l   +  LM— )}   du 

sV  J777 


is   convergent   for  all    values  of     k   .     Additionally,   the  integrable 
singularity  at     u  =  k     presents  no  difficulty  to  evaluation.     After 
substitution  and  a  change  of  variable  the  integral   becomes: 


1 

TT 


f  „     H  +  it»l<!  x)Jo<x'  -  Jl(!  «l».!«)l*»,WJ.(l  x)  -  Vx)Vo(?  *> 
[J0(|x)-J0(x)f  +[Y0(|x)-   Y^x)]2 


(x.n(l   +  1     — )}     dx  (5.16) 


^     "    & 


ka; 
For  small   values  of    x  >   integral    (5.16)   can  be  approximated  by 


p 
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4       br2     a   -,       1   b       2   „  2  1       2   .  b 

— -  -  H—  hTJ  -  w  —  x  —  In  x ^  x  —  £n  —  x 

ttx       a  tt    dx  J        2   a       TT TTX  2       it  a 


}dx 


(1 


-  i[x  Jin  x]£  +   1)   -  -£n  ^ 
ttl Jva  tt        a 

<*  !»2 


)     dx 


This  last  integral  is  finite,   thus  the  integral    (5.15)  is  convergent  as 
u  +  0   .       For  very  large  values  of    x  ,   (5.16)  is  proportional  to 


oo  4 

TTX 


1- 


+  [|sl"n(|x-  ^)cos(x-J)  +  0(1)  +  [additional  terms]]  1 


T7£ 


TT^ 


y[cos(|x-  J)-  cos(x-J).+  0(l)]S  [s1n(|x-  J)  -s1n(x- J) +  0(1)]'} 


x     {£n(l  + 


<V^F^ 


-)}  dx 


_4_ 

TTX 


[1    -T 
J  _£_  (2  +  0(-)) 

M  TTX     V  VX 


y— ]     {An(l  + 


i  AS^ 


:)}  ^ 


■J 


Lka 


[1 


_4_ 

TTX 


—  (1+oW) 

TTX  X" 


—]       (J)    dx 


where  we  have  used  the  expansion  of    £n(l  +  z)     with     z  =  A/x  ;  z  «  1; 
z  ?  -1    .     Additionally  the  convergence  of 


sin  mx  (— )   dx       and 


00 

t 


M 


cos  mx(-)   dx 


has  been  utilized.     We  are  now  able  to  conclude  that  the  integral 
in   (5.15)  is  convergent  as     u  ■*  ™  . 
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Combining  all   the  foregoing,   it  is   apparent  that  the  integral 

appearing  in  (5.15)   is   convergent  for  all   values  of    k     and  it  is  appar- 
ent from  taking  the  limit  as     ka  ■*  0     of  (5.16)   that 

oo 

lim     f  K(1)(u)   {£n(l  +       M      )}   du  =  0 


ka  +  0 


JJ~7 


Additionally,   the  term  involving  the  summation  in   (5.15)  is  very  much 
smaller  than  one  and  therefore  ignorable.     The  exponential   factors  are 

also  negligible  when  large  values  of     z     are  concerned.     Therefore 

+        -2 
[x  (k)]         in  the   low  frequency  limit  becomes 


[x+(k)]"2    =  .-  In  \  [J  (ka)]e"lka  for    ka,kb  «  1  (5.17) 

a      o 


The  reflected  TEM  mode  field  is  now  evaluated  by  substituting  (5.17) 
into  (5.14).     The  results  are 

oo 

u                i     -,  *\             If  nHoj)    ,   /k,\     -ik(z+a)-iwt  /c  10n 

H^  (p.z.t)   b -^— '-  JQ(ka)  e  du  (5.18) 

Reflected  TEM  ^  L> 

By  the  convolution  theorem  and  well  known  integral    (5.8.1) 

t' 
%  (p,Z,t)=-X-f     M(t'-T)    g(x)   dx  ;      z<0  (5.19) 

2reflected  TEM  -«> 


where 
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g(x)  =   \ 


W^ 


2      1  L 


ti-| 


0  ;       otherwise  ' 


and 


t  - 


C 


If  we  again  ignore  for  the  moment  the  frequency  restraints  and  allow  a 
unit  step  function  voltage  excitation,  then  the  reflected  field  is 
given  by  the  integral  of    g(r)    .     Carrying  out  the  integration  produces 
the  following: 


H,     p.z.t      = 

<t>2 

reflected  TEM 


§J 


f) 


o     irp  Jin  — 

a 


J.sin-'f 


TT 


—       c 


3   ^         ^    a      \ 

_  TT  <    X    <  —     / 
C  —      —  C 


T    >      J 

—      C 


^ 


T   =    t 


and  the  electric  field  is 


Ep  (p,z,t)  =  - 
Reflected  TEM     ° 


-    HA  (p.z.t) 

2 

reflected  TEM 


e_     ""(J). 


In  the  case  of  a  step  input  the  reflected  electric  field  is 

2a      |z| 
essentially  the  negative  of  the  incident  field  for     t  _>  —  +  1— L 


c      / 
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(see  equation   (2.6)).    The  current  reflection  coefficient  is  -J   (ka)  for 
ka  <<  1    ,  which  is   consistent  with  transmission  line  theory.     For 
allowable  inputs  of  the  form  (1   -  e"  '    )  with     T     large  enough  so  that 
the  low  frequency  restriction  is  not  violated,   the    reflected 
fields  are  given  by  the  above  equations, multiplied  by  (1   -  e  ), 

f  =  t-4--f  . 

c        c 


-96- 
6.     Radiation  from  the  Pulse  Excited  Cylindrical  Monopole  Antenna 

The  objective  of  this  chapter  is  to  examine  the  radiation  from 
a  pulse  excited  cylindrical  monopole  antenna  in  the  light  of  results 
obtained  from  the  preceding  study  of  a  coaxial  waveguide  antenna.     We 
have  determined  that  in  the  far  zone  the  radiation  appears  to  originate 
from  the  edges  of  the  open  ended  structure  located  in  the  aperture 
plane.     The  edge  starts  radiating  when  the  leading  edge  of  the  incident 
current  pulse  on  the  feed  line  impinges  upon  it.     Radiation  from  the 
inner  edge  is  opposite  in  sign  to  radiation  from  the  outer  edge. 

We  begin  by  considering  the  approximate  far  zone  fields  from 
the  annular  slot  antenna.     This  antenna  has  been  the  subject  of  several 
authors   [28-32].     The  structure  of  the  annular  slot  antenna  is  that  of 
a  coaxial   line  terminated  in  an  infinite  ground  plane  as  shown  in 
Figure  6-1.     Note  that  the  center  conductor  is  not  hollow. 

We  will   follow  the  approach  of  Ref.    [32]  and  replace  the  elec- 
tric field  in  the  aperture  by  an  equivalent  magnetic  current  and 
develop  the  fields  in  the  half  space     z  _>  0     from  this  equivalent 
source  and  image  theory.     The  equivalent  magnetic  current  is 


)     =  4    x  E'   a'     =  -E'   a! 

m  z         p       p  p     <J> 


Since  this  source  is  just  above  a  perfectly  conducting  ground  plane, 
the  electric  vector  potential   is  given  by 


-97- 


2a  2b 


Figure  6-1.     The  annular  slot  antenna 
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e       f    /■     2E'   a'   e 

^  "   ■  ■>    ->,, p'dP'd*' 

0     a  k-r'| 


where 


and 


e  E     =  -  V  x  F 
o 


r-r'l   =  [p2+  p'2  -  2pp'    cos(cD-(D')  +  (z-z')2]1/2 


Now  it  is  clear  that  the  electric  vector  potential   and  the  fields  are 
independent  of  the  coordinate     <j>  ,  from  the  cylindrical  symmetry  of 
the  structure.     Therefore  we  can  set     <j>  =  0    without  loss  of  generality, 
We  also  know  that  the  only  nonzero  component  of     F     is  the     $     com- 
ponent due  to  the  symmetry  of  the  equivalent  magnetic  source.     From 
geometry  we  find 


a,,   =  a     cos  4>'   -  a y  sin  <t>' 

<p  y  * 

since     &  =  0   ;  a     =  a,  for  this   formulation.     Therefore 

y       <f> 


£.         27  br  .ikR 


F  =  "  2?  a<j, 


E   ■    cos  <J>'  ^—  p'd^'dp' 


0     a 


with 


R  =   (p2+  p'2  -  2pp'    cos   (J)'   +   (z-z')2)1/2 


From  Maxwell's  equations  and  the  Lorentz  gauge  it  is  an  easy 
matter  to  show  that 
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H  ■  ^  [-  V(V-F)   +  k2   F] 
ko 


Since     r     has  only  a     <})     component  which  is  independent  of  the  angle     <{> 


V  •    F     =     0 


and 


iu*     br  Y  ikR 

H<t> =  iw  F<t> =  ~  tt  J  Ep  p'dp'  J   cos  ♦'  Sr dt,>'  (6J) 

a  0 


Equation   (6.1)  is  exact — no  approximations  have  been  imposed.     If    E   ', 
the  electric  field  intensity  in  the  aperture,  were  known  exactly,   then 
the  integration  could  be  carried  out  and  the  result  would  be  exact. 

If  we  approximate  the  aperture  field  as  being  that  of  the  inci- 
dent field  in  the  aperture,  i.e., 

E'(p,0,u)   =  -5H- 
p'   £n  | 


then 


«/    i   br  ?* 


*"  i  i  A 


In  the  far  zone  the  distance  R  becomes 


R  =  (p2  +  p'2  -  2PP'  cos  <J>'  +  (z-z')2)1/2 


=   (r2+  p'2-  2p'r  sin  6  cos  $'   -  2rz'cos  6  +  z'    )1/2 

=  r[l  +  2^"  -  %~  sin  6  cos  <f>'  -^  cos  6  +  ^y]1/2 
r  r 
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r[l   -  £-  sin  0  cos  (J)1] 


since     r  »  p1     and     r  >>  z'   , 


ikR  ei'<(r  "  P ' s i n  0  cos  cj)1) 


We  have  retained  more  terms  of  the  approximation  to     R     in  the  phase 
factor  than  in  the  amplitude  factor,  since  small   variations  in    R 
cause     large  changes  in  the  value  of  the  exponent.  The  amplitude 

of  the  function  is   relatively  insensitive  to  variations  in     R     as 
r  ■+  <»    whereas  the  phase  is  not. 

Upon  incorporation  of  the  well  known  expansion  [28,  p. 230] 


n=_oo 


eikp'sinecos  +•   =    j  Mji  jn(kp,sin  e)ein*' 
equation   (6.2)  becomes 


.,      b  2tt 

H*  ■  -  "*r  ^^f—        cos  ♦'  I   (-i)  Jn(kp'sin  e)e      d*'dpl 

*  (£n  -)r    Ja  ^  n=-~ 

..      b  2^ 

=  _       o  M  e 2  cos^.^.   j  (kp'sin  e)dp' 

ZlT  £n  2.      r  1 

*n  a  a  0 

(jjeoV(co)     eikr 


*na 


B—    H_—        ^(kp'sin  e)dp' 


u£QV(co)  eikr  "?  d(Jo(kp' sin  9)) 


b         r  k  sin  6 


*na  a 


to    V(«o)     e^      J(kb  sin  6)  -  JQ(ka  sin  9) 
%    Inl         r  Sin  6 
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The  approximate  time  behavior  of  the  far  zone  field  follows 

directly  from  the  inverse  Fourier  transform  of  (6.3) 


H.(r,e,t)  =  -i-   fa    f-J 

*  ^W  %   Un  ?-)r  sin  9 


-i«(t-3 

e  V(w)[Jo(kb  sin  9) 


-  J  (ka  sin  e)]dw      (6.4) 

When  this  integration  is   carried  out  by  the  convolution  theorem,  the 
impulse  response  of  the  annular  slot  antenna  is   found  to  be  given  by 
equation   (5.10)  and  the  step  response  is   found  to  be  given  by  equation 
(5.11). 

From  this   result  it  is   apparent  that  in  the  annular  slot  antenna 
the  radiation  also  appears  to  emanate  from  the  edges  in  the  aperture 
plane.     The  radiation  starts  when  the  incident  current  on  the  conduct- 
ing surface  reaches  the  edge.     One  edge  is   located  at     p  =  b     and  the 
other  at     p  =  a   .     For  an  incident  impulse  of  voltage  the  outer  edge 
contributes  radiated  magnetic  intensity  in  the  amount 


f 


1 


0 


1 


I  < 


'I    V/(b  sin  e)2-T2 


x  <  -  —  sin 
—      c 


-sin  9  <  x  <  —  sin 
c  —      —  c 


\ 


\ 


b     . 
;       t  ^  -  sin 


J 

(6.4.1) 
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and  the  inner  edge  contributes   radiated  magnetic  intensity  in  the 
amount 

/ 


1 


0 
-1 


yo    irr  sin  6  iln  - 

a 


M  sin  Q)Z-  tZ 


t  <  -  -  sin 
—      c 


-  ;     -  4  sin  e<  x  <  -  sin  e  V 

it-  C  —        —   C  f 


\     x  >  —  sin 
—  c 


with 


t-t-t 


(6.4.2) 


We  make  the  assertion  at  this  point  that  the  tip  and  the  base  of 
a  cylindrical  monopole  antenna  (Figure  6-2)   give  off  radiation  in  a 
manner  identical   to  that  of  the  annular  slot  antenna.     In  this 
case,  the  radiation  from  the  tip  must  be  delayed  in  time  by  an  amount 
which  equals  the  time  it  takes  light  to  travel   the  length  of  the 
antenna,  h/c  .     We  further  assume  that  the  currents  on  the  conductors 
are  reflected  in  the  same  manner  as  the  currents  on  the  conductors  of 
the  open  ended  coaxial  waveguide.     That  is,  the  total   current  must  go 
to  zero  at  the  tip  and  the  base. 

In  Appendix  E  we  develop  the  radiated  fields  of  a  very  long 
hollow  cylindrical   antenna.     Transcribing  those  results,  we  have  that 
the  magnetic  intensity  of  such  an  antenna  is 


»>,»,.)  ..fr-SfaL-,.- 

v  yo  sin  9  lx\  \ 

a 


ikr    JQ(ka  sin  e) 


r 


sin  6 


This   result  is  identical   to  the  frequency  domain  radiation  from  the 
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Figure  6-2.     The  cylindrical   monopole  antenna 
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tip  of  the  center  conductor  of  the  annular  slot  antenna  determined  in 
(6.3)  above.     Since     ka  <<  1    ,  the  particular  form  of  the  end  does  not 
have  a  significant  effect  on  the  field  intensities.     We  therefore  take 
these  results  to  be  valid  for  all   thin  cylindrical   antennas  regardless  of 
the  end  configuration  (with  the  obvious  exception  of  those  with  end  loading.) 

From  the  above  discussion,  it  is  clear  that  the  contribution  to 
the  radiation  from  the  tip  of  the  center  conductor  is  identical   for 
h  =  0  and  for  extremely  large  values  of    h  .     We  assume  this  result  to  be 
valid  for  all  cylndrical  monopoles  with  n  >  0   .     Since  the  base  of  the 
antenna  is  excited  at     t  =  0   ,  when  the  incident  excitation  voltage 
reaches  the  ground  pi ane;  appropriate  mathematical   factors  must  be  in- 
corporated into  the  field     intensity  equations  to  account  for  the 
finite  time      h/c     required  for  the  excitation  current  to  reach  the 
tip  of  the  antenna,     (c    is  the  speed  of  light). 

Since  the  electric  field  at  the  junction  of  the  feedline  and 
the  antenna  must  be  continuous,  when  the  flowing  charges   reflected 
from  the  tip  reach  the  junction,   charges  of  opposite  sign  must  travel 
back  along  the  antenna  [8],     If  we  consider  that  the  annular  thickness 
(b-a)  «  A   ,  then  image  theory  is  applicable  and  it  is   readily  seen 
that  the  current  relaunched  onto  the  antenna  from  the  base  is  equal 
in  magnitude  and  opposite  in  sign  to  the  current  incident  from  the 
antenna. 

No  radiation  besides  the  initial   radiation  due  to  the  current 
incident  from  the  feedline  originates  from  the  base  of  the  antenna. 
The  current  reflected  from  the  antenna  tip  is  cancelled  by  the  image 
current  at  the     z  =  0     plane. 
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To  determine  the  total   contributions  to  the  far  zone  radiation, 
use  is  made  of  Figure  6-3.     Proper  accounting  must  be  made  of  the 
radiation  from  the  tip  of  the  antenna  and  that  reflected  from  the  ground 
plane.     It  is  readily  seen  from  this  figure  and  the  above  discussion 
that  the  radiated  magnetic  intensity  in  the  frequency  domain  is 


H.(r.e^.  /foJtfsLs^: 


ikr 

[J  (kb  sin  6)  -  J  (ka  sin  e) 


e      ° 


ikh(l  -  cos  0)         ikh(l  +  cos  6) 

x      {§ ^ }  3 

1  +  e2ikh 


kh  4  (2n+l)  j      ;         n=0,l,2,3,---  (6.6) 


The  case  for    kh  -  (2n+l)  *■    will  be  treated  later. 

The  time  behavior  follows  from  the  inverse  Fourier  transform  of 


(6.6) 


Mr.e.t)  =  —     f  H.(r,e,w)  e~ia)t  dw  (6.7) 

-00 

If  the  excitation  voltage  is  harmonic  time  dependent,  then 


V(w)=  Ref/Fir  A  6(a)-  u  )] 


From  (6.6)   and  (6.7),  we  have 
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Figure  6-3.       Effect  of  image  plane 
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H.(r,6,t)  =  f-S.    A 

♦  K     (£n  5r  sin  8 


Re 


oj  b 
J  (-?-  sin  e) 


o     c 


co  a                      c 
-  J0(-f  sin  Ql* —^ 


oj  h  oj  h 

i-Hl-  cos  9)       i  -Ml  +  cos  9) 


+  e 


3! 


1  +  e 


x      {e      o  c'} 


Upon  simplification  this  becomes 


H.(r,6,t) 


(e~  A  cos  oj  (t-f) 
0  0  c 

Mo  (Jin  ?-)r  sin  e 


oj  b 
Jo(-T-sin  0) 


wna  oj  h 

[J„(-?-sin  e)][cos(-2-cos  e)] 


ox   c 


COS 


OJQh 


oj  h 


b  i  a  ;     -2-?  (2n+l)  J  ;     n=0,l,2,3, 


(6.8) 


This  expression  shows  that  the  radiation  pattern  is 


F(e) 


ojb  ojh  oja  ojh 

Jj-4-sin  9)cos  -f--  [J  (4-  sin  e)][cos(-2-  cos  6)] 


0     c 


ox  c 


sin  e 


(6.9) 


The  radiation  pattern  given  by  (6.9)   can  be  compared  with  the  pattern 
derived  when  the  antenna  is   assumed  to  be  ^ry  thin  with  a  filamentary 
sinusoidal   current  distribution  [33].     For  that  case 
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co  h  ^ 

COS   — COS  (u)   -cos    0) 

F(8)  =  C     .     -       0C (6.10) 

sin  9 


In  the  limit  when  the  feedline  diameters  are  "thin"   (i.e.,  A  »  a,b)  . 
the    patterns         given  by  (6.9)  and  (6.10)   are  identical.     Equation  (6.8) 
shows  that    the  radiation  along     6  =  0     is  always  zero,  as  it  must  be. 

To  get  expressions  for  the  cases  where     kh  =  (2n+l)  j  , 
n=0,l,2,3,...,  we  must  consider  the  fact  that  ohmic  losses  are  present 
in  any  real   system  allowing  us  to  use  a  complex  value  of    k     or  we  can 
take     k     to  be  complex  for  the  sake  of  mathematical  expediency  and  let 
its  imaginary  part  vanish  at  a  later  stage  in  the  analysis.     Thus, 

to 
k  =  k1  +  ikp       ;  k2  «  k,   =  - 

With  complex    k     the  factor 

eikh(l  -  cos  9)        ikh(l  +  cos  6)  _  CQs[kh  Q] 

,  ~     2ikh  cos  kh 

cos[(k,+  ik2)h  cos  e] 
cos(k,+  ik2)h 

cos(k,h  cos  e)cosh(kph  cos  6)  -  i   sin(k,h  cos  e)sinh(k2h  cose) 
cos  k,h  cosh(k2h)  -  i   sin  k-.li  sinh(k2h) 

In  the  limit  as  k,h  ->  (2n+l)  j  ,  n=0,l  ,2,3,  •  •  • ,  this  expression  is 
seen  to  be  finite  for  kp  t  0  and  the  magnitude  is  only  limited  by 
the  fact  that  the  conductivity  of  the  material  is  not  infinite 
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lim  (cos(kh  cos  e)}  _^       cos[(2n+1)f  cos  9] 

kjh  +  Un+DJ       cos  kh  -i[sin(2n+l)J-]sinh(k2h) 

n=0, 1,2,3, •••     and    k2  =  e  >  0 

The  magnetic  intensity  becomes 


JE~      A  sin(cj  t  -  — ) 

H.(r.e.t)     -    /S.     ° Si 

♦  *  yo       (in  2-)r  sin  9 

a 

-JQ(o3o  |  sin  e)[cos((2n+l)|-  cos   9)] 


sin[(2n+l)J]  sinh(k2h) 

The  input  current  to  the  antenna  peaks  at  the  values  of     kh 
given  by     kh  =  (2n+1)5-  ,  n=0,l,2,3,"  •    .     The  explanation  for  this  is 
that  the  monopole  can  be  considered  as  a  simple  extension  of  the  feed- 
line.     The  tip  of  the  antenna  would  correspond  to  the  open  end  of  a 
transmission  line.     Any  odd  number    of     quarter  wave  lengths  from  the 
tip  along  the  line  toward  the  source  would  correspond  to  a  point  where 
the  impedance  is   a  minimum  (a  point  of  series  resonance)  and  the  cur- 
rent is   a  maximum.    When     kh  =  (2n+l)5-  ,  the  feed  point  is  an  odd 
number  of  quarter  wavelengths   from  the  open  end  of  the  extended  trans- 
mission  line. 

Taking  the  inverse  Fourier  transform  of  (6.6),  we  find  the 
radiated  magnetic  intensity  for  a  voltage  impulse  of  excitation 
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TT^ 5 ' 
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_     ( 
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;     elsewhere 
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;       0  <  n     <  -  sin  6 
-    n  -  c 


;       elsewhere 


JL  =  t-— --HI  -  cos  6 ) 

n  c      c  c 


I  J      0<Tn<-|sin 

— «-  —    n  —  c 


I  ■*  yif^r 


elsewhere 
^.t-r-fid  +  e-eJ-M     (6.12) 


For  a  unit  step  of  voltage  excitation  the  radiated  magnetic 
intensity  is  the  integral   of  the  preceding  expression. 
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Figure  6-6.      Radiated  magnetic  intensity  of  annular  slot  antenna 

n  b 

Ordinate:    H     =  H,(   ' 

Antenna  Dimensions:      h  =  0     and     a  =    .  8b 


r  Jin  -  sin   e)(volts) 
oo  a 
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The  magnetic  intensity  for  arbitrary  time  dependence  can  be  ob- 
tained by  convolution  of  the  input  excitation  voltage  with  the  impulse 
response  (6.12)   or  by  convolution  of  the  time  derivative  of  the  input 
voltage  with  the  step  response  (6.13).     As  expected,   the  radiated  field 
for  the  infinite  monopole  is  the  same  as   that  for  the  finite  length 
monopole  for  the  time  interval  before  the  current  pulse  reaches  the  tip 
of  the  antenna.     The  results,  equations   (6.12)   and  (6.13),   are  easily 
extended  to  the  infinite  cylindrical   antenna  by  letting     h  -*  °°  .     The 
radiated  fields  for  the  annular  slot  are  found  by  letting    h  -*•  0     in 
expression   (6.6).     It  is  also  worth  observing  that  the  results   (6.12) 
and  (6.13)   are  based  on  the  assumption  of  no  ohmic  losses  as  the  cur- 
rent pulse  travels  back  and  forth  between  base  and  tip  of  the  antenna. 
The  attenuation  can  be  found  using  the  conventional   transmission  line 
approach  for  small  attenuations  on  good  conductors.     The  determination 
of  the  time  domain  fields  in  the  case  of  a  lossy  conductor  is  com- 
plicated since  the  resistance  per  unit  length  is  proportional   to  the 
square   root  of  the  frequency.     The  incident  pulse  supplies  the  total 
of  the  radiated  energy,  ohmic  losses,   reflected  energy  and  energy 
stored  in  the   induction  fields.     The  stored  energy  of  the   induction 
fields  is  the  source  of  energy  for  the  second  and  subsequent  pulses 
launched  down  the  feedline  toward  the  source,  the  energy  radiated  after 
time     2h/c,  and  the  ohmic  losses  which  occur  after  time     2h/c  . 

The  theoretical  results  obtained  in  this  chapter  are  in  good 
agreement  with  the  experimental  and  numerical  work  of  Schmitt,  et  al 
[6],  Palciauskas  and  Beam  [10],  King  and  Harrison  [8],  Burrell  [34], 
and  Lamensdorf  [14]. 
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Figure  6-7  is   a  comparison  of  the  theoretical   results  obtained  in 

this  work  with  the  experimental   results  of  Schmitt,  et  al.  [6].     Reference 
[6]  states  that  the  second  pulse  measured  in  the  experiment  should  be 
moved  .74  ns  toward  the  origin  and  its   amplitude  should  be  about  3/2  as 
much  as  observed.     This   is  necessary  to  account  for  experimental   devia- 
tions  associated  with   locating  the  receiving  probe  a  finite  distance  from 
the  antenna.      In  Figure  6-7  these  corrections   are  applied  to  the  second 
and  subsequent  pulses. 

It  is  also  constructive  to  include  the  following  observation. 
When  a  cylindrical   dipole  antenna  is  infinitely  thin  and  center  fed,  it 
is  common  to  consider  it  as  an  open  ended  transmission  line  that  has 
been  spread  or  opened  out.     A  sinusoidal   antenna  current  distribution 
with  current  nodes  at  the  ends  is  assumed  and  the  radiated  fields  are 
obtained  on  that  basis.     If  it  is   further  assumed  that  the  current  at 
the  antenna  feed  point  is  equal   to  the  current  existing  at  this  point 
on  an  equivalent  length  of  open  ended  transmission   line,   then  the  field 
expressions  developed  from  the  assumed  sinusoidal   current  distribution 
are  identical   to  the  field  expressions  obtained  in  this  work  in  the 
limit  of  yery  thin  antennas. 

To  demonstrate  this  point,  we  first  observe  that  the  radiated 
magnetic  intensity  created  by  a  filamentary  sinusoidal   current  distribu- 
tion of     IQsin  k(h  -   |z|)     is   [33] 

il  eikr 
H*  =  2tt°  sin  6  [°os  kh  "  cos(kh  C0S  e)^ 
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The  driving  point  current  is     I0sin  kh   .     The  other  symbols  have  the 
usual  meanings   consistent  with  their  previous   use  within  this  chapter. 

It  is  also  common  knowledge  that  the  input  impedance  of  an 
open  ended  lossless  transmission  line  of  length     h     is 

Z  =  -i  Z  cot  kh 
o 

where     Z       is   the  characteristic  impedance  of  the  line.     If  we  drive 
such  a  line  by  a  perfect  voltage  source     V(u)),  the  source  current  is 
readily  obtained  by  Ohm's   law.     Since  the  antenna  driving  point  cur- 
rent is  assumed  to  be  equal   to  the  current  which  would  exist  at  the 
input  of  an  equivalent  length  of  open  ended  transmission  line,  we 
equate  the  driving  point  current  to  the  source  current.       Therefore, 


I  sin  kh  =       ■  /(")  ■■ 
o  -i Z  cot  kh 

and  ° 

I     =         V(")  , 
o       -iZ  cos  kh 
o 

By  substituting  this   value  of     I       into  the  field  intensity  equation, 
we  derive  the  general   result  for  wery  thin  antennas, 

u  [»  a  ,a         -VH  elkr  n       cos(kh  cos   6)  ■, 

Vr'6'w)  "  2^Z  rsin  9  t1 ZSTTh ]  (6.14) 


Outside  of  a  discrepancy  in  sign,  it  is  an  easy  matter  to  show  that 
this  general    result  reduces  to  equation  (6.6)  when  the  limiting  form 
for  small   arguments   replaces  the  Bessel    functions  of  equation   (6.6) 
and  the  characteristic  impedance  of  a  lossless  coaxial    line  is 
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substi  tuted,  i  .e. , 


o       2^JeQ  a 


The  discrepancy  in  sign  is  the  result  of  not  specifying  a 
voltage  reference  point  and  it  is  not  significant  relative  to  the 
results . 

The  approximate  transient  behavior  of  a  very  thin  dipole 
antenna  fed  by  a  transmission  line  with  characteristic  impedance 
Z      is   found  by  taking  the  inverse  Fourier  transform  of  equation 
(6.14). 
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Conclusions 


The   results  presented  in  this   report  provide  a  comprehensive 
analysis  of  the  open  ended  coaxial  transmission  line  or  coaxial  wave- 
guide antenna.     The  fields   radiated  from  such  structures  can  be  com- 
puted without  a  priori    assumption  of  the  aperture  field.     The  current 
reflection  coefficient  for  the  fundamental  mode  is   found  to  be 
-J  (ka)     for    ka  <<  1    . 

The  conventional   approximation  for  the  fields   radiated  from 
the  coaxial   aperture  antenna  gives  good  agreement  with  the  results  ob- 
tained through  the  Wiener-Hopf  method  for  thin   radii   in  the  feedline 
(kb.ka  «  1) . 

When  the  excitation  is  such  that     A  <<  b-a,    the  radiation 
pattern  obtained  by  the  conventional   analysis  method  is  shown  to  be 
valid  for  the  predominantly  forward  direction  only.     We  determine  that 
the  classical   radiation  pattern  obtained  by  conventional  methods  should 

I/O 

be  multiplied  by  the  factor      |cos   el''       for     A  -*■  0   .     The  resulting 
pattern  is  one  in  which  the  radiation  at  high  frequencies  is   confined 
to  the  forward  direction  as  expected  from  geometrical   optics  predic- 
tions. 

A  model   for  the  cylindrical  monopole  or  dipole  antenna  is 
developed  based  on  the  observation  that  radiation  appears  to  emanate 
only  from  the  base  and  tip  of  the  structure  in  pulse  excitation.     This 
model   is  used  to  determine  the   radiated  electromagnetic  field  of 
cylindrical  monopole  antennas   for  both  harmonic  time  dependent  and 
arbitrary  time  dependent  excitation  voltages. 
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The  mathematical   techniques   required  to  obtain   results   using 
this  model  are  extremely  simple  and  in  most  cases  hand  calculations 
are  sufficient.     The  results  obtained  show  good  agreement  with  experi- 
ment and  are  useful   for  both  analysis  of  the  radiation  field  when  the 
driving  voltage  is  given   and  for  synthesis  of  a  driving  voltage  to 
produce  a  specific  radiation  field. 

Future  research  can  be  aimed  at  removing  the  restriction  that 
the  width  of  the  annular  region  of  the  coaxial  waveguide  antenna  is 
very  much  smaller  than  the  wavelength  of  the  exciting  source.     The 
results  could  also  be  extended  to  the  thick  antenna  case  (i.e., 
kb.ka  >.  1).     Additional    research  effort  on  the  transient  behavior  of 
the  cylindrical  monopole  could  be  directed  toward  characterizing  the 
effect  produced  by  different  end  cap  configurations,  determination  of 
frequency  limitations,  accounting  for  dispersion  and  describing  the 
effect  of  higher  order  mode  components  of  the  aperture  field  on 
the  transient  radiation. 
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Appendix  A  -  Analytical   Determination  that  H  Modes  are  not  Present 
on  Structure 

In  this  appendix  we  will  show  that  H  modes  are  not  present  on 
the  coaxial  waveguide  antenna  excite-d  by  a  TEM  mode. 
It  is  obvious  that  symmetry  precludes  excitation  of  all  H  modes  except 
HQn   (n=l,2,3,---)  modes. 

By  the  edge  conditions   [19]     (see  Figure  2-1) 
we  know  that  the  longitudinal   current  must  approach  zero.     Thus  the 
sum  of  the  incident  current  and  the   reflected  (longitudinal)   current 
must  be  zero.     Therefore  the  reflected  current  on  each  conductor  has  a 
longitudinal   component  equal   and  opposite  to  the  incident  current.  All 
of  the  reflected  current  is  longitudinal  as  it  must  be  to  completely  can- 
cel    the  incident  current.     To  have     H        modes,  a  circumferential   com- 
ponent of  current  should  be  present  on  the  structure.     Since  no  current 

source  is  available  to  launch  a  circumferential   current,  the     H      modes 

on 

cannot  be  present. 

The  above  paragraph  notwithstanding,   let  us  assume  that  some 
current  is   launched  in  the  circumferential   direction  upon  reflection 
of  the  incident  current.     It  is  shown  in  Chapter  2  that  the  longitudi- 
nal  current  in  each  conductor  is  equal   and  opposite.     In  keeping  with 
this,  we  assume  that  the  reflected  currents  on  inner  and  outer  conduc- 
tors must  also  be  equal   in  magnitude  and  directed  in  opposite  circum- 
ferential  directions.     It  will  now  be  shown  that     HQn  =  o     (n=l,2,3,- •  •)• 
The  applicable  scalar  Helmholtz  equation  for  this  problem  is 
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[3       +   3       +13     +  k2       1    ]  E   (p.z.u.)     =     0  (A.l) 

9z^       dpd       p   dp  p         * 


E   (p,z,w)  =0  at       p  =  a,b   ;     z  <  0 


E  (p,z,w)  =  0  at       p  =  0       ;^°1Z£ 


with  edge  condition: 

E  (p,z,io)  -+  0  at      p  =  a,b   ;     z  -*  0 

Fourier  transforming  (A.l)  gives: 

^+ilp"^-^E#(p.a.»)-     0  (A.2) 

Using  the  regions  of  Figure  2-2,  with  H.     replaced  by     E,     and 

the  asymptotic  behavior  of    E  (p,z,aj)     as     z  •*■  ±°°  ,  we  find  that 

E     (p.a.oj),  E     (p,a,a))     and     E     (p,a,w)     are  analytic  in  the  strip 

shown  in  Figure  2-4.     Solving  equation   (A.2)  with     E     =  E       ,  E       ,  or 
E,     respectively,  yields 

E^     =  AK^yp)  (A. 3) 

E^     =  Bl^yp)   +  CK^yp)  (A.4) 

E^     =   DI^yp)  (A.5) 

(A. 3)  satisfies  the   radiation  condition  and  (A.5)   goes  to  zero  at 
p  =  0   .       At    p  =  a   , 
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E,    (a, a, to)  =     E     (a,a,u) 


E,    (a,a,aj)  =     E,    (a, a, to)   =  0 
(t>3  (j>2 


(A.6) 
(A.7) 


^pE^(p,a,u>)] 


p=a        K       r3 


p=a 


(A. 8) 


or 


E*  (a,a,o>)  =     Bl^ya)  +  CK^ya) 


C  = 


E^  (a,afu>)  -  Bd^ya)) 


(A. 9) 


^CpE^(p.o,u)] 


p=a 


(A. 10) 


From  (A. 5)  we  have 


^pE .  (p,o.u)] 


■  Ya 


p=a 


E^  (a.a.u)   I0(ya) 
TTTal 


(A.ll) 


where  use  was  made  of 


E     (a, a, to) 

^3 


(A. 11.1) 


Subtracting  (A.ll)   from  (A. 10)   and  using  (A. 8),  it  follows   that 
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^{pE-^p.a.u))] 


_  ^pEL(p'a»u^ 
p=a       K       T2 


p=a 


B  +  Ko(Ya)       lo{ya) 

^j  +  Y  aE^(a.a.U)[1qira-  +  y-^y] 


Define  the  left  hand  side  of  the  above  equation  as     G"(a) 


G"(a)  =  - 


E     (a, a, co) 


K^IY+  K^yaJI^Ya) 


(A. 12) 


At     p  =  b 


E.    (b,a,w)     =  E     (b,a,w) 


E"   (b.a.co)     =  E~  (b.a.to)     E     0 
<J>1  <J>2 


^CpE^p.a.a))] 


=  d^pEd>  (P»a'u^ 
p=b       ap       '2 


p=b 


or 


E+   (b.a.a.)     =  Bl^yb)   +  CK^yb) 
E+   (b,a,w)   =  Bl^yb) 


KjIybT 


(A. 13) 
(A.14) 
(A. 15) 


(A. 16) 


^{pE     (p.a.eo)] 


p=b 


=  Yb[BIo(Yb)   -  CKQ(Yb)] 


B  +  Ko(Yb) 


(A.17) 
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d^-pE<j>1(p'a*a,)] 


p=b 


+  Ko(^b) 

-Yb   E^(b,a,M)   K^y 


(A.18) 


where  we  have  used 


A    = 


E     (b,a,oj) 

Jn 


(A. 19) 


Subtracting  (A. 17)   from  (A. 18)   and  using  (A. 15)   it  follows   that 


g^tpE"  (p,a,u>)] 


p=b       dp       h 


p=b 


i^WI 


Define  the  left  hand  side  of  the  above  as     E"(a)   .     Then 


E"(a) 


-  B 


(A. 20) 


If  we  substitute   (A. 20)  into  (A. 12)  we  get 


G'(a)   -  E"(a) 


KJJZJ=  K^yaJI^ya) 


(A. 21) 


By  equating  the  values  of     C     given  in   (A.  16)   and  (A. 9)  ,   it  follows   that 


^(ya)     + 


Ii(Yb)       I,(Ya) 


E*3U*a'a))  =  K^BT  E^3(a'a'U))   "  B[k{w  ■  l<7(YaT]  Kl(Ya) 


Upon  use  of  equation  (A. 20)   the  above  relation  reduces  to 

+  Mya) 

El  ^'"'^   "ITlYFr^  U,a,a3)  =  E"(°t)[Ki(Ya)I1(Yb)-I1(Ya)K1(Yb)]. 


Substitution  of  this  equation  into  (A. 21)  gives   the  companion  equation: 
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Uyb)       E*/b'a'w) 
G   W   "  E  (a)  T^i7  =     ^(ybjl^ya)  <A'22> 

From  the  definition  of  G"(a)  we  note  that  it  is  proportional   to 

H     (a,a,(jj)  -  H     (a.a.w),  and    E"(a)     is  proportional   to    H     (b,a,u>) 

z3  z2  zl 

-  H     (b,a,oo)     with  the  same  proportionality  factors.     By  well   known 

z2 
boundary  conditions  we  know  the  surface  currents  on  the  outer  conduc- 
tors  are,   respectively: 


■apx[Hz  (a, a, to)  -  Hz  (a,a,w)]az  =  +1^  a^ 


a  x[H     (b.a.w)  -  H     (b,a,a))]a     =  -I     a 
P       ^1  ^2  *b 

In  this  problem  we  have  assumed  that    b-a  <<  ^x—    and  therefore  the 
parallel  plate  case  is  approached.     Under  these  conditions  with  TEM 
excitation,  the  currents  on  the  plates  can  be  taken  as  equal  in  mag- 
nitude and  opposite  in  direction,  with  a  great  degree  of  reliability 
Therefore     G"(a)  =  -E"(a)     and  equation  (A. 21)  becomes 

E     (a.a.oj) 
_^3 

I^YaKK^yb)   +  K^ya)] 


G"(a)   =  — (A. 23) 


Simultaneously  (A. 22)  becomes 

E     (b,a,w) 

G"(a)  =  ]- (A.24) 

K1(yb)[I1(Yb)  +  I^ya)] 

We  may  factor  the  denominator  of  the  right  hand  member  in  (A. 23)   and 
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(A.24),   using  the  techniques  of  Appendix  D. 

K1(Yb)[I1(yb)  +  J^ya)]  =  N+(ct)N"(a) 
with  N+(a)   =  N'(-a)  %  (|a|)"1/2         as     a  -  » 

with  | Im  a|    >  -kp 

Similarly 

I1(Ya)[K1(yb)  +  K^ya)]  =  M+(a)  M~(a) 

with  M+(a)  =  M"(-a)  %  (|a|)"1/2         as     a  -  » 

with  | Im  a|    >  -kp 

-1/2 
It  is  also  necessary  to  point  out  that     G~(a)  ^  a    '       as     a  ■»■  °°  . 

-1/2 
This  is  determinable  from  the  edge  condition  that     H  (p,z,u))  ^  z     ' 

as     z  ->•  0"     at    p  =  a,b   . 

Substitution  of  these  factored  forms  into  (A. 24)   and  (A. 23) 

and  algebraic  operation  yield  the  entire  or  integral   functions   (A. 25) 

and  (A. 26) 

E     (b,a,u)) 

G"(a)   M"(a)   =  — ^ (A. 25) 

M+(a) 

E     (a,a,u) 

G"(a)   N"(a)  =  -~z (A.26) 

N+(a) 

Each  side  of  (A. 25)  or  (A.26)  is  the  analytic  continuation  of  the 
other  side.     Both  sides  of  (A. 25)   and  (A.26)   are  analytic  in  the  strip 
-kp  <  t  <  kp   .     Since  each  side  of  (A. 25)   and  (A.26)   is  bounded  for 
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all    a,   they  must  also  be  constant  (Liouville  theorem).     Examination 
of  the  asymptotic  behavior  of  the  above  expressions  shows  that  each 
side  is  identically  zero.     E     (a,a,uj)  •>  0     as     |a|   ->°°  ,  since 

+  +  + 

E.(a,z,w)  -*■  0     as     z  ■»  0  .     Therefore     E.    (b,a,w)  =  E.    (a.a.u)  =  0   . 
Each  of  the  coefficients  A,B,C,D     can  be  expressed  by  a  linear  combina- 
tion of     E,    (b,a,w)     and  E,   (a,a,w).     This  conclusion  is  evident  by 
observation  of  equations  (A. 11.1),   (A. 19),  and  the  simultaneous  solu- 
tion of 


E+  (a.ct.o))  =    Bl^ya)  +  CK^ya) 
E+   (b,a,(0)    =   Bl^yb)   +   CK^yb) 


Hence     E,    ,  E.       and     E,     are  identically  zero  and     E,(p,z,w)     does  not 
<J>1       4*2  93  <J> 

appear  in  any  region  of  the  coaxial   antenna  structure. 
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Appendix  B.       Determination  of  Branch  Cuts   for    y 

In  this  appendix  we  will  develop  the  branch  cuts  in  the  com- 
plex a  plane  that  define  the  region  of  analyticity  of  Y  ,  a  multi 
valued  function  of    a 


Y 


-J777 


Our  starting  point  for  the  work  in  this  appendix  will   be  equations 
(2.9)   and  (2.9.1): 

H.    (p,a,w)  +0       as       p  -►  °°  (2.9.1) 

The  solution  to  equation   (2.9)  satisfying  (2.9.1)  was  found  in  Chapter 
2, 

H^   (p,a,u>)   =  A(a)   K^yp)  (2.15) 

The  asymptotic  value  of  the  modified  Bessel    function  [20]  is  known  to 
be 

MYp)     "   s/^Y?    e"YP     as     p  *  °° 


Hence 


GO 

HA  (p.z.eo)     «    — ■ f  A(a)   e"YP  e-iaZ  da 

*1  2Ap"      i 


2/yp 

as       p  ■*  °°  (B.l) 
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From  the  Sommerfeld  radiation  condition  we  know  that  the  field  inten- 
sity must  be  an  outward  traveling  or  evanescent  wave  for  large 


=y^ 


r  =  Vp  +  z 

Following  the  manner  of  Reference  [21,  p. 20],  we  note  that  if 
we  require  that 

Re  y  >  0  (B.2) 

Im  y  <  0  (B.3) 

then  the  radiation  condition  will  be  satisfied  by     HA  (p,z,w)    .     This 

♦l 
may  be  verified  by  examination  of  (B.l).     The  branch  cuts  will  be 

chosen  in  a  manner  which  ensures  that  the  conditions   (B.2)  and  (B.3) 

are  satisfied. 

The     a    plane  is   viewed  as  a  two-sheeted  Riemann  surface  with 
the  sheets  connected  along  the  branch  cut.     In  each  of  the  sheets  the 
function    y     is  single  valued.     The  sign  of  the  function  on  one  sheet 
is  the  negative  of  the  sign  on  the  other  sheet.     The  value  of    y 
becomes  discontinuous  only  if  a  branch  cut  in  the     a    plane  is  crossed. 

Let  us  define  the  branch  cut  of    y    such  that     Re  y  >  0     in 
the  top  sheet  and     Re  y  <  0     in  the  bottom  sheet.     Thus  the  two  sheets 
are  joined  together  by  the  curve  given  by     Re  y  =  0   ■ 

The  branch  point  is  obviously     k     and  we  will,   for  computa- 
tional purposes,  take     k     complex 

k  =  k,  +  ik2     ,         k,   »  k2  ,     k,,k2  >    0 
We  may  therefore  write 
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Y2  =  ex2  -  k2  =  (a  +  it)2  -  (k]  +  ik2)2 


=  (a2  -  x2)   -  (k2  -  k2)  +  2i(ax  -  k^)  (B.4) 


Note  that     Re  y  >  0     in  the  entire  top  sheet,   therefore  if  we 
specify    y2  =  re19     and    y  =  r1/2el9/2    with     0  <  9  <  2tt  ,   then 

Re  y  >  0     only  if     |e|   <  tt  .     Thus  the  branch  cut  must  be  given  by 

?  2 

j  ©  |   =  tt     or    Re  y    —  0     and     Im  y    =  0   .     The  branch  cuts  which  divide 

the     a    plane  can  be  easily  identified  by  examination  of  Figure  B-l 

and  are  seen  to  be  a  portion  of  the  curve     or  =  k,k2   . 

2  i'tt 

Since  the  branch  cut  is  given  by    y    =  re        it  follows  that 

on  the  branch  cut  y=r  '  e  *'     or    Re  y  =  0  .     The  graph   (B-2)   depicts 

the  branch  cuts  for  the  top  sheet. 

For  the  branch  cuts  in  the  bottom  sheet,  it  is  only  necessary 
to  use  the  opposite  sign  of    y     everywhere  in  Figure  B-2. 

It  is  also  worthwhile  to  state  that  for  this  branch  definition 

y  =     /r,  r2  e  where         8  =  * — 

n       x     .       3tt  tt  ^3tt       _.        tt  „  _  ..  tt 

2  —  *h  T"    *     ~  7  —  *2      T"  "  2"  5" 


n — 2~  n — 2" 

Y  =  -i  vk  -  a  or        </k  -  a    =  i'y 
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1      1  Re  y     and  Im  y     <  0 


Figure  B-l.      Regions   of  complex    ex    plane 
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Figure  B-2.      Branch   cut  for  multivalued  function  y 
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Appendix  C.       Useful  Expressions   and  Expansions 

In  the  main  body  of  this  work, extensive  use  is  made  of  ordinary 
and  modified  Bessel  functions  and  Hankel   functions.     Listed  below  are 
explicit  expressions,  asymptotic  approximations  and  limiting  forms  for 
some  of  these  functions  [20,24]    of  integer  order. 


(i)     Relations  between  modified  and  ordinary  Bessel    functions: 

-  -x-v  iri  ■*■  iri 

I   (z)  =  e  J  (z  ec       )       ;  (-it  <  arg  z  <f) 


v 
vtri      / ,  \       75-  rri 


Kv(z)  =  \  Tri  e2"  H^Cze2"      );         (-ir  <  arg  z  <\) 


H(1)(z)  =  J   (z)  +  i   Y  (z) 

V  V  V 


H^2)(z)  =  JV(2)-1YV(Z) 


(ii)     Wronskians: 


Iv(z)  Ku+1(2)  ♦  Iv+1(z)  yz)  -  I 
Wz)  Vz>  "  Jv<z>  Wz>  ■  h 


(iii)     Series: 


J(z)  =  1  -  * — 7  +  - — * - — y-  + 

0  (1!)*       (2!T  (3!T 
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lz2  (lz2)2 

Y   (z)  =  ^(in  §•)  +  C}  J  (z)  +  ^    ^ 5-  -  |    4      a      -h   ... 

0  ^         2  °  "     (l!)2       2    (2!)2 

(iv)     Limiting  forms   for  small   arguments   (when    v     is   fixed  and     z  -»■  0): 
Jv(z)  %   (lz)v/r(v+l)  ;         v  *   -1,-2,-3,--- 


Y   (z)  -v-  -iH(1)(z)  %  iH(2)(z)  ^  £  in  z 

O  O  0  7T 


Yv(z)  *  -iH,(1)(z)  %  iH^,2)(z)  %  -  ^-  (|)"v;     Re(v)   >  0 


Iv(z)  *  (lz)v/r(v+l)  ;         v  1  -1,-2,-3,.- 


K  (z)  *  -in  z 
o 


kv(z)  ^^-r(v)(|.)-v 


(v)     Asymptotic  Expansions   (when     v     is  fixed  and     |z|  •*■*>) 


Jv(z)  =  J^{cos(z  "  ?w  "?H  +  OC-plp)   ;     (|arg  z|    <  tt) 

Yv(z)  =7l7{sin(z  "  ?vit  "£)J  +  0(JZT)    ;     (|arg  z|    <  ff) 
I   (z)  vfiL-  ;  (|arg  z|   <  J) 


Kv{z)^Jke~z    '         ('arg  zl  <£> 
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(vi)     Derivatives: 

U_zvl  (z)]  =  zv  I     ,(z) 
dzL     vv  '  v-r 

(vii)  Recurrence  Relations: 


I     ,(z)  -  I  .,(2)  =  —  I  (z) 

v-r   '        v+r   '        z     vv   ' 


K     ,(z)   -  K  Al(z)  =   -  %■  K  (z) 

v-r  '       v+r  '  z    vv  ' 


(viii)  Other  Relations: 


I     (z)  =  I   (z) 

-V  V 


K     (z)  =  K  (z) 

-V  V 
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Appendix  D.       Factorization  Procedure 

The  purpose  of  this  appendix  is  to  furnish  a  factorization  pro- 
cedure for  the  Wiener-Hopf  kernels  which  arise  throughout  this  work. 
The  factorization  procedure   used  extensively  in  this  work   and  developed 
in  this  appendix  was  originated  by  Bates  and  Mittra  [23].     The  mathe- 
matical   development  of  this  procedure  is   first  presented  and  then  it  is 

applied  to  the  specific  expressions  to  be  factored. 

-Yh 
We  begin  this  appendix  with  a  factorization  of    e  with 


Y 


n — ? 

ya  -  k       using  the  well  known  theorems  of  [4], 


Theorem  A:     Let     f(a)     be  an  analytic  function  of     a  =  a+ix     regular 

-P 
in  the  strip    t_  <  x  <  x+  ,  such  that     f(a)  <  C|o|       ,  P  >  0     for 

|c|  -*■  «  ,  the  inequality  holding  uniformly  for  all     x     in  the  strip. 

Then     f(a)     can  be  decomposed  such  that 


f(a)   =   f+(a)  +  f"(a) 


with 

oo+ic 
f+(a)   =  i     f         ^-  dB  ;         x     <  c   <  T  <  T+ 


oo+id 

[ 

-<»+-id 


r<«>  -  -  sir  |    w^- de     ^     x.  <  x  <  d  <  T+ 


where  f  (a)  is  analytic  in  the  portion  of  the  a  plane  defined  by 
x  >  x  and  f~(a)  is  analytic  in  that  portion  of  the  a  plane  de- 
fined by     x  <  x+   . 
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Theorem  B:       If    Jin  g(a)     satisfies  the  conditions  of  Theorem  A,  im- 
plying that     g(a)     is  an  analytic  function  of    a  ,  which  is   regular 
and  nonzero  in  the  strip     x_  <  x  <  x+     and    g(a)  ♦  1     uniformly  as 
a  ♦  ±°°     in  the  strip,  then     g(a)     can  be  factored  such  that 


g(a)  =  g  (a)   g"(a) 


with 


«H-1C 


g+(a)  =  exp     ^    j       ^M  d6     .       T_  <    c  <  T  <  T+ 


-c°+1C 

aH-id 

9a 
B- 


g"(a)  =  exp     -  i  &n  9(g)  dB  ;  x     <  x  <  d  <  xx 


-aH-id 

The  function     g  (a)     is  analytic  in  that  portion  of  the     a    plane 
defined  by     x  >  x     and     g"(a)     is  analytic  in  that  portion  of  the 
plane  defined  by     x  <  x+   .     The  principal  branch  of  the  logarithmic 
function  is   used  throughout.     It  is  also  worthy  to  note  that  if 
g(a)  ^  t— r      and     in  g(a)  ^  -£n|a|      as     |a|   ->-°°  ,   Theorem  B  can  still 
be  applied.     The  integrals  are  convergent  in  the  sense 

T+ix 


11m  f     {     }  dB 

T  ■*  °°       J  . 


-T+i  x 
If    f(a)  and  g(a)     are  even  functions  of    a     then 

f+(-a)   =   f"(a)         or         f+(a)   =   f"(-o) 
and 
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9  (-a)  =  g"(a)  or         g  (a)  =  g"(-a) 

The  function     e~Y     ,  which  we  wish  to  factor  is  interpreted  as 


fll/E"?    ^    a2>k2     ^    e+ih7^T7    fQr    k2>a2 
£n  g(a)  =  -yh  =  -h(a2-  k2)(a2-  k2) 


Examine 


The  first  factor  is  an  entire  or  integral   function  of    a  .  Therefore 
we  only  need  to  decompose  the  second  factor     1/y   .     Since     y     is  mul- 
tivalued, we  will   use  the  branch  cuts  for  the  top  sheet  of  the  two 
sheeted  Riemann  surface  of  the     a    plane   (see  Appendix  B). 
From  Theorem  A,  we  have 

°°+ic 
?M  =  ?~T     I  *  I       -k2<c<t  (0.1) 

-oo+ic      7e  -  k   (6-a) 

The  integration  path  for  equation  (D.l)  is  diagrammed  in  Figure  D-l. 
The  integral   (D.l)  will  be  evaluated  by  integration  along  a  contour  in 
the  upper  half  plane,  which  does  not  cross  the  branch  cut  located  in 
that  half  plane  (see  Figure  0-2).     By  Cauchy's  theorem  the  value  of 
the  integral   (D.l)   is 

f+(a)  =  — l—  ~  ^T     (      * (°-2) 

r> — 7    ""    J       r? — 7 

Ja-  kc  P}+  P2  V|3  -  k^(B-a) 

The  contributions  to  the  integral   from  the  semicircular  portions  of 
the  path  at     IB]-*00      are  zero  by  Jordan's   lemma. 
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Im#l 


ik< 


(3  plane 


\ 


Re/3 


^x 


-IC 


•  a 


\"ik2 


k  ■  k]  +  ik2 


Figure     D-l 
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Figure   D-2 . 


iIm/3, 


Re/3 


Figure   0-3. 
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For  the  limiting  case  of  no  losses,  k«  *  0   ,   the  path     P,+  P? 
is  adjacent  to  the  coordinates  as  shown  in  Figure  D-3.     Under  these 
circumstances  the  remaining  integral   of  (D.2)   is  expressed  as  the 
five  separate  integrals  which  follow: 


0 


EFT 


idy 


dx 


i  Jy2+  k2  (iy-a)  0  i  v/k2-  x2  (x-a) 


+    lim 


-y  irei9de 


r  *  °    i    y/2rkeie+r2ez19(k  +  re1'9  -  a) 


dx 


-iJkZ-  x2   (x-a) 


00 

f 


idy 


0     -1  v/yZ+  kZ  (iy-a) 


By  reversing  the  limits  on  the   last  two  integrals  and  evaluation  of 
the  third,  we  can  reduce  the  above  set  of  integrals  to: 


dx 


+    1 


_dy_ 


0   (x-a)  7k2-  x2  0      v/y2+  k2  (y+ia) 


A  further  reduction  is  achieved  by  a  change  of  variable.     Let 

x  =  k  sin  6     and    y  =  ik  sin  6   .     Incorporation  of  this  change  yields: 


tt/2 
1     ' 

IT 

0 


d6 


k  sin  6  -  a 


-loo 


d6 


k  sin  6  -a 


Note  that  use  has  been  made  of  the  identity,       sin(-iy)  =  -i  sinh  y 
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From  Reference  [20,  p. 78]  we  have  the  following  integral 


a+b  sin  z 


(  2  a  tan  \  +  b  2         2 

— * — 3-T77"  arctan[ — 5 — ?  i/?  3  5    (a     >  b  ) 

(aZ-bZ)1/Z  (aZ-b2)1/Z 

1  a  tan  |  +  b  -(bZ-a2)1/2        ,       - 

i?^?7?"n[atan^b+(bJ-a2)'^;<b    "  ' 


Through  this  integral  we  are  able  to  complete  the  evaluation  as   fol- 


lows: 


r-  atan  j-+k" 
2     .2.1/2       arctan      2      *  ,/2 


u(ac-  k*) 


r) 


arctan 


J0 


6  -ic 

'-  atan  ^  +  k 


r-a  tan  |  +  k  -  (k2-a2)^2 
£n  5 * — s-r  in 


Tr(k2-  a2)1/2        [-atan  £+k+(k2-a2)\ 


^2     .2*1/2 

(a  -  k    ) 

(a2  >  k2) 


■a  tan  |+ k -(k2-^)31"1 


■  Tool 


V 


atan  |-+  k+(kZ-c^)^_ 


(k2  >  a2) 


v 


Since     tan(-iy)  =  -i   tanhy       and       lim    tanh  y  =  1    ,   the  above  relation 

y  ■*■  °° 
simplifies  to: 


( 


7r(a2-  k2) 


1 
ir(k  -a  ) 


w  htan(7^>]     ;   (°2  >  k2) 

L  va  -  k 


.2     2.3s 


2       2x1/2, 


an 


r+ia+  k  -  (kfc-a  P  1         T-a+  k  -  (kc-  a  ) ,/t" 
[+ia+k    +(k2-a2)^J      'L-a+ k  +  (k2-a2)1/2  . 


(k2  >  a2) 
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The  trigonometric  identity 


z1-z2 
arctan  z,  +  arctan  z-  =  arctan[ ] 

t  r  2,z2 


was  applied  to  the  top  expression. 

Recognizing  the  relation  [25,  p. 62] 


tan^z-^lnC   j   ~+  }zz  J 


and  applying  it,  we  are  able  to  achieve  great  simplification 


rY  -   (a+k)-,       .     f  2  N  .  2x 
K  *  (a+kj]       '     (a    >  k  } 


iry 


_L    ,n[Y-   K^j]        ;      (k2  >  a2) 

iry         Ly  +   (a  +  k)J  ' 


2         2  2         2 

Note  that  the  same  result  is  found  for     a     >  k       and  k     >  a     .  When 

the  above  is   reinserted  into  (D.2),   the  decomposition  is   completed 

and 


f+(a)  - 1  [i  + 1  top  :  iat  Etaj 


-  w  -  i-it^SHfn 


iry 


J_   to[ia!j;i  +  Y]  =  i        [Y+   (a-kjj 
Try  (a+  k)   -  y         ""T  y  -   [a-  k)-1 


The  factorization  of     g(a)  =  e"Y       can  now  be  easily  completed.   From 
the  above  we  have 


-147- 

«■  ^  ■  ♦  ¥  «»^Mf^ 

Because     g(a)     is  even  in     a  ,  it  follows  that 
in  g  (-a)  =  Jin  g"(a) 


and 


Finally 


«••"<•>  ■♦*£  '"^4f^ 


+  iin 

9+(a)  =  g-(-a)  ■  e      *      ^  t  fe'-k)3  (D'3) 


This   result  is  consistent  with  that  of  Noble  [4,  pp. 20-21],   and  is 
shown  to  be  correct  by  the  following: 


g(a)  =  g  (a)  g"(a)  =  e 


_  %-qftSI3  *  -t^})]} 


=  e       5  =     e"hY  (D.3.1) 


where         in  z  =    in\z\  +  i<j> 

The  only  possible  singularity  of    g   (a)     in  the  upper  half  plane  is  a 
branch  point  singularity.     It  is  obvious  that  this  singularity  is  not 
present  and  that    g  (a)     is   regular  in  the  upper  half  plane,  since 
substitution  of    ±i|y|      for    y     has  the  same   result.     Similarly,   for 

g"(a)  , 
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ffcfc- 


:ih  M{£n  e  i4>]        -h<j> 


=  e 


<f>  =  2  tan 


-1    lY 


a- 


Theorem  C:    [23]     Let    G(a)     be  an  analytic  function  of    a    in  the 
region       x_  <  a  <  t+    with  the  following  conditions: 

1.  G(a)     is   regular  in  the  region     x_  <  a  <  x+  . 

2.  G(a)  1  0   ;     6(a)  =  G(-a)   (even). 

3.  G(a)  ^  cav  e"   'a'      as     |  Re  ct|     ■*  °°      where     v     and     h     are 
real  constants. 

Then     G(a)     can  be  represented  within  the  strip  by 

G(o)  =  G+(a)  G'(a) 

where     G  (a)     is  analytic  and  nonzero  in  the  half  plane     t  >  t_     and 
G"(a)     is  analytic  and  nonzero  in  the  half  plane     t  <  x+  .     The 
expressions  for    G  (a)     and     G"(a)   are  given  by 

G+(a)   =  G-(-a)  =   {v^GW  (1+  ^)W2} 

{B+(a)     exp[^l£n[l^}]-^]}  (D.4) 


with 
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(a) 


B  (a)  =  exp 


°°-id 
1  f       aF(g) 

2rT      J        B(B-a) 
-oo_i  d 


de 


;       |  Im  otl  <|d|;    It  I    >   Idl 


(D.4.1) 


(b) 


F(B)  =  h(62-  k2)1/2  -  v  £n(B2-  k2)1/2  -  in  c  +  in  G(B)  (D.4.2) 


(c) 


2       2   1/2 
Y  =   (a  -  k  )  '       with  branch  of  Appendix  B, 


(d) 


k  =  k,+  ik«     ;     k,   >>  k?   ;       k«,k,   >  0     and     |t_|   <  k~   ;     t+  <_  k~ 


We  will  now  show  the  validity  of  the  above  theorem.     We  begin  by  observ- 
ing that 


'-id 


B  W=expU      j        ^-^    d6 

L  -oo-i(J  J 


If    B     is   replaced  by     -B     and  use  is  made  of  the  fact  that 
F(B)  =  F(-B)     then     B+(-a)     becomes: 


B  (-a)  =  exp 


and 


[-oo+id  -i 

1        f       aF(B)dB 
2iri     J         B(B-a) 
«H-id  J 


B+(-a)   =  B"(a) 


We  next  find  the  product 
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B+(ct)  B"(a)  =  exp 


1 
2~FT 


r  aF(B)   dB 
BlB-a) 


where     c     is   the  contour  shown  in  Figure  D-4.     Observe  that  the  asymp- 
totic behavior  of  the  integrand  can  be  used  to  show  that  there  is  no 
contribution  to  the  integral    from  the  vertical  segments  of  the  contour 
at     |  Re  8 1  =>  °°     .     By   Cauchy's  theorem,   the  value  of  the  integral 
exponent  in  the  above  expression  is 

2iri  I    (residues  of  the  integrand,  from  the  poles  at 
3  =  0     and     8  =  a  ) 


=     -F(0)  +  F(a) 
=     hy  +  ihk  +  In 


I 


G(a)   , 

2  v/2 
G(0)(l-\)       • 


Therefore 


B+(a)   B"(a)   =  G(ct>  U&X+  ihkJ 

2  Vfn 


c  v/o 
G(0)[1  •-  2_] 
IT 


From  (D.4),  we  formulate 


G+(a)   G"(a)   =  G(0)(1   -  **■) '**   B+(a)    B"(a) 


7' 


«p|¥  -^W  +  ^^5$  -  i»] 


Substituting  for     B  (a)   B"(a)     and  using  equation   (D.3.1),  we  find 
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Im/3 


-00 


°°— *a  Re/3 


>a 


v. 


^TTd- 


Figure   D-4. 
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G  (a)   6"(a)  =  G(a)    ,  which  establishes  Theorem  C. 

Precise  knowledge  of  the  asymptotic  behavior  of    G  (a)  =  G"(-a) 
as     |a|   ■*  °°    within     t  >  x_     is  necessary  for  successful   use  of  the 
Wiener-Hopf  method.     As  our  next  task,  we  will   determine  this  asymptotic 
behavior.     Notice  that     B  (a)     is  bounded  as     |ct|   -»- °°  ,  therefore  by 
letting     |a|  ■*  °°    in  (D.4),  we  find  that 


G+(a)   =  G-(-a)     *     (|a|)v/2     1 


ilh. 


An 


m    e 

0(1 -Ma 


L±MO 

Y-(a-k) 


i   ^.£n|^| 
=   (|a| )         e 


(D.5) 


for  |a|  ■*• «°  within  t  >  t_.  In  the  manner  of  [23]  the  factoriza- 
tion formula  of  Theorem  C  will  now  be  converted  to  a  form  convenient 
for  numerical  work.     In  this  development,  we  concentrate  on  changing 


B   (a)  =  exp 


+°°-id 
1       f      aF(g) 

B(B-a) 
-oo_id 


dB 


to  a  more  tractable  expression. 

Taking  the  logarithm  and  integrating  by  parts  yields 
°°-id 


in  B+(c)   ■  jJt-    {       F(B)    #}  f 


-oo_id 
ao-id 


=  F(B)  An(^) 


■id 


A(B)  in(£j2L)   dB 


where 
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A(g)  =  gh  ,  Bv  ,  1         dG(p) 

2Tr(k2-62)1/2         2ri(k2-B2)       2-rri   G(B)     dB 


If  we  let     0  ■»■  -B   ,   then 


£n  B+(a)   =       |       A(-B)  in  (§±2.)   dB 


This   integral  will   be  evaluated  on  a  closed  contour  in  the  upper  half 
plane  by  the  residue  theorem.     Therefore 

Zn  B+(a)  =  -   I    A(-B)  £n(l  +  |)   dB  -  2tti   [(residues  at  the  poles 
i  of  integrand  inside 

contour  ) 

The  path     p     and  the  contour  are  given  in  Figure  D-5.     The  contribu- 
tion from  the  semicircular  path  at     |B|   ■+ °°     is  zero. 
We  now  assume  that 

(1)  G(ct)  has  at  most  one  branch  singularity  in  the  upper  half 

2       2  1/2 
plane   located  at     a  =  +k     in   the  form    y  =  (a  -  k   )  '    . 

(2)  G(a)  has  only  simple  poles  in  the  upper  half  of  the     a 
plane,   located  at     a  =  +p      with     | Im  p   |     >  t+   ,  n=l,2,3,-.- 

(3)  G(a)  has  only  simple  zeros   in  the  upper  half  of  the     a 
plane  located  at    a  =  z      with     | Im  z J    >  x+  ,  m-1,2,3,--*   . 

Examination  of  A(b)  shows  that  its  only  possible  poles  within 
the  contour  are  those  of  g-a  ^iq\  A  •  The  poles  of  this  function 
occur  at  the  poles  and  the  zeros  of    G(B)    .     The  residue  contribution 
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Figure   D-5. 
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Figure   D-6, 
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from  the  poles  of    G(3)     is     (p~r)     and  from  the  zeros  of    G(3)     it 
is   (^-y).     Since     G(3)     is  an  even  function,  when     3     is   replaced  by 
(-3)  the  contributions   from  the  poles  and  zeros  become  the  negative 
of  the  above.     We  may  now  write 

in  B+(a)  =  -   f  A(-B)  in(l  +  |)d3  -  ]      in[l +£-]  +  ][     4n[1  +  2_] 
I  B  n=l  pn     m=l  zm 


Set 


u=(kV)1/2  =  KB2-k2)1/2 


where  the  branch  on  the  upper  Riemann  sheet  is  used  (see  Appendix  B) 
B.,{U2.  k2)V2=(k2.   u2,l/2 


Using  this  substitution,  we  have 


Jin  B  (a) 


A(-B)  an(l  +  #)dB  =  C(u)  *n[l *a  *  i^]du 

6  J,  (k2-uZ)1/Z 


where 


C(u) 


1     d 


2tt 


2     2»l/2. 


♦      v  i     a     j.       G((k._/)i/.)] 

2tti  u       Ztti    du 


(D.6) 


,2     ,2. 


(see  Figure  D-6). 

From  Appendix  B  we  know  that  on  the  branch  cut     (&*"-  kc)  =  0   , 

2       2 
therefore  (3  -  k  )     is  pure  imaginary  and  the  variable     u     is  pure 

real.     So  the  introduction  of  the  variable  change  maps  the  contour    p 
into  a  contour  along  the  real   axis.     The  contour     p'     must  be  in- 
dented around  any  poles  of  the  integrand  which   appear     on  the  real 
axis.     The  indentations  are  consistent  with  pole  movement  from  the 
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real   axis  if  the  medium  is   assumed  slightly  lossy.     Note  also  that  we 
have  accounted  for  a  negative  sign  by  taking     p'     in  the  opposite 
direction  to  the  path  which     p    maps  into  under  the  variable  change. 

Examining     C(u),  we  note  that  the  second  term  has  a  pole  at 
the  origin,  and  the  third  term 

1     d     r.     pm.2     2>l/2n  1 d    r,,,2     2x1/2, 

"  77T  HIT  Un  G((k  -u  )  '    )]  = * — n   rw    ^r?  G((k  -u  )       ) 

Z7T1   du  (2ir1)  G((r-uZ)'^)   du 

could  possibly  have  a  pole  at  the  origin,   and  at  other  locations  along 
the  real   axis,  which  follows   from  the  fact  that     u     is   real.       The 
contribution  to     in  B   (a)     from  integrating  over  the  semicircular  path 
around  the  origin  is 


\  [R0-  v]  £n[l  ♦  f] 


where     R      is  the 
o 


lim     u  tiln  G(k2-u2)1/2)]  (D.7) 

u  -v  0       au 

1     N 
The  remaining  poles  will   contribute     +  j    I     Rn^   '  where 

n=l 

R„(a)  =     lim     (u-pj     iOn  G((k2-u2)1/2)]  Un[l  +        a       3} 


V     du 


u+p            n       au  [J- 

Hn  s/k  -  u 

where     p       is  the  nth  ordered  pole  of    C(u)  on  the  real  axis. 
rn 

We     may  calculate  the  integral   along  the  straight  line  portions 
of  the  path  as 
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f     K(u)     *n[l 


vk  -u 


where 


h         1 


ITT, 


and 


K(u)  =  £-3^  CUu)  +  L(ue,n)] 


L(u)  =  ^jOG((k2-u2)1/2)] 


(D.8) 


(D.9) 


The  integral  is  interpreted  as  a  principal   value  type,   and  it  is 
denoted  by  the  bar  on  the  integral  sign.     The  principal   value  of  the 
logarithmic  function  is  understood  to  have  been  used  throughout. 
In  summary,  we  have  found  that 

00 

in  B+(a)  ■     I       K(u)     4n[l  +         a      ]  du  +  i(R  -v)  £n[l  +  £] 
6*0-  Jtt  ° 

+  jl  Rn(a)   "     I     ^  +  H-]  +  £       *n[1  +  F"] 
n=l  Hn         m=l  m 


and  the  value  of    G  (a)     can  be  written  as 


R„/2 


6+(a)  =  J&(6)   (1  ♦  ft  °   '■  B       [1  +  ^-]"1   n     [1  +  2L] 


n=l 


'n        m=l 


m 


exp 


oo 

f 


-Bl+    +     K(u)  Hn[1  +         a        ]du 


6+0"1 


AT7 


(D.9.1) 
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If  we  perform  a  similar  derivation  in  which  we  use    B"(a) 

instead  of    B  (a),  as  the  starting  point,  the  result  would  be: 


In 


B"(oc)  =     {       K(u)  Jln[l  +  — z— ;]du  +  1{R  -v)     An[l  +  £l 

-  fl  Rn(a)  ♦  j    £n[l   -  JL]  -  j     inCl   -  §-] 
n=l  rn         m=l  m 


Similar  to  (D.9.1),  it  follows  that 


G"(o)  =  •GIUT  (1  +  g) 


R     /2         00  00 

*  o       n    [i  -^-]"1    n  (l  -?-) 


n=l 


'n        m=l 


m 


x    exp[^iln(T;  j^j)   -±£W 


6+0 


•J  k  -u 


1 


♦  £  I  V«) 


(D.9.2) 


From  (D.9.1)   and  (D.9.2),  and  knowing  that     G+(a)  =  G"(a)   ,  we 
recognize  that  the  factorization  is: 


G+(a)  =  G"(-a)  =     G(0)   (1  +  f)  °  '     H     (1  +  ^-)_1 


n=l 


x    H  (1  +  §~)  exp 

m=l  m 


^  "tHfcft1 "  ™ 


oo 

1      K(u)  An[l  +  — ^i—  ]du  +  i    I  Rn(|a|)l  (D.10) 

6->0  A  -  a 
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To  complete  the  theoretical  portion  of  this  appendix,  we  cite  a 
theorem  which  provides  for  the  expansion  of  a  class  of  functions  into 
infinite  products   [26, p.  136]. 

Theorem:     Let     f(a)     be  a  function  with  its  only  zeros  at     a-,  .a^.a.,,  ••  • 

where       lim    a      is  infinite  and  let     f(a)     be  analytic  for  all     a     . 

n  ♦  °° 
Then     f(a)     can  be  expressed 

f'(0)  /a 

awAv     °°  a/a,, 

f(o)  ■  f(0)  e     f(0)       n     {1  -^}e      n     ;         a     y  0  (D.ll) 

n=l  an  n 

If    f(a)     is  both  even  and  entire  with  its  only  zeros  occurring  at 

±a,  ,±a2>±a2>' ' '     and       lim    t  a    =  t  7  •   then     f(a)     can  be  written 

n  ■*  °° 


f'(0)  aa  aa 


f(0)      Z     M.oj  «nTr    Z     ri  *  fiU  Q"niT 


f(o)  =  f(0)  e     nuj       n     {1-|-}  enn    n     (1  +  |_)  e  niT       (D.12) 

n=l  an  n=l  an 


The  exponential   factors  for  the  infinite  products  ensure  convergence 
of  the  separate  products. 
Consider 

za  za  22 

/i       z  \  _     niT  .    M       za\  ^nir  „   ,      z  a      .       „-*.«» 
(1-— )e         ^(1-— )e      ^  1  -  -j-y    ;      n 

Thus, 

za 


r     z  a 
n=l  an  n=l  n  tt 


n    (i  .  |_)  e™   -v  1  -  J,  fzT     as     n  * 


r        1 

Because      J    "V    conver9es  as     n  *  °°  »  ^e  infinite  product  also  con- 
n=l  n 


verges  as     n  +  » 
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We  will  now  apply  the  foregoing  factorization  formulas  to  factor 

the  expressions: 

{I0(ya)   [K0(Yb)   -  KQ(Ya)]  }_1 
and 

{KQ(Yb)   [Io(Yb)  -  IQ(Ya)]  }_1 

Taking  the  top  expression  first,  we  define 

f(1)(a)     =     I0(Ya)     =     J0(a  Jk2-  a2) 

G(1)(a)     =     KQ(Yb)   -  KQ(Ya)  =   i   JCH^b  y/kT?) 


*{:h*J^n 


o 

where  use  has  been  made  of  the  well  known  relation  between  the  modified 
and  the  ordinary  Bessel   functions  [20, p. 370], 
It  is  readily  seen  from  (D.12)  that 

iaa    ^  .   iaa 

f(1)(a)  =  J  (ka)     n     (1  -  A)e"  ™      n  O+Ti^e    ™ 


o 


n=l 


W  n=l        'W 


where     a     =  ±i|a  |    ;     |a  |   =  y/(~)  -  k     ;     z       is  the  nth  ordered  zero 

of    J  (z)     and     z    ^  (n  -  xh     as     n  ■*  °°     .       The  value  of    ka    must 
o  n  4 

be  restricted  so  that     |a  |   1  0  ,  therefore  take     ka  «  l    . 

Theorem  C  is  used  to  factor    G*     (a).     From  the  asymptotic 
value  of  the  modified  Bessel  functions,  we  have 
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«<,,w^«-b|a|-/&«-a|a| 


.„-l/2e-a|a|     £[(1,1/2  ,+/».  (1)1/2] 


as     a  -*■  » 


Thus    v  =  -  ^-    and    h  =  a.     From  (D.9), 


L^Cu) 


■  IjOn  G((/^?)1/2)>  ^[In  i  JCH^(ub)  -H^(ua)] 


-b  H^^ub)  +  a  H^dia) 
H^(ub)  -  H*1}(ua) 


From  (D.7) 


R    ■     lim    u  -g—  L^   '(u)  ■  0     ,         L(u)  has  no  other  poles  on 

u  -*■  0  . .  ,  _ 

the  real   axis,  so     R     =0 

for  all     n   . 
By  (D.8)  and  the  analytic  continuation     H^2)(ze_Tri)  =  -eiv7ri  H(1)(z) 


*»>(„)  . 


a        ,     <-bHJl'(ub)  +  aNJl'(u«) 


-bHJ2)(ub)  +  aH|2)(ua) 
HoZ)(ub)  "  H0Z)(ua) 


a  _  1 


12-+2b[Y1(ub)Jo(ua)-J1(ub)Yo(ua)]  +  2a[Y0(ua)J0(ub)  -  J1(ua)Y0(ub)] 
[J0(ub)   -  J0(ua)]Z  +  [YQ(ub)  -  YQ(ua)]2 
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Factorization  formula  (D.10)  gives 


Thus 


G(1)+(a)  ={|i[H^(kb)   -  H^(ka)]}1/2  expC^  ln[£{£$.] 

CO 

.  ika    +    j    K(l)(u)     teCI*       l°l      1du3  (D.13) 

6^0+  A2-  u2 


+  -1 

X(o)  =  X  (a)  X"(a)  =  {l0(Ya)[Ko(Yb)   -  KQ(Ya)]   } 


X+(a)=  *— (D.14) 

0  n=l  MOfil 


ex     '     is  introduced  in  the  "plus"  function  and    e~x^a'     in  the  "minus1 
function  to  ensure  algebraic  behavior  of  each  of  these  functions  as 
|a|   -*-«>.     Note  also  that     x(-°0  =  -x(°0     since     x  (-<*)  =  X~(<*)    •     We 
determine  the  value  of    x(°0     from  the  asymptotic  behavior  of  the 
denominator  of  equation  (D.14). 

By  (D.5)  we  know  the  asymptotic  behavior  of    G*   '   (a)     is 


i  ^*n 


G^V)     *       |a|-1/4    e      * 


2a 


as       a 


in  the  upper  half  plane.     The  asymptotic  behavior  of  the  infinite 
product  in  the  denominator  of  (D.14)  can  be  found  from  the  well   known 
property  of  the  gamma  function  [4, p. 41], 
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n    [i  +     o   ]  e-a/an  =  e"Ca/a  — i 
_t  an+bJ  „/a  . 


r(T+  l) 


n=l 


Introduce  Sterling's  formula, 


1 


■<!♦!♦  i) 


r(o)  *  e"a  a       ^v^T  {1   +   (12a)"1   + 


}   as     a  ■*  °°     , 
|  arg  a|    <  tt 


Combining  these  two  relations  gives   the  useful   result 

~  P(r*D 


I    I1"  nHi  Cl  +  a^F]  e 
a  ■*■  n=l 


a    t  ,.  an       „   "  la 


JTn 


rl  b-, 

v  a  ' 


exp  (SJl  -  C  -  in  M-]}   (D.15) 
a  a 


In  (D.15)     C     is  Euler's  constant  (0.5772   .    .    .   ).     Applying  this   re- 
sult we  find  the  asymptotic  behavior  of  the  infinite  product  in  (D.  14) 
to  be 


eta 
ct     \_  niT 


1 
4~ 


,aa  I  .  1*i 


lim        n(1+      «     je"^     ,o|     Vexp{^Cl-C.£n|^|  +  ^]} 
a  +  °°     n=l  I  nv 

|  arg  a|  <tt 

where     |aj   '■*■  (n  -  j)  j  as     n  ■*  *  . 

By  collecting  the  above  results,  we  determine  that 


x(a)  -gfl-C 


m 


aa 


TT 


+  ^  +  m 


2a 


} 


(D.16) 


and  that     x+(a)  =  x~(-a)     behaves  as     |a|   '       as     |a|   -»-  °°     in  the 
region  of  the  a  plane  with     x  >  -k~   . 
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Turning  to  the  other  expression  to  be  factored, 

K0(Yb)[I0(Yb)-  Io(Ya)] 
we  define 

f(2)(a)  =  Io(Yb)   -  Io(Ya)  -  Jo(b  A2-a2)   -  JQ(aJ~*T7) 

G(2)(a)  =  Ko(Yb)     =     i  f-CH^^b/?!?)] 

Using  (D. 12),  we  have 

m  io/a+bx 

f(2)(a)  =  [JQ(kb)   -  JQ(ka)]     n  (1   -  jjl-j)e     m  "*" 

^  iot/  a+b  v 

x    n     (1  ♦  T^-r)eli^^~ 


m=l  i\%\ 


where     am  =  ±i  la  I      , 


-mx2       ,2 


ami    =  /(-C-)     -  k  and     z       is  the  mth  ordered  zero 

of    J0(z)  -  Jo(f  t) 


The  value  of    kb     is  restricted  to     kb  «  1    ,  so  that     |a  |   ^  0   . 

The  asymptotic  behavior  of  the  zeros  of  [Jn(z)   -  Jn(|"  2)^  are 
determinable  by  considering  the  asymptotic  behavior  of  the  Bessel 
function  for  large  arguments  and  considering  that  the  difference  be- 
tween the  values  of  a    and    b    does  not  affect  the  magnitude.     Differ- 
ences between  the  values  of    a     and    b  ,  however,   cannot  be  neglected 
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in  phase  factors.     Through  these  considerations  and  trigonometry,  we 
find  that  the  most  rapidly  varying  roots  are  those  given  by 


1         O     K 

\  *  (m  "  J]  Fb"       as      ■  *  -    5 


Since 


m  =  1,2,3," 


8(2l(o)^y^    e"blal         as     a -oo 
1 


then    v  =  -  j    and    h  =  b 


L(2)(u)  =  d_Un(i|)Hn)(Ub)]=     i 


H*"(ub) 


0        u  +  0       W 


(2) 


H^dib)       H<2)(ub) 


tt      2TnLu(T) 


Hj"(ubl    "  H<Z,(ub) 


-  b  n       2  i     _       i 


TT 


TTUb       T2 


j£(ub)  +  Y^(ub) 


G<2>+(«)  =  [Jl   H^(kb)]1/2  exp[l^£n[l^|]  -  J 

00 
+    -f     K(2)(u)     An[l  +  ■       la'      ]     du 


kb 


6-K)n 


vk  -  u 


(D.17) 


Collecting  the  results,  we  get 
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Y(o)  =  Y+(a)  Y"(a)  =  {  K  (Yb)[I   (yb)  -  Uya)]}"1 


where 


Y+(a)  = 


e#(o) 

~77\  /  »  ig(a+b) 

GUj+(g)  v/j   (kb)   -  J  (ka)     n     (1  +       °L_)  e     "^ 


m=l 


" 


(D.18) 


By  (D.5),  we  know  that 


icdb  0n|2a| 

c(2)+/   i  „    ,    ,-1/4    "  T"an|T' 
Gx       (a)  'V;  |a|  e 


as     a    ♦  °° 


in  the  upper  half  plane  and  from  the  asymptotic  behavior  of  infinite 

products   (D.15),  we  know 

00  i    /a+bv 

lim  n  n  a     \Jivir~2~~'       .    ,-1/4 rot(a+b)ri     P 

ii  ,       +  n rv6  %     a  exp^    )X — HI  -C 

|a|-*»          m=r  1latnl                                                i2tt 

|  arg  ct|  <tt  "* 


£n 


g(a+b) 
2tt 


*1») 


'%!  -  (m  "  }}  FF 


as     m  •*■  <» 


Therefore 


♦(a)  = 


g(a+b)  r. 


C  -  Jin 


g(a+b) 


ZF 


i-n-   .   2b     „    .2a 


♦t*^*"^        (D-19» 


and       Y+(a)   =  Y"(-a)  %   |a 


1/2 


as     |g|  ■*  °°    in     t  >  -k2  . 


-167- 
Appendix  E  -  Radiation   from  the  Infinite  Cylindrical  Antenna 

In  this  appendix  we  shall  determine  the  radiation  from  a  hollow 
infinite  cylindrical   antenna.     An  integral  equation   for  the  current  on 
such  an  antenna  will  be  formulated  and  solved.     The  fields  will   then 
be  determined  from  the  current.     The  integral  equation  considered  in 
this  appendix  was   first  formulated  and  solved  by  Hallln  [35,36].     He 
was  interested  in  the  reflected  current  on  a  cylindrical   antenna,  when 
the  incident  current  wave  comes  from  a  source  that  appears  to  be  an 
infinite  distance  from  the  end.     Although  the  integral   equation  is  the 
same  as  that  found  by  Hall^n,  the  application  used  in  this  appendix  is 
novel. 

Consider  Figure  E-l,  which  is  an  illustration  of  the  actual 
geometry  and  its  equivalent  representation,  which  is   amenable  to 
analysis.     Since  the  antenna  is  made  of  perfectly  conducting  material, 
its  surface  is  an  equi potential   and 


E  (r,t)  =  -v(t)  6(z)  when     p  =  a,  -«  <  z  < 


In  the  frequency  domain  this  becomes 


Ez(r,u>)  =  -v(u>)  6(z) 


In  view  of  the  rotational  symmetry,  it  can  easily  be  shown  that 

9^A  (r) 

E(r.co)  =  +io)[A_(r)  +  \ K-]  =  -v(u>)  6(z) 

z  z  vr      dz 


r    is  the  vector  from  the  origin  to  a  point  on  the  surface  of  the 
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Figure  E-l.     Infinite  cylindrical   antenna  and  equivalent  repre- 
sentation using  gap  generator  as  source 
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antenna  and    A  (r)     is   the  well  known  vector  magnetic  potential,  which 

satisfies  the  relation     B  =  V  *  A     . 

The  symmetry  of  the  structure  also  requires  that     A  (r)  ■  A  (z)> 

thus  the  axial   component  of  the  vector  magnetic  potential  satisfies  the 

differential  equation 


92A  (z)         -  .2 

S—  +  kZA,(z)  =  +  I£-  v(«)     6(z)  (E.l) 

3z 

We  also  have  from  symmetry  that    A  (z)  =  A  (-z). 

The  solution  to  (E.l)   obtained  in  a  straightforward  fashion  is 

,    /    >.   i  k  |  z  | 

A,(z)  =  C,cos  kz  +  C,sin  kz  +  *^W* 

Z  1  Z  Zu) 

From  the  symmetry  condition  that    A(z)  =  A(-z)     and  the  requirement 
for  only  outgoing  waves  on  the  infinite  structure,   C,  =  Cp  =  0     and 

i    /    \   i  k  I  z  I 

A7(z)   =  kv(")e      '  (E.2) 

z  Zco 

The  vector  potential  will  now  be  determined  by  an  alternative  fashion, 
whereby  use  is  made  of  the  free  space  Green's  function  and  the  surface 
current.     The  general   formula  is: 

*(r)-£    jK(r')£— - -    dS' 

with 

|r  -  r'|=y(z-z')2+  p2+p'2-  2pp'cos(4>-4>') 
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In  the  case  of  the  infinite  cylindrical  antenna  with  the  obser- 
vation point  on  the  cylindrical  surface,  this  general   formula  becomes 


Az(z) 


o 
4tt 


f  111!) 


2tt 


2tt 


ikr 


d<t> ■  dz ' 


(E.3) 


where 


r=   J(z-z')Z+  4a2sin2(^-) 


and      <j>        is  taken  as  zero.     Since  the  vector  potential  is  independent 
of    $  ,  the  particular  choice  of    <p     is  immaterial.     Since  (E.2)  and 
(E.3)  are  both  equivalent  expressions  for  the  vector  potential,  they 
may  be  set  equal  to  each  other.     The  results  are 


2tt 


27       dz'  Hz') 


ikr 


2-rrk  v(co)  e 


ik|z| 


uoo) 


2*f±      V(u,)  eiklzl 


(E.4) 


where     -°°  <  z  <  °°     . 

The  equation  (E.4)  is  valid  for  all  values  of    z  ,  since  it  is 
based  on  the  mathematical  model   of  the  actual   antenna  structure. 
However,  in  the  actual  antenna  we  observe  that  the  current  on  the 
feedline  after    t  =  0  ,  the  time  when  the  incident  excitation  reaches 
the     z  =  0     plane,  is  a  TEM  mode  wave.     The  current  wave  can  be  ex- 
pressed as: 

•ikz 


I(z)  =  1(0)  e' 


z  <  0 


-171- 

1(0)     is  the  magnitude  of  the  current  at  the  plane     z  =  0   .     In  view 
of  the  above,  it  is  justified  to  represent  the  antenna  feedline  system 
with  the  following  set  of  equations: 

/fo  -'   *-ikz 


27T    ikr  2tt    7T  v(o))e'^  ;       z  >  0 

1Kr  V  MQ 


,1    [  dz-  I(z')  j    V-d*'  = 

-«>  0  I  g(z)  j       z  <  0 


I(z)  =   1(0)  e'ikz  ;       z  <  0       (E.5) 


We  restrict  the  analysis  to  the  consideration  of  TEM  mode  feedline  cur- 
rent only.     The  function     g(z)    for     z  <  0     is  unknown.   The  solution  to  the 
integral  equation   (E.5)   is  the  current  on  the  infinite  cylindrical 
antenna  when  there  is  an  outward  (from  the  origin)   travelling  current 
wave  on  the  feedline. 

The  current  is  now  expressed  in  terms  of  its  spatial  Fourier 
inverse 


I(z')  =  -L     f  e"iaz'   1(a)  da 

-00 

1(a)  =   I+(a)  +  I_(a) 


wi  th 

oo 


1      f   ,#_..»    Jaz' 


T(a)  =  —      I(z')  e,u*     dz' 


I"(a)  =  —  f  I(z')  eiaZ'    dz' 


JTH 
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-k?z  k?z 

Note  that     |I(z')|<Afl  as     z  -  »     and     |  I(z')|    <.  B  e  ^         as 

z  -*■  -»    where     kp     is  the  imaginary  part  of  the  propagation  constant. 
1(a)     is  obviously  analytic  in  the  strip     -kp<  x  <  k„     of  the 
a  =  a  +  ix     plane.     I   (a)     is   analytic  in  the  half  plane     x  >  -kp     and 
I~(a)     is  analytic  in  the  half  plane     x  <  kp. 
From  (E.5)  it  follows  that 
0 


I'(o)  ■ 


1 


KO)  e 


-ikz  _+iaz 


dz 


mi 


fit  (a-k) 


for    x  <  k«   .     Additionally 

1 


r+/    \         If     2tt         m    »     ikz     iaz  , 
G  (a)  =  — — ;    v(o))  e        e        dz 


J 


6    /Mo 


"o     v(a>) 


for    x  >  -k. 


Upon  putting  these  results  back  into  (E.5),   we  obtain 


2tt 


2tt 


^2T   J 


■iaz 


1(a)   da 


.ikr 


d<fr'  =  -J- 


1 

2lf 


■iaz 


e_lc"  G(a)   da 


or  equivalently 


1 


dz' 


■iaz'rT+ 


Cr(a)  - 


11(0) 


2tt 


-00  — oo 


-]da 
m  (a-k)         j, 


ikr 


d<t>' 


■  j«-iaz[y?^{^*6-(.)]* 


(E.6) 
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2tt 


We  digress  at  this  point  to  find  an  equivalent  form  for 

,ikr 

'— —  d<j>'    .     It  is  well  known  that  the  free  space  scalar  Green  func- 

0  0 

tion  in  spherical   coordinates  satisfies   (V  +  k   )  G(f,r"' )  =   -6(r-~r'), 


and 


|r-r'| 


In  cylindrical   coordinates  the  Green  function  satisfies 


and 


(V2+  k2)   Gfp.p'.W.z-z')   =   -  6(P'P')   6(W)   fi(z-z') 


G(p,p',^',z-z')  =-U      I     e1"1^-*') 

4n    m=-«> 


r,-1a(W)    Im(YP<)   Km(YP>)   da 


2       2  1/2 
In  the  above  equations       y  =   (a  -  k   )   '      and     p<       is  the 

,  > 

(  ^fS+r  )   of    p     and     p'    .       These  two  Green's  functions   are  equiva- 

lent  three-dimensional    representations.     Therefore, 
ik  r-r' 


|r-r'| 


.  1     r    eim(<M>')     (  e-ia(z-z')  j  (       )K  (       } 

m=  -oo  ' 


da 


m=-«> 


and 


2tt    . 


ik|  r-r' 
r-r 


«•   =2      [e"1^2"2')    Io(YP<)   Ko(YP>) 


da 
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We  now  evaluate  the  above  equation  at     p  =  p'   =  a     and  sub- 
stitute it  into  (E.6). 


OO  00 


if" 

— OO  — OO 


e  LI  (a) 


il(0) 
*^T(a-k) 


-]  da 


-laz 


c  IS.  Mii  +  G-(a)]da 


Jo       (a+k) 


'  e-i6(z-z,)Io(Ya)Ko(Ya)d6 


(E.7) 


Interchanging  the  order  of  the  integrations   on   the  left  hand 
side  of  the  above  equation,   and  recognizing  that 


je-i(a-e)z'  dz,  =  2iT6(a_g) 


equation   (E.7)   transforms  into: 


00  OO 


'  da     f  e-i3z[I+(a)   -     "W     ]  I   (Ya)K  (Ya)  S(a-B) 


dB 


00 

r 


-co 


e-iaz  j. 


2tte 


M(4+G"(a)] 


u      Ta+k" 
o 


da 


Integration  over  the     3     variable  produces  the  result: 

00 

2     f  da  e"iaz  {[I+(a)   -  — liM_]  [I  (yt)  k  (Ya)] 

J                                         v^T(a-k)         °  ° 


r.  /%     v(a>)      .    G'(a)     ,    . 


(E.8) 
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Let     I(Ya)   K0(ya)  be  factored  into  the  product     L  (a)   L~(a), 
where  L  (a)     is  analytic  in  the  half  plane     t  >  -k,   and     L"(a)     is 
analytic  in  the  half  plane     x  <  kp. 

After  substituting  the   factored  form  into  (E.8)  we  derive  the 
following  result: 


2   f  da  e"iaZ  Ha)   {[I+(a)   -       UW   ]  L+(a)  - 

i  J7Zln-\t\ 


«£¥(a-k) 


--L-[i/5    -^-+G"(a)]}     =     0 
L'(a)        V    2Uo     (a+k) 

The  integrand  is  regular  in  the  strip     -k?  <  x  <  k?     of  the       plane, 
Since  it  must  vanish,  it  gives  the  result 

[iV)  .     1KQ)  ]  L+(o0  .  £M+  i  EL      v(M) 

v^F(a-k)  2L'(a)  V   2yo  (o+k)L"(a) 


2       2 
Multiplication  of  both  sides  by     a  -  k      yields   the  integral   or  entire 

function 


L+(a)   I  +  (a)(a2-k2)  -  2lMl^l  L+(a)   .  £W    (o2.k2) 

^  2L"(a) 

+  i    J]%     v(M)(a-k)  =     p(a)  (E9) 

2Po       L"(o) 

The  left  hand  side  is  analytic  in  the  upper  half  plane   (x  >  -kp) , 
the  right  hand  side  is   analytic  in  the  lower  half  plane  (x  <  kp)     and 
each  side  is   the  analytic  continuation  of  the  other.     We  divert  our 
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attention  from  (E.9)  for  the  moment  to  consider     L  (a),   L~(a)     and 

their  asymptotic  behavior  as     a  •*  «     in  their  respective  half  planes 

of  regularity. 

Earlier,  we  defined  the  functions 


L+(a)   L"(a)   =   IQ(Ya)   KQ(Ya) 


By  the  methods   of  Appendix  D,   this   factorization  can  be  carried  out 

with  the  results 

iaa 


G+(3)(a)   /TTkal    n     (1  +  TTTTr)  em7T 


L+(a)=L"(-a)   = 


m=l 


m 


>(a) 


G+(3>(a)  -  q  iH^ka)]1'2  expff  *n[l±{^}]  -  *£ 


(       K(2)(u)     £n[l   +       lal         ]  du 

j  +  n — 7 

6+0  Vk  -  u^ 


(u)  "  -  I  '  ■  ™£    ,2/...,   .  „2,- 


TT 


J^(ua)  +  Y^(ua) 


♦(a)  -£[1  -C-ln|2*J    +£n|^|+ij] 


m' 


z    o       2 
I3-.I    =   J  (— )  "  k  where     z_     is  the  mth  ordered  zero  of 


m 


J0(z) 
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+  -1/2 

L  (a)    ^    |a|  as      |a|   +«     in  half  plane     t  >  -kp 


L"(a)    ^    |a| 


1/2 


| ot|   ■*  °°     in  half  plane     x  <  k? 


as       a    ->■  °° 


Before  returning  to  (E.9),  we  note  from  the  edge  conditions 
[21]  that  the  current  goes  to  zero  at  the  terminus  of  the  hollow  con- 
ductor.    As     I   (a)     and     G~(a)   represent  transforms   of  current  we  can 
deduce  the  asymptotic  behavior  of  these   functions   for    a  -+  °°     by 
examining  the  behavior  of  the  respective  currents   as  we  approach   the 
ends  of  the  antenna.     By  considering  the  antenna  length  to  be  finite 
and  changing  the  variable  to  place  the  origin  at  the  terminal  ends, 
it  is  easily  shown  that     I   (a)     and    G"(a)   vanish  as     |a|   ■*  °°      in 
their  respective  half  planes  of  analyticity. 

Application  of  the  asymptotic  behavior  of  the  various   func- 
tions to  (E.9)  shows  that 


L+(a)    I+(a)    (a2-k2)    -  "WW    L+(a) 


<   |a|     as     a  -*  °°     , 
t  <  -k„ 


and 


2     ,2 


TTE 


G"(g)    Csf=J&   i    i  /    ""o  vM(a-k) 
L'(o)         ~ 


2y 


o     L"(o) 


<     |a|        as     a  -»■  °°     , 
x  <  k0 


From  the  extended  form  of  Liouvi  lie's   theorem,   P(a)     is  a  polynomial 
of  degree   less  than  or  equal   to  one  [4,   p. 38].     We  may  therefore  write 
P(a)  =  C,+  C2(a-k),  where     C,     and     C2     are  constants  to  be  determined. 
By  setting     a  =  k     in  the   left  side  of  equation   (E.9),  we  get 
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c  i2I(0)k  L*(k) 

1  '  ffi 


Additionally, 


C 


1         r     -   \+(    \   t  +  /    W     u\        il(0)   L"l"(g)(a+k) 
— TT  +  U  =  L  (a)   I   (a)(a*k)  +  — ^ ^— '- '- 

01  K         2  /2T  (a-k) 

By  letting     a  ■*  °°     in  the  above  expression,  we  find     C?  =  0   .     Thus 
P(a)     is  a  constant. 

From  (E.9)  and  the  value  of     P(a)     just  determined,  we  have 

I+(ct)  =  .        il(0)       .  i2I(0)k  lV) 

JR  (a-k)         •5¥(a2-k2)L+(a) 

Recall   that     1(a)  =   I   (a)  +  I~(a),  so  the  transform  of  the  current  is 


i(tt)  ,  .  «iWfc  i*M 

^7(a  -r)L+(a) 


i2I(°)k   L^k>    l'^ (E.10) 

^  (a  -k^)I0(Ya)K0(Ya) 


Inverse  Fourier  transformation  of  the  above  yields 


I(Z)   s  ■  lK0)k   L+(k)    7    e-iaZ  L-(a)   da (EJ1) 

i(a  "k   )I0(Ya)Ko(Ya) 

It  is  also  clear  that  if  we  close  the  contour  in  the  upper  half  plane, 

-ikz 

the  integration  produces     I(z)  =   l(0)e         ;   z  <  0   ,   as   required. 
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The  vector  magnetic  potential   can  now  be   found  from  the  general 


formula, 


A_  =  -£ 
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.00  —00 
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This  equation  is  valid  for     p  >  a   ,   z  >  0   .     Note  that  the  magnetic 
flux  density  is  given  by 
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9A 


therefore  the  magnetic  intensity  is 

il(0)k   L+(k)     J     e"iaZ  K^yp)    L"(a)   da 
H^p.z.iu)- -g 7T-^ (EJ2) 


Equation   (E.12)  is  exact  for  an  infinite   length  antenna.      If  we  assume 
that  it  also  holds   for  an  extremely  long  antenna  and  employ  the 
methods  of  steepest  descent  as  was   done  in  Chapter  4,   the   resulting 
equation   for  the  magnetic  intensity  is  given    by 


H 


t (n\~ikr        J  (ka  sin  6)         ,  +  ,.  x 
,     _     \       l(0)e                o                  '         L  (kj  /,-   ,,» 

(p,z,w)   =      i'  ttz   s v- L (E - 13) 


^,',u,=-j^ sTrTe 


L"(k  cos  e) 


Earlier  we  assumed  that     1(0)     was   the  incident  current  which  is 
reflected  from  the  discontinuity  at  the     z  =  0     plane  [8].     The  re- 
flected current  on  the  feedline   from  pulse  excitation  is  known  through 
experiment  to  approximate  a  reflected  incident  current  pulse.     From 
the  edge  conditions  for  this  structure  it  is  clear  that  the  current 
on   the  outer  conductor  vanishes  at     z  =  0   .     If  the  current  on  the 
outer  and  inner  conductors  are  to  be  equal   in  magnitude  and  opposite 
in  direction,  then  the  current  on  the  inner  conductor  must  also  vanish 
at     z  =  0   ,   the  point  of  discontinuity.     We  therefore  conclude  that  a 
pulse  of  current  is  launched  onto  the  antenna  and  an  oppositely 
traveling  pulse  is   launched  back  upon  the  center  conductor  of  the  feed 
line.     The  magnitude  of  this   current  is   assumed  to  equal   the  magnitude 
of  the  incident  pulse. 
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If  the  incident  field  is  given  by  equation   (2.7)  of  Chapter  2, 
then 

I(0>  ■  "2#  ;*4 

v     o     Jin  — 

3 

The  negative  sign  is   required  since  the  reflected  current  travels  in  a 
direction  opposite  to  that  of  the  incident  current.     Incorporation  of 
the  above   into  (E.13)   gives   the  final    form  of  the  magnetic  intensity, 

f^Z  ../    \  J^r    J   (ka  sin   8)         ,+/ln 
H   (p,z,W)  =   -     U^4  V-      °  sin    9  -^ (E-14) 

*  N/Mo  Jin  ?.    r  S1n    e  L  (k  cos   9) 


If 


lim 


l!ikl 


ka+O  L"(k  cos  9) 

then  (E.14)  obviously  reduces  to 

/^i     »/   \     Jkr    J~(ka  s"in  e) 

3 

for  an  antenna  with     kh  >>  1     and     ka  <<  1    . 

To  show  the  validity  of  this  limit  we  rewrite  it  from  the  ex- 
pressions for     L  (k)  and     L"(k  cos   8)    , 


ii (i+ w)emiT         G+<3)(k)  eH~k  cos  e) 

lim      — r—r r—r 

ka-0      "  M       k  cos  gvJ-ika  cos  9)/(rmr)  GH3)(-k  cos  9)  e^k) 
J1  Vl      TTr~1     ,e 


m=l 
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exp 


lim 
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lim      e  xp 
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-]  du 


This  final   limit  was  investigated  in  Chapter  5  where  it  was  shown  to  be 
unity. 
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